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Summary 


An  Imperfect,  laminated,  circular,  cylindrical,  thin  shell,  simply 
supported  or  clamped  at  the  boundaries,  and  subjected  to  a  uniform  axial 
compression  and  torsion  (individually  applied  or  in  combination)  is  analyzed. 

The  analysis  is  based  on  nonlinear  kinematic  relations,  linearly  elastic 
material  behavior,  and  the  usual  lamination  theory.  The  laminate  consists 
of  orthotropic  laminae,  which  typically  characterize  fiber  reinforced  com¬ 
posites.  Two  types  of  formulation  have  been  developed;  one  is  referred 
to  as  the  w,p- formulation,  based  on  Donnell- type  of  kinematic  relations. 

The  governing  equations  consist  of  the  transverse  equilibrium  equation  and 
the  in-plane  compatibility  equation.  These  two  equations  are  expressed  in 
terms  of  the  transverse  displacement,  w,  and  an  airy  stress  resultant  function, 

F.  The  other,  referred  to  as  the  u,  v,^- formulation,  is  based  on  Sanders 
type  of  kinematic  relations.  The  governing  equations  for  this  case  consist 
of  the  three  equilibrium  equations.  These  three  equations  are  expressed  in 
terms  of  two  ln-plane  displacement  components  u ,  v,  and  the  transverse  dis¬ 
placement  component, w .  Donnell's  type  of  shell  theory  approximation  can  be 
treated  as  a  special  case  in  the  u ,  v ,  w -formulation. 

Some  results  are  generated  for  certain  geometries  (isotropic  and  lami¬ 
nated)  and  these  serve  as  bench  marks  for  the  solution  scheme.  Results  are 
also  generated  for  composite  cylinders  by  changing  several  parameters.  The 
scope  of  these  parametric  studies  is  to  establish  the  effect  of  (a)  geometric 
imperfections,  (b)  lamina  stacking,  (c)  ln-plane  and  transverse  boundary  con¬ 
ditions  and  (d)  load  eccentricity  on  the  critical  conditions.  Moreover,  dynamic 
critical  loads  are  obtained  for  certain  configurations  under  axial  load  (sud¬ 
denly  applied). 
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CHAPTER  I 
INTRODUCTION 


Shell-like  structural  configurations  find  wide  uses  in  complicated 
aerospace  structural  systems.  Their  use  requires  sophisticated  analyses 
in  order  to  answer  questions  associated  with  their  behavioral  response 
to  external  loads  and  extreme  temperature  environments.  In  the  past  forty 
years  or  so,  numerous  investigations  addresses  themselves  to  several  specific 
questions  of  shell  behavior,  and  the  answers  to  these  questions  have  tremen¬ 
dously  enhanced  our  understanding  of  their  behavior.  All  of  this  was  done 
primarily  for  metallic  construction  of  these  configurations.  In  particular, 
attention  was  paid  to  the  degree  of  approximation  involved  in  the  use  of 
various  kinematic  relations  (which  led  to  several  linear  and  nonlinear  shell 
theories),  to  the  discrepancy  between  theory  and  experiment  for  the  buckling 
of  shells  (post-buckling  analyses  and  imperfection-sensitivity  studies),  to 
the  use  of  stiffening  for  shell  configurations  (including  eccentricity  effects) 
to  the  effect  of  support  conditions,  cutouts,  foreign  Inclusions  and  others. 
Moreover,  as  the  size  of  shell-like  structures  increased  and  as  the  computa¬ 
tional  capability  improved,  large  computer  codes  became  available,  for  the 
analysis  of  the  configurations. 

In  the  recent  few  years,  the  constant  demand  for  lightweight  efficient 
structures  led  the  structural  engineer  to  the  use  of  nonconventional  materials, 
such  as  fiber-reinforced  composites.  The  correct  and  effective  use  of  these 
materials  requires  good  understanding  of  the  system  response  characteristics 
to  external  causes  (loads,  properties  of  the  environment,  etc.).  Several 
research  programs  have  been  initiated  in  order  to  evaluate  the  physical 
properties  of  such  materials.  The  main  emphasis  in  these  studies  is  placed 
on  the  characterization  of  physical  properties  (finding  the  constants  in  the 


constitutive  relations  and  how  the  environment  affects  them).  In  addition, 
there  are  several  efforts  related  to  failure  criteria  and  failure-related 
effects,  such  as  scissoring  and  delamination. 

In  1975,  R.  C.  Tennyson  (1)  made  a  review  of  previous  studies  on  the 
buckling  of  laminated  cylinders.  According  to  Tennyson's  (1)  review,  perhaps 
one  of  the  earliest  stability  analyses  of  homogeneous  orthotropic  cylindrical 
shells  was  published  by  March  et  al.  (2)  in  1945.  After  that  time,  several 
theoretical  analyses  limited  to  orthotropic  shell  configurations  were  performed 
by  Schnell  and  Bruhl  (3),  Thielemann  et  al.  (4),  and  Hess  (5).  In  these  studies, 
simply  supported  end  conditions  were  partially  satisfied.  The  general  linear 
theoretical  solutions  to  anisotropic  cylinders  were  presented  by  Cheng  and  Ho 
(6)  (7),  Jones  and  Morgan  (8),  Jones  and  Hennemann  (9)  and  Hirano  (10).  Several 
papers  were  involved  in  the  comparison  of  the  efficiency  and  accuracy  between 
Flugge's  linear  shell  theory,  which  was  employed  by  Cheng  and  Ho  (6)  (7),  and 
other  shell  theories  (such  as  the  work  done  by  Tasi  (11),  Martin  and  Drew  (12) 
whose  theory  was  based  on  Donnell's  equations,  and  the  work  done  by  Chao  (13), 
whose  analysis  was  based  on  Timoshenko's  buckling  equations).  Stiffened  com¬ 
posite  cylinderical  shells  have  been  analyzed  by  Jones  (14).  Terebushdo  (15) 
and  Cheng  and  Card  (16).  Theoretical  analyses  of  the  effect  of  initial  geo¬ 
metric  imperfection  based  on  anisotropic  shell  theory  have  been  published  for 
the  loading  cases  of  pure  torsion  (17)  axial  compression  (18)  and  combined 
loads  (19)  (20).  Moreover,  several  computer  codes  (21-32)  (based  on  finite 
elements  and/or  differences)  that  deal  with  the  analysis  of  stiffened  shell 
configurations  have  been  modified  in  order  to  account  for  laminated  shell 
construction.  These  codes  do  serve  their  purpose,  and  that  is  that  they  are 
very  good  analytical  tools.  On  the  other  hand,  it  is  very  difficult,  if  not 
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possible,  to  use  these  codes  for  parametric  studies  or  for  evaluating 
the  applicability  and  limitations  of  various  shell  theories.  In  this 
report,  the  following  are  presented: 

(1)  The  mathematical  formulation  and  derivation  of  the  governing 
equations,  based  on  Donnell-type  (33)  nonlinear  kinematic  relations  and 
in  terms  of  the  transverse  displacement  component  and  an  Airy  stress 
(resultant)  function,  defined  in  the  text. 

(2)  The  mathematical  formulation  and  derivation  of  the  governing 
equations,  based  on  Sanders '-type  (34)  nonlinear  kinematic  relations  and 
in  terms  of  the  three  displacement  components  (small  strains  but  moderate 
rotations  about  in-plane  axes). 

(3)  Solution  schemes  for  both  formulations.  The  solution  methodology 
for  the  first  formulation  includes  post-limit  point  behavior,  while  the  so¬ 
lution  methodology  for  the  second  formulation  refers  only  to  the  pre- limit 
point  behavior  and  it  is  employed  to  estimate  critical  static  conditions 
(limit  point  loads).  The  listing  of  the  related  computer  codes  are  presented 
in  the  Appendices  of  this  report. 

(4)  Some  numerical  results  are  generated  (and  presented  herein)  with 
two  objectives  in  mind,  (a)  Some  serve  as  bench  marks  for  the  solution 
schemes  and  (b)  some  limited  parametric  studies  are  performed  in  order  to 
assess  effects  of  boundary  conditions  and  of  the  lamina  stacking  sequence, 
for  axially- loaded  laminated  cylindrical  shells. 

In  closing,  this  report  should  be  viewed  as  the  first  in  a  series  of 
reports  dealing  with  the  behavior  of  geometrically  imperfect,  stiffened  and 
laminated,  thin,  circular,  cylindrical  shells,  supported  in  various  ways 
(all  possible  extreme  cases  of  transverse  and  in-plane  boundary  conditions) 
and  subjected  to  static,  as  well  as  suddenly  applied,  destabilizing  loads. 
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CHAPTER  II. 


MATHEMATICAL  FORMULATION  AND  SOLUTION 
METHODOLOGY 

The  governing  equations  are  derived,  with  all  necessary  steps  shown  in 
detail,  in  Appendix  A.  The  geometry  is  a  thin,  circular,  geometrically  im¬ 
perfect  cylindrical  shell.  The  construction  consists  of  an  orthogonally  and 
eccentrically  stiffened  laminate  (each  lemina  is  orthotropic).  Note  that  a 
laminated  geometry,  an  eccentrically  stiffened  metallic  configuration  and  a 
metallic  shell  are  all  special  cases  of  the  construction  used  herein.  The 
stiffeners  are  uniform  in  geometry  and  with  constant  close  spacing,  which 
allows  one  to  employ  the  "smeared"  technique.  The  boundary  conditions  can 
be  of  any  transverse  and  in-plane  variety.  This  includes  free,  simply-sup¬ 
ported  and  clamped  with  all  possible  in-plane  combinations. 

The  loading  consists  of  transverse  (uniform  lateral  pressure)  and  eccentric 
in-plane  loads,  such  as  uniform  axial  compression  and  shear.  Eccentric 
means  that  the  line  of  action  of  these  loads  (applied  stress  resultants)  is 
not  necessarily  in  the  plane  of  the  reference  surface. 

In  the  derivation  of  the  governing  equations,  the  usual  lamination  theory 
is  employed.  Moreover,  thin  shell  theory  (Kirchhoff  -  Love  hypotheses  with 
two  different  approximation)  and  linearly  elastic  material  behavior  one  assumed. 
The  primary  assumptions  are  listed  in  Appendix  A.  On  the  basis  of  these 
general  assumptions  two  sets  of  field  equations  are  derived.  One,  referred 
to  as  the  w,F  formulation,  is  based  upon  Donnell-type  of  kinematic 
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relations.  For  this  case,  the  governing  equations  consist  of  the  transverse 
equilibrium  equation  and  the  in-plane  compatibility  equation.  These  two 
equations  and  the  proper  boundary  conditions  are  expressed  in  terms  of  the 
transverse  displacement  component,  w,  and  an  Airy  stress  resultant  function, 

F  .  The  second,  referred  to  as  the  u,  v,  w  -  formulations,  is  based  on 
Sanders'  type  of  kinematic  relations,  those  corresponding  to  small  rotations 
about  the  normal  and  moderate  rotations  about  in-plane  axes.  The  governing 
equations,  for  this  case,  consist  of  the  three  equilibrium  equations,  expressed 
in  terms  of  the  displacement  components  u,  v,  and  w.  Also,  the  proper  boundary 
conditions  are  expressed  in  terms  of  u,  v,  and  w.  In  this  formulation,  the 
Donnell  approximation  is  a  special  case  of  the  more  general  Sanders'  kine¬ 
matic  relations. 

The  solution  methodology  is  an  improvement  and  modification  of  the  one 
employed  and  described  in  Refs.  36  and  37.  For  details  the  reader  is  referred 
to  Appendix  A.  A  brief  description  of  the  solution  scheme  is  given  below  and 
only  for  the  wj  -  formulation. 

1) .  First,  a  separated  form  (fourier  series  type)  is  assumed  for  the 
dependent  variables  ,  w(x,y)  and  F(x,y).  In  addition  the  initial  geometric 
imperfection  is  also  expressed  in  a  similar  form. 

2) .  Next,  these  expressions  are  substituted  into  the  compatibility 
equations.  Use  of  trigonometric  identities  and  use  of  the  orthogonality  of 
the  trigonometric  functions  reduces  this  nonlinear  partial  differential 
equation  (compatibility)  into  a  system  of  (4k  +  1)  nonlinear  ordinary  dif¬ 
ferential  equations.  Furthermore,  use  of  the  Galerkin  procedure  in  connection 
with  the  equilibrium  equation  (in  the  circumferential  direction)  yields 

(2k  +  1)  additional  nonlinear  ordinary  differential  equations  in  the  (6k  +  2) 
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dependent  (on  x)  functions  needed  to  describe  the  response  of  the  system. 

Thus,  through  these  steps  the  two  nonlinear  partial  differential  equations 
are  reduced  to  a  set  of  nonlinear  ordinary  differential  equations. 

3) .  The  nonlinear  ordinary  differential  equations  are  reduced  to  a 
sequence  of  linear  systems  by  employing  the  generalized  Newton's  method 
(Ref.  38).  Iteration  equation  are  derived,  through  this,  based  on  the 
premise  that  a  solution  to  the  nonlinear  set  can  be  achieved  by  small  cor¬ 
rections  to  an  approximate  solution. 

4) .  Finally,  the  field  equations  (linearized  iteration  equations)  and 
the  corresponding  boundary  terms  (linear  set  of  equation)  are  cast  into  finite 
difference  form  by  employing  the  usual  central  difference  formula. 

Finally,  a  computer  program  has  been  written  (see  Appendix  B  for  Flow 
Charts  and  Program  Listings)  for  generation  of  results.  The  solution 
algorithum  is  a  modification  of  the  one  described  in  Ref.  43.  This  modification 
is  fully  described  in  Appendix  C. 
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CHAPTER  III 


DESCRIPTION  OF  STRUCTURAL  GEOMETRY 

Three  basic  configurations  are  used  in  generating  results.  The  consist 
of  a  four-ply  laminated  cylinder,  an  isotropic  cylinder  and  an  orthotropic 
cylinder.  All  configuration  are  geometrically  imperfect  but  the  imperfection 
in  either  symmetric  or  (virtually)  axisymmetric. 

The  laminated  geometries  considered  in  the  present  study  are  variations 
of  the  one  employed  in  (44).  This  reference  reports  experimental  results  for 
a  symmetric  angle-ply  laminate,  subjected  to  uniform  axial  compression  and 
torsion.  In  addition  some  isotropic  and  orthotropic  configuration  are  also  used. 

III.l  Laminated  Geometry 

For  the  laminated  geometries,  five  different  stacking  combinations  of 
the  4-ply  laminate  are  used  in  the  study. 

First,  the  common  geometric  and  structural  features  are:  each  lamina  is 
orthotropic  (Boron/Epoxy;  AVCO  5505)with  properties 

E„  =  2.o6?0X  /o'kN/m'i&x/O*  pSi),  Vn  -  0.21  ; 

Eaa  =  0.IU2X  Jo'kH/riW  l.)  i  R*t90.£C*.(l£i*.); 

G,2  -  0.0W2X/0*  kN/m'  ( 0,6SX/0*p  Si);  L-S?l  C#t.  {/Sin.); 


(  /  ) 


hpiy  *  o.o /U62  cnio.  oo$*7n.) 

(.hpiy  -  hk~  h)t-i  ;  k  *  / .  2 ,  3 ,4 ;  fda*  ) 


The  five  different  stacking  combinations  are  denoted  by  I  -  i,  i 


1,  2.. 5,  and  correspond  to 

I-  I  :  4S’/-4S’/-4S’M$‘  i  1-2  •  4S’/-4S°/4S/-4S  ;  13  -  -  12 
I  -  4 :  9o‘/io/io’/0  ‘  ;  I-ff:  0’/30’/6o’/?o‘  U) 

Where  the  first  number  denotes  the  orientation  of  the  fibers  of  the  out¬ 
most  ply  with  respect  to  x,  and  the  last  of  the  innermost.  Geometry  1-1  is 
a  symmetric  one  and  it  corresponds  to  that  of  (44).  Geometries  1-2  and  1-3 
denote  antisymmetric  regular  angle-ply  laminates,  while  geometries  1-4  and 
1-5  are  completely  asymmetric. 

111. 2  isotropic  Geometry 

The  isotropic  cylinder  has  the  following  geometric  and  structural  fea¬ 
tures  (aluminum  alloy) 

E  =  7MX/07kN/ml  ( /o.s x/o‘ psi ) ;  v=o.$ 

R  -  /0./6CW.  (4in)  j  L/x  =  /  ;R/h  ~/ooo  (3; 

111. 3  Orthotropic  Geometry 

Finally,  the  properties  of  the  orthotropic  configuration  are  (single 
0°  -  ply  shell  made  of  the  Boron/Epoxy  material) 

Bxx  =  2 Ml / /0*  kN/ri  (30  X/06  psi)  ;  i)My  5  0.2 1 
Eyy  =  0.  IU2  XlO^kN/on  C  -2.7  X/o6  psi ) 

Gxy  -  0.04422  X/0*  KN/mL  0.  i$X/o*PSi) ,  R  =  !  90S  CM.  (7 Sin.) 

L  -  381'OOn.  (/Sin.)  ;  i-  0.0S3SSCM.  (0.02i2in.)  (4. 


III. 4  Imperfection  Shapes 

Two  Imperfection  shapes  are  used  in  the  study,  one  which  is  symmetric, 
and  one  which  is  virtually  axisymmetric 


SytuuJu'c  :  W °Oc.y)  S  SA  font 


(5) 


a*uys»u&<<  ••  w  (x.y;  =  Ski-  +0.1  s0th-~Cfi4*ilr 


(6) 


where  5  is  a  measure  of  the  Imperfection  amplitude.  Note  that 

symmetric  imperfection,  Eq.  (  5),  5  *  w°  _  /h,  while  for  the 

max 

axisymmetric  imperfection,  Eq.  (6  ),  l  *  w°  _  /l.lh. 

max 
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(virtually) 


e 


e 


e 


9 


CHAPTER  IV 


NUMERICAL  RESULTS  AND  DISCUSSION 

Numerical  results  are  generated,  for  the  geometries  described  in  the 
preceeding  chapter,  using  the  W-F  formulation,  for  two  load  cases:  (a)  uni¬ 
form  axial  compression  and  (b)  torsion.  The  loads  are  applied  individually 
and  in  combination.  The  results  consist  of  finding  pre-  and  post-limits  point 
behavior,  as  well  as  critical,  conditions  for  static  and  dynamic  (sudden- 
some  results)  application  of  the  loads. 

The  generated  results  serve  a  multitude  of  purpose.  Some  results  serve 
as  bench  marks  for  the  solution  methodology  and  the  computer  code.  These 
results  are  compared  with  already  known  and  accepted  numbers.  Some  results 
correspond  to  parametric  studies,  tftlch  are  performed  in  order  to  enhance  our 
understanding  of  the  behavior  of  laminated  shells.  The  effects  of  lamina 
stacking  on  critical  conditions  is  studied.  Furthermore,  the  effect  of  in¬ 
plane  and  transverse  boundary  conditions  on  critical  loads  is  evaluated  for 
some  geometries.  Moreover,  the  imperfection  sensitivity  is  fully  assessed 
for  all  geometries.  Dynamic  critical  loads  are  obtained  for  very  few  geome¬ 
tries.  Most  of  the  generated  results  are  presented  in  tabular  and  graphical 
form.  All  generated  results  are  not  presented,  herein,  for  the  sake  of  bre¬ 
vity.  The  conclusions,  though,  are  based  on  all  generated  data. 

IV.  1.0  Axial  Compression 

Several  studies  are  performed  for  this  load  case.  Each  one  of  these 
studies  is  described  and  discussed  separately. 


For  this  study,  the  load  is  applied  through  the  reference  surface  (which 
is  the  inidsurface  of  the  laminate)  and  the  boundary  conditions  are  SS-3 
(classical  simply  supported).  The  imperfection  shape  is  symmetric,  Eq.  (5  ) 

Table  4-1  shows  critical  loads,  N^  (limit  point  loads),  for  each  geo- 

xx 

metry  and  various  values  of  the  imperfection  amplitude  parameter  ,  5.  It 
also  presents  the  range  of  n-values  used  in  finding  critical  loads,  and  the 
n-value  corresponding  to  the  critical  condition.  These  results  are  also  pre¬ 
sented  graphically  on  Fig.  4.1. 

Geometry  1-1  is  the  one  reported  in  (44).  According  to  this  reference, 

the  classical  (linear  theory)  critical  load  is  165  lbs. /in  (N  )  and  the 

“cl. 

experimental  value  is  106  lbs. /in.  Note  from  Fig.  4.1  that  through  extra¬ 
polation  N^  at  5  =  0  is  approximately  equal  to  148  lbs. /in.,  which  is  10% 
lower  than  the  reported  [44]  classical  value. 

The  results  for  geometries  1-2  and  1-3  are  identical.  Both  geometries 
are  antisymmetric.  This  is  reasonable  since  (a)  the  imperfection  shape 
is  symmetric  with  respect  to  a  diametral  plane  and  (b)  the  axially- loaded 
cylinder  does  not  distinguish  between  a  positive  45°  direction  and  a  nega¬ 
tive  45°  direction. 

Moreover,  for  virtually  the  entire  range  of  g -values  considered,  the 
1-2(3)  geometry  seems  to  be  the  weakest  configuration,  while  the  asymmetric 
configuration  corresponding  to  1-5  is  the  strongest.  The  order  of  going 
from  the  weakest  to  the  strongest  is  1-2(3),  1-1,  1-4  and  1-5.  Note  that 
1-5  is  a  geometry  for  which  the  0°  -ply  is  on  the  outside.  Now  since  buck¬ 
ling  occurs  in  an  inward  transverse  displacement  mode  (w  is  positive),  then 
the  outside  layer  is  in  compression  and  it  is  reasonable  to  expect  the  strong' 
est  configuration  to  correspond  to  1-5,  the  fibers  of  the  outer  ply  are  in 
the  longitudinal  direction. 


Furthermore,  the  difference  between  1-4  and  1-5  geometries  is  the  order 
of  stacking  (one  is  the  reverse  of  the  other).  Their  behavior,  then,  can  be 
compared  to  the  behavior  of  orthogonally  stiffened  metallic  shells  with  out¬ 
side  and  inside  stiffening.  Geometry  1-5  is  comparable  to  outside  stiffening, 
while  geometry  1-4  to  inside. 

Figs.  4.2  and  4.3  present  typical  equilibrium  paths  for  all  geometries. 
Fig.  4.2  corresponds  to  geometry  1-1,  while  Fig.  4.3  to  geometry  1-4.  As 
seen,  the  response  is  in  terms  of  plots  of  applied  load  versus  average 
end  shortening,  e^.  It  includes,  pre-limit  point  behavior,  limit  points 
and  post-limit  point  behavior,  for  each  5 "value.  The  entire  curves  corre¬ 
spond  to  the  same  wave  number,  n.  This  n-value  is  the  one  that  yields  cri¬ 
tical  conditions  (the  one  at  the  instant  of  buckling).  If  a  clear  picture 
of  post-limit  point  behavior  is  desired,  one  should  show  the  plots  that  cor¬ 
respond  to  other  wave  numbers .  This  would  possibly  reveal  that  the  post¬ 
limit  point  curves  cross  each  other,  as  in  the  case  of  isotropic  shells  (46). 

Finally,  for  the  two  asymmetric  configurations,  1-4  and  1-5,  critical 
dynamic  loads  are  calculated  of  the  entire  g -range  (see  Fig.  4.4).  These 
are  obtained  by  employing  the  criteria  described  in  (46,  39),  and  they 
correspond  to  lower  bounds  of  critical  conditions  when  the  axial  compression 
is  applied  suddenly  with  Infinite  duration.  According  to  this  criterion 
and  methodology  for  estimating  critical  dynamic  conditions,  when  5  *  0 
(perfect  geometry)  the  static  and  dynamic  critical  loads  are  the  same.  As 
the  Imperfection  amplitude  increases  the  dynamic  loads  are  smaller  than  the 
static  loads.  For  these  geometries,  1-4  and  1-5,  and  0<g<2.0,  the  dynamic 
critical  load,  is  never  smaller  than  607.  of  the  corresponding  static 
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IV.  1.2  Effect  of  Boundary  Conditions 

The  effect  of  both  transverse  and  in-plane  boundary  conditions  are 
assessed. 

Results  are  also  generated  for  the  isotropic  geometry  (aluminum  alloy) 
and  various  in-plane  boundary  conditions.  These  serve  as  bench  marks  for 
the  solution  scheme,  and  the  results  are  presented,  in  part,  on  Table  4.2 
and  Fig.  4.5.  For  this  geometry  the  shape  of  the  imperfection  is  taken  to 
be  axisymmetric,  Eq.  (  6  ).  On  Table  4.2,  the  n-value  that  corresponds  to 
the  critical  load  is  given  in  brackets.  Note  that  for  small  g -values  (see 
Fig.  4.5),  the  trend  is  exactly  that  suggested  by  Hoff  and  Ohira,  indepen¬ 
dently  (see  (47)),  i.e.,  the  weakest  configuration  is  SS-1,  the  next  one 
SS-2,  while  SS-3  and  SS-4  yield  the  classical  results.  Note  also  that, 
through  extrapolation,  (as  g  —  0) ,  the  present  results  agree  with  those  of 
(47).  For  SS-1  the  ratio  of  critical  load  to  classical  load  is  0.55,  for 
SS-2  0.68,  and  for  SS-3  and  SS-4  0.98.  Clearly  here  (isotropic  case)  the 
geometry  for  boundary  conditions  SS-1  and  SS-2,  is  not  very  sensitive  to 
geometric  imperfection,  while  for  SS-3  (primarily)  and  SS-4,  it  is.  Note 
that,  for  small  g-values,  the  v  =  const,  in-plane  boundary  conditions  (SS-3 
and  SS-4)  yield  a  stronger  configuration.  For  higher  5 -values  the  stronger 
configuration  corresponds  to  u  =  const,  in-plane  boundary  conditions  (SS-2 
and  SS-4). 
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Tattle  4.2  Effect  of  In-Plane  Boundary  Condition  on  Critical  Load 
(Isotropic  Geometry,  Simply  Supported  Case). 


SS-1  1  SS-2 


Note  that,  no  attempt  is  made  here  to  find  the  shape  of  the  imperfection 

that  yields  the  lowest  critical  load.  For  the  case  of  the  laminated  shell, 

the  imperfection  amplitude  parameter,  g,  is  varied  from  0.05  to  two.  The 

first  number,  0.05,  corresponds  to  a  virtually  perfect  geometry  shell,  while 

the  second  number  (two)  denotes  an  amplitude  in  the  neighborhood  of  two  shell 

thicknesses  (this  is  considered  very  large  for  thin  construction). 

In  order  to  establish  the  imperfection  sensitivity  of  the  laminated  shell 

and  the  effect  of  boundary  conditions  on  the  limit  point  load  (critical  load) , 

geometry  1-5  is  employed,  along  with  a  symmetric  type  of  imperfection,  Eq.  (  5  ). 

As  already  established,  geometry  1-5  yields  the  strongest  configuration  for 

SS-3,  by  comparison  to  all  other  geometries  (I-i,  i  =  1,  2,  3,  4). 

Table  4.3  Effect  of  Boundary  Conditions  on  Critical  Loads.  (Laminated  Geometry  1-5) 


SS-1 

n«7 


0.50  26. 


SS-2 

n«8 


32.39 


N  ,  kN/m  (lbs /in) 
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SS-4 

n-9 
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Table  4.3  lists  critical  loads  for  various  boundary  conditions  and  5 -values 
(f  =  w°  /h;  for  this  case) .  The  value  of  n  denotes  the  number  of  full  waves 
around  the  circumference  at  the  instant  of  buckling.  These  results  are  shown 
graphically  on  Figs.  4.6  and  4.7.  A  number  of  observations  are  made.  First, 
for  low  g-values  (see  Fig  4.6)  SS-3  and  SS-4  yield  stronger  configurations 
than  SS-1  and  SS-2.  For  higher  values  of  5,  SS-2  and  SS-4  yield  stronger 
configurations  than  SS-1  and  SS-3.  Another  way  of  stating  the  same  thing  is 
that  for  low  5 -values  the  v  =  const,  in-plane  boundary  condition  yields  a 
stronger  configuration,  while  for  higher  g-values  the  u  =  const,  in-plane 
boundary  condition  yields  higher  critical  loads.  This  conclusion  is  the 
same  for  isotropic  geometries.  On  the  other  hand,  for  the  clamped  case, 

CC-2  and  CC-4  (u  =  const.)  yield  stronger  configurations  than  CC-1  and  CC-3 
for  the  entire  5 -range  considered.  Another  observation  is  that  for  SS-1 
and  SS-2  the  geometry  is  not  as  sensitive  to  initial  geometric  imperfections 
as  it  is  for  SS-3,  SS-4,  and  CC-i  (i  =  1,  2,  3,  4)  [see  Figs.  4.6  and  4.7]. 

It  is  also  worth  mentioning  that  a  comparison  between  the  values  at  5  =  0 
between  SS-1  and  SS-4  is  reminiscent  of  what  happens  in  the  isotropic  case 
(the  critical  load  for  SS-1  is  virtually  half  the  value  of  that  for  SS-4). 


IV.  1.3  Effect  of  In-plane  Load  Eccentricit 


Next,  the  effect  of  load  eccentricity  is  assessed.  In  all  configurations 
for  which  results  are  generated,  the  shell  midsurface  is  taken  as  the  reference 
surface.  Then  it  is  assumed  that  the  uniform  axial  compression  is  applied 
eccentrically,  which  induces  a  bending  moment  at  the  boundary,  M  *  E  N^ 

[see  Eqs  A-35  &  A-3Q.  Note  that  this  load  eccentricity  affects  only  the 
simply  supported  boundary  conditions. 
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Results  are  generated  and  presented  for  the  isotropic  geometry,  ortho - 
tropic  geometry,  and  laminated  1-1,  1-4  and  1-5  geometries,  using  a  symmetric 
imperfection  shape  Eq.  (  5  ),  and  classical  simply  supported  boundary  condi¬ 
tions  SS-3. 

These  results  are,  in  part,  presented  on  Tables  4. 4-4. 6. 

One  might  expect  a  negative  edge  moment  (corresponding  to  positive  load 
eccentricity)  to  have  a  stabilizing  effect  on  an  axially-load  cylindrical 
shell,  regardless  of  the  construction.  Contrary  to  this,  the  generated  re¬ 
sults  do  not  support  the  expectation.  For  small  eccentricities  (-0. 5<E/h<0.5) 
and  isotropic  geometry  (see  Table  4.4)  the  response  seems  to  be  insensitive 
to  the  eccentric  application  of  the  load.  This  is  true  for  both  imperfection 
shapes  [axisymmetric  and  symmetric,  Eq.  (5  )  &  (  6  )3. 

Table  4*4  Effect  of  Load  Eccentricity  (Isotropic  &  Orthotropic) 


Imperf . 
Shape 
& 

Geometry 

\  e/h 

V- 
u  \ 

N  ^  in  kN/m  (lbs/in.) 
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Table  4.5  Effect  of  Load  Eccentricity  (Laminated  I-l  Geometry) 


s' 
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in  kN/m  (lbs/in.) 
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E/h  ■ 

0.5 

E/h  • 

•  0 
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nm 

Axiaym. 
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cm 
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ns 
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SB 
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■m 
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SS-4  boundary  conditlona  and  n  ■  6 


Table  Effect  of  Load  Eccentricity  (Laminated  1-4  and  1-5  Geometries; 
Symmetric  Imperfection;  SS-3  boundary  conditions).  ? 


% 

in  kN/m  (lbs/in.);  n  ■  8 

• 

| 
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For  very  large  eccentricities  (|E/h|  >  12),  positive  eccentricity  has  a 
stabilizing  effect,  tfiile  negative  eccentricity  has  a  destabilizing  effect. 

In  the  intermediate  range  an  irregularity  is  observed.  It  was  suspected  that 
one  possible  reason  for  this  behavior  may  be  attributed  to  the  Poisson  effect. 

As  the  load  is  applied,  quasistatically,  the  midportion  of  the  shell  moves  out¬ 
ward  because  of  the  Poisson  effect;  it  reaches  a  maximum  expansion,  before  the 
load  reaches  its  critical  value,  and  then  an  inward  motion  takes  place,  and 
finally  at  and  after  collapse  this  inward  motion  continues.  This  sequence  of 
events  and  the  corresponding  stabilization  or  destabilization  of  the  load  ec¬ 
centricity  is  heavily  dependent  on  the  value  of  Poisson's  ratio  or  the 
term  in  the  extensional  stiffness  matrix.  For  instance,  some  data  are  gene¬ 
rated,  for  the  isotropic  geometry  (£  "  0.5;  SS-3  and  axisyymnetric  imperfection) 

l 

but  with  v  ■  0.1.  The  limit  point  loads,  N  (critical  load)  for  three  values 
of  eccentricity  (E/h)  are:  3.305  kN/m  (18.88  lbs/in)  for  E/h  *  +  0.5;  2.76 
kN/m  (15.81  lbs/in.)  for  E/h  -  0;  and  2.745  kN/m  (15.68  lbs/in)  for  E/h  =  -0.5. 
This  clearly  shows  that  positive  eccentricity  has  a  stabilizing  effect.  This 
observation  is  also  true  for  the  orthotropic  geometry  (see  Table  4.4)  for  which 
the  value  of  A^  is  small  by  comparison  to  Aj^.  On  the  other  hand,  for  v  =  0.3 
and  the  laminated  geometries  for  which  the  values  of  A^  are  the  same  order 
of  magnitude  as  An*  it  cannot  be  said  that  positive  eccentricity  has  a  stabi¬ 
lizing  effect  (see  Tables  4.5  and  4.6).  In  reality,  for  these  geometries  no  de¬ 
finite  conclusion  should  be  drawn  regarding  stabilization  through  load  eccen¬ 
tricity  (or  applied  edge  moment).  It  is  worth  observing,  though,  that  for  all 
laminated  geometries  (see  Tables  4.5  and  4.6),  whatever  the  effect  is,  it  does 
diminish  with  increasing  amplitude  of  imperfection. 
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IV.  2.0  Torsion  with  and  without  Axial  Compression 

For  this  particular  load  case,  in  addition  to  the  axisymmetric  shape  for 
the  geometric  imperfection,  two  additional  shapes  are  employed  in  the  studies. 

These  additional  shapes  correspond  to  approximations  of  the  linear  theory 

.  n 

(see  Appendix  D)  buckling  modes  for  positive  and  negative  torsion  for  all 
five  geometries. 

In  particular.  Appendix  C  deals  with  solutions  to  the  linearized  buck¬ 
ling  equations  for  the  case  of  pure  torsion.  To  this  end,  the  Galerkin  pro¬ 
cedure  is  employed  and  the  following  approximation  is  employed  for  the  buck¬ 
ling  modes 

(AinC*2#  +Bin4*n2n 


L 

atm 


(7) 


Because  of  orthogonality,  only  one  n-value  is  needed  In  Appendix  D  a  ten- 
term  approximation  (M*5)  is  obtained  for  all  five  geometries.  By  stu-yin? 
the  results,  one  two-term  approximation  for  positive  torsion,  w°(+) ,  and  one 
two-term  approximation  for  negative  torsion,  w°(-),  for  all  five  geometries 
are  used  in  this  study.  The  various  coefficients  are  first  normalized  with 
respect  to  Eq.  (7  ),  and  then  adjusted  such  that  the  maximum  amplitude 
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The  generated  results  for  this  case  are  presented,  in  part,  both  in  tabu¬ 
lar  and  graphical  forms.  The  discussion,  though,  and  the  related  conclusions 
are  based  on  all  data. 

First,  Table  4.7  shows  values  of  critical  torsion,  N  ,  for  the  two  asym- 

xy 

metric  imperfection  shapes,  Eqs.  (  8  )  and  (  9  )  (corresponding  perfect  geo¬ 
metry  buckling  modes  for  positive  and  negative  torsion)  and  several  values  of 
the  imperfection  amplitude  paraowter.  The  torsion  is  applied  in  both  directions 
and  the  critical  values  are  recorded.  The  corresponding  minimizing  value  of  n 
(number  of  full  waves)  is  shewn  in  parenthesis. 

Note  that  the  linear  theory,  perfect  geometry  critical  values  (from  Ap¬ 
pendix  d)  for  geometry  1-1  are  39.9  lbs. /in.  for  positive  torsion,  and  -75.5 
lbs./ln.  for  negative  torsion.  Moreover,  the  experimental  results  obtained 
from  (44)  for  this  geometry  (1-1)  are  26.5  lbs. /in.  for  negative  torsion. 

Note  that  the  construction  (orientation  of  the  plies)  is  such  that  the 
configuration  is  much  weaker  when  loaded  in  the  negative  direction,  regardless 
of  which  of  the  two  imperfection  shapes  is  used.  Furthermore,  when  w°(+)  is 
present  the  configuration  is  somewhat  sensitive  for  positive  torsion  (see  second 

column  at  5  ■  0.10,  N  ■  35.32  sensitive  for  negative  torsion  (see  third 

xy 

column).  On  the  other  hand,  when  w°(-)  the  reverse  is  true,  l.e.  the 
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configuration  is  insensitive  for  positive  torsion  (fourth  column)  and  rather 
sensitive  for  negative  torsion  (last  column).  Note  that  the  experimental  values 
(+  26.5  lbs. /in.  and  -65.72  lbs. /in.),  compare  well  with  the  theoretical  values. 
Note  that  the  tested  specimen  (44)  is  of  unknown  imperfection  shape  and  ampli¬ 
tude  . 

Next,  Table  4.8  presents  critical  shear  stress  resultants  (and  the  mini¬ 
mizing  n-value  in  parenthesis)  for  all  five  geometries  and  an  imperfection 
shape  similar  to  the  positive  torsion  buckling  mode  of  the  perfect  geometry, 

Eq.  (8  ).  These  results  are  shown  graphically  on  Fig.  4.8).  Note  that  the 

strongest  configuration  corresponds  to  1-5,  while  the  weakest  to  the  symmetric 
geometry  1-1 .  This  conclusion  holds  true  for  the  imperfection  shape  used,  w°(+). 

It  is  worth  observing  that  the  regular  angle-ply  antisymmetric  geometries, 
1-2  and  -1-3,  yield  virtually  the  same  strength  for  positive  torsion  and  w°(+). 
Moreover,  geometry  1-4  is  much  weaker  by  comparison  to  the  other  asymmetric 
geometry  (1-5)  but  not  as  weak  as  the  symmetric  geometry.  These  observations 
are  reminiscent  of  the  old  external  versus  internal  positioning  of  the  ortho¬ 
gonal  stiffeners  controversy  concerning  metallic  stiffened  configurations.  In 
relation  to  this,  in  the  case  of  orthogonally  stiffened  complete  spherical 
shells  subjected  to  uniform  pressure  (see  Ref.  48)  it  is  observed  that  the 
weakest  conf iguration  corresponds  to  zero  (or  close  to  it)  stiffener  eccen¬ 
tricity,  and  the  strength  of  the  stiffened  sphere  increases  as  the  eccentricity 
increases  in  either  direction  (inward  or  outward).  Thus,  one  can  conclude  from 
Fig.  4.8  that  all  five  configurations  are  imperfection  sensitive,  but  not  as 
sensitive  as  they  are  for  the  case  of  uniform  axial  compression  (See  Fig.  4.1). 
This  conclusion  is  in  line  with  the  behavior  of  metallic  cylindrical  shells 
with  or  without  stiffening  members. 
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In  Ref.  44,  experiments  are  conducted  for  geometry  1-1,  to  determine  the 
interaction  curve  that  separates  the  stable  from  the  unstable  region  between 
uniform  axial  compression  and  torsion.  Because  of  this,  numerical  results  are 
obtained  fro  geometry  1-1  and  two  imperfection  shapes.  One  is  virtually  axi- 
symmetric,  Eq.  (6  ),  and  one  similar  to  the  (positive  torsion)  perfect  geo¬ 
metry  buckling  mode,  Eq.  (8  ).  The  theoretical  interaction  curves  are  gene¬ 
rated  for  several  values  of  the  imperfection  amplitude  parameter,  by  the 
following  steps.  First,  the  critical  value  for  pure  torsion  is  obtained. 

Then,  starting  with  zero  torsion  and  several  values  of  the  applied  shear  stress 
resultant,  but  smaller  than  the  critical  pure  torsion  the  corresponding  cri¬ 
tical  axial  compression  is  obtained.  In  each  combination  a  study  of  the  ef¬ 
fect  of  n  is  performed.  The  results  are  presented  in  tabular  form  on  Tables 
4.9  and  4.10  and  graphically  on  Figs.  4.9  and  4.10. 

The  data  of  Table  4.9  are  plotted  on  Fig.  4.9  and  of  Table  4.10  on  Fig.  4.10 

On  both  figures  the  experimental  (44)  Interaction  curve  is  shown  by  the  dash¬ 
ed  line.  Not  knowing  what  the  imperfection  shape  and  amplitude  of  the  tested 
cylinder  are,  these  plots  may  suggest  a  reasonable  comparison  between  theory 
and  tes  t . 
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IV.  3.0  CONCLUSIONS 


All  of  the  conclusions  are  based  on  the  generated  results,  which  are 
obtained  by  the  W,  F-formulation.  No  results  have,  as  yet,  been  generated 
by  the  u,  v,  w- formulation. 

From  all  results,  one  may  list  the  following  as  the  most  noteworthy  ccn- 
c lus ions . 

1.  Buckling,  for  all  configurations,  is  of  the  violent  type  (snap  through 
buckling  through  limit  point  instability). 

2.  For  SS-3  boundary  conditions  and  axial  compression  with  zero  eccentricity, 
the  strongest  configuration  corresponds  to  the  asymmetric  congiguration, 

1-5,  while  the  weakest  configuration  corresponds  to  the  antisymmetric 
configurations,  1-2  and  1-3. 

3.  Again  for  SS-3  and  axial  compression,  the  dynamic  critical  loads  (lower 
bounds,  \rtien  the  corresponding  static  loads,  but  their  values  are  never 
smaller  than  60%  of  the  static  critical  loads. 

4.  The  average  end  shortening  (for  axial  compression),  corresponding  to  the 
limit  point  for  the  same  |-value,  is  smaller  for  the  asymmetric  geometries 
(1-4,  1-5)  than  for  the  symmetric  (1-1 )  and  antisymmetric  (1-2  and  1-3) 
geometries  by  almost  a  factor  of  three. 

5.  For  the  isotropic  geometry  (SS-i  boundary  conditions) 

5a:  For  the  perfect  configuration  and  very  small  imperfections,  the  effect 
of  in-plane  boundary  conditions  is  such  that  SS-3  and  SS-4  (v  =  const.) 
yield  stronger  configurations  than  SS-1  and  SS-2  (N  =  -F  =  0) 

>*y 

For  higher  values  of  the  imperfection  amplitude,  % ,  SS-2  and  SS-4 
(u  *  const.)  yield  stronger  configurations  than  SS-1  and  SS-3 


5b: 


6.  For  the  laminated  geometry,  the  effect  of  in-plane  boundary  conditions  for 
SS-i  is  the  same  as  for  the  isotropic  geometry.  For  clamped  boundaries, 

CC-2  and  CC-4  (u  =  const.)  yield  stronger  configurations  that  CC-1  and 
CC-3,  for  the  entire  5 -range. 

7.  For  both  geometries,  1-5  and  isotropic,  the  sensitivity  to  initial  geometric 
imperfection  is  dependent  upon  the  in-plane  boundary  conditions  for  SS-i. 
When  v  =  const  (SS-1  and  SS-2),  the  geometries  are  not  very  sensitive.  On 
the  other  hand,  when  u  =  const  the  geometries  are  very  sensitive. 

8.  As  far  as  the  effect  of  load  eccentricity  on  critical  loads  is  concerned, 
no  general  conclusion  can  be  drawn.  But  whatever  the  effect  is  (stabi¬ 
lizing  or  destabilizing  for  a  given  geometry) ,  it  diminishes  with  in¬ 
creasing  value  of  the  imperfection  amplitude  parameter  (£ -values). 

9.  When  loaded  in  pure  torsion,  the  strongest  configuration  corresponds  to 
geometry  1-5  (asymmetric) ,  while  the  weakest  corresponds  to  the  symmetric 
geometry  1-1,  for  the  imperfection  shape  corresponding  to  the  positive 
torsion  buckling  mode,  w°(+). 

10.  Geometry  1-1  is  weaker  when  loaded  in  the  positive  direction  than  when 
loaded  in  the  negative  direction  regardless  of  the  imperfection  shape 
(for  all  that  were  employed). 

11.  When  loaded  in  pure  torsion,  laminated  shell  configurations  are  sensitive 
to  initial  geometric  imperfections,  but  not  as  sensitive  as  when  loaded 
in  axial  compression. 

12.  Comparison  between  theoretical  predictions  (corresponding  to  various  im¬ 
perfection  amplitudes  and  shapes)  and  experimental  results  is  reasonably 
good. 
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APPENDIX  A 


MATHEMATICAL  FORMULATION 

A.  1.0  Introduction 

The  governing  equations  are  derived,  in  this  section,  for  the  following 
geometry  and  loading.  The  thin,  circular,  cylindrical  shell  is  assumed  to 
be  geometrically  imperfect.  The  construction  is  laminated  (each  lamina  is 
orthotropic)  and  in  addition,  the  shell  is  orthogonally  and  eccentrically 
stiffened.  The  stiffeners  are  uniform  and  with  uniform  close  spacing, 
which  allows  one  to  employ  the  "smeared"  technique.  The  boundary  conditions 
can  be  of  any  transverse  and  in-plane  variety.  This  includes  free,  simply- 
supported  and  clamped  with  all  possible  in-plane  combinations.  The  loading 
consists  of  transverse  (uniform  lateral  pressure)  and  eccentric  in-plane 
loads,  such  as  uniform  axial  compression  and  shear.  Eccentric  means  that 
the  line  of  action  of  these  loads  (applied  stress  resultants)  is  not 
necessarily  in  the  plane  of  the  reference  surface.  In  the  derivation  of  the 
governing  equations,  the  usual  lamination  theory  is  employed.  Moreover, 
thin  shell  theory  (Kirchhoff-Love  hypotheses)  and  linearly  elastic  behavior 
are  assumed.  The  primary  assumptions  are  listed  below: 

(1)  The  shell  is  thin  (total  smeared  thickness  is  much  smaller  than 
the  initial  average  radius  of  curvature-cylinder  radius). 

(2)  Normals  remain  normal  and  inextensional. 

(3)  The  strains  are  small,  the  rotations  about  the  normal  are  small 
and  the  rotations  about  in-plane  axes  are  moderate. 

(4)  The  Imperfection  shape  is  such  that  the  initial  curvature  is  small 
[R|w°.ii  «U  i  ■  x,y]. 

(5)  The  stiffness  are  along  principal  directions. 

(6)  The  s tiffener- laminate  connections  are  monolithic. 


(7)  The  stiffeners  do  not  carry  shear;  shear  is  entirely  trans¬ 
mitted  by  the  laminate  . 

(8)  The  stiffeness  are  torsionally  weak  and  thus  they  do  not  con¬ 
tribute  to  the  shell  twisting  stiffness  (the  equations  and  related  pro¬ 
grams  can  easily  be  changed  to  accomodate  the  case  of  torsionally  strong 
stiffeners) . 

On  the  basis  of  these  general  assumptions,  two  sets  of  field  equations 
are  derived.  One,  referred  to  as  the  w,  F  -  formulation,  is  based  on 
Donnell-type  of  kinematic  relations.  The  governing  equations  consist  of 
the  transverse  equilibrium  equation  and  the  in-plane  compatibility  equation. 
These  two  equations  and  the  proper  boundary  conditions  are  expressed  in  terms 
of  the  transverse  displacement  component,  w,  and  an  Airy  stress  resultant 
function,  F.  The  second,  referred  to  as  the  u,  v,  w  -  formulation  is  based 
on  Sanders'  type  of  kinematic  relations,  those  corresponding  to  small  rotations 
about  the  normal  and  moderate  rotations  about  in-plane  axes.  The  governing 
equations  for  this  case  consist  of  the  three  equilibrium  equations.  These 
equations  are  expressed  in  terms  of  the  three  displacement  components,  u,  v 
and  w.  Also,  the  proper  boundary  conditions  are  expressed  in  terms  of  u,  v, 
and  w.  The  corresponding  Donnell  approximation  appears  as  a  special  case  of 
the  more  general  Sanders'  kinematic  relations.  The  derivation  along  with  all 
necessary  relations  are  presented  separately  for  each  formulation. 

A.  2.0  The  w,  F  -  Formulation 

The  geometry  and  sign  convention  for  this  formulation  are  shown  on  Figs. 
A.l  and  A  .2. 

The  topics  of  kinematic  relations,  stress  and  moment  resultants,  governing 
equations,  boundary  conditions  and  solution  procedure  are  treated  separately. 
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A*  2.1  Kinematic  Relations 

Let  w°  be  measured  from  the  perfectly  cylindrical  surface  to  the 
refer  surface  of  the  laminated  shell.  Let  w  denote  the  transverse 
displacement  component  of  reference  surface  material  points  and  be 
measured  from  the  undeformed  surface.  Let  u  and  v  denote  the  usual 
in-plane  displacement  components  along  the  x  and  y  directions  respective¬ 
ly. 

The  Donnell-type  (33)  kinematic  relations  are  given  by 

=  ~  H  k*m 

~  ~  Z.hyy  (A-l  ) 

Y.y  =  Y’-2Zk,j 

where  the  superscript  "o"  denotes  reference  surface  strains  and  the  k's 
denote  the  reference  surface  changes  in  curvature  and  torsion.  Note  that 
the  positive  z-direction  is  inward  (see  Fig.  A.l). 

According  to  Donnell  the  e°'s  and  k's  are  related  to  the  displacement 
components  by 


€j^  =  U,K  4-£  */,x  4  K/,k  vs/.; 


(A -2) 


foy  =  W.y  4  lr,t  +  *>y  4  VS/.X  M/.y  4  */,y 

Xtt-  ~  (  )>x  = 

3  ^>y,y  =  (  W,y),f  =  Wiry  (/A  -3,) 

“  'V.xy 

A*  2.2  Stress-strain  Relations 

Each  lamina  is  assumed  to  be  orthotropic  and  the  directions  of 
orthotropy  (1*2)  make  an  angle  6  with  the  in-plane  axes  (x,y). 
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The  orthotropic  constitutive  (it  is  assumed  that  the  generalized 
Hooke's  law  holds)  relations  for  the  kth  lamina  are  given  below.  Note 


that  for  an  n-ply  laminate  k  varies  from  one  to  n,  and  the  first  ply 
(or  lamina)  is  on  the  outside,  while  the  nth  ply  is  on  the  inside 
(see  Fig.  A.l). 
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where  2e^2  ■  y^  and  1»  2  are  the  ortho  tropic  directions. 

Since  one  is  interested  in  relating  the  stresses  to  the  strains 
in  the  xy  frame  of  axes,  the  usual  transformation  relation  for  second 
order  tensors  are  employed  (see  Ref.  35  for  details)  and  the  transformed 
constitutive  equations  (for  the  kth  ply)  become 
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Next,  the  stress-strain  relations  for  the  stiffeners  are 


=  ^ si  (A  -Q) 

9yyst  =  Br  €yy 
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where  Egt  and  Er  denote  the  Young's  moduli  for  stringer  and  ring  material 
respectively.  Note  that  according  to  the  smeared  technique  assumptions, 
stiffeners  do  not  transmit  shear. 

A.  2.3  Stress  and  Moment  Resultants 

Instead  of  dealing  with  stresses,  it  is  more  convenient  in  thin 
shell  and  plate  theory  to  deal  with  integrated  stresses.  This  leads  to  the 
introduction  and  definition  of  stress  (N^)  and  moment  (M^)  resultants. 

For  a  stiffened  laminate  these  are 
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where 

l  and  l  are  the  stringer  and  ring  spacings  (respectively),  A.  denotes  the 
x  y  1 

proper  stiffener  cross-sectional  area  with  Ax  denoting  stringer  area  and 

A  ring  area,  and  h  h  denote  the  outer  surface  and  inner  surface  coordinate 
y  on 

of  the  laminate  (see  Fig.  A. 1).  Note  also  that  the  above  definitions  lead  to 
the  sign  convention  shown  on  Fig.  A« 2 

Substitution  of  Eqs.  a-5  and  A»8  for  the  stresses  in  Eqs.  a-9  and  a-10 
prior  substitution  of  Eqs.  A-l  for  the  strains  in  Eqs.  a-5  and  a-8  and 
performing  some  minor  mathematical  operations  lead  to 
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positive  z)  and  I  ,  I  are  the  stiffener  second  moment  of  areas  about 

X  y 

c  J  c 

centroidal  axes. 

After  performing  the  indicated  operation  [Eqs.  A_n  and  a-12]»  one 
may  write 
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0 


“ff(  ZyjQAy)  O 


0 


(A- He) 
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04-/$-; 


with 


O') 


B»  = 


2  Q.  (kk-D 

k  =  l  0 

«  -  W  .  a  ,  * 

2  Q*  (A«-A*.,) 


Di*  = 


*  -(k)  , 1  v 

2  Qii  (A. -A*.,) 


Since,  in  the  derivation  of  the  field  equation  for  this  formulation 
the  dependent  variable  are,  w  and  a  stress  function  F  (through  which  the 
stress  resultants  are  derived),  then  it  is  convenient  to  express  the 


mement  resultants  in  terms  of  the  N^'s  and  the  h  's. 


Starting  with  Eqs .  A-13 ,  one  may  write 


'N„ 

Nyy 

=  tfq] 

-  ts^i 

'k„ 

kw 

Nyy. 

M. 

2K.) 

(.A-/6) 


From  this,  one  can  solve  for  the  strain  vector,  or 


i 

T> 

ICO 

-t- 

Xxx 

=  lAj" 

Nyy 

Kyy 

.Nyy. 

2  Hx)/i 

(.A-/7) 


Another  form  for  this  equation,  Eq.  (17),  is  the  following 


where 


IV 1 

'n„ 

jCxx 

• 

€yy 

ii 

*5 

Nyy 

+  i%i 

Kyy 

Vxy 
**  J 

.***. 

(4-/8) 


[%]  =  i  A:ii 


l6it]  -  [A^Bij 


(A -IQ) 


r 


Next,  substitution  of  Eqs .  A-18  into  the  expression  for  the  moment 
resultants,  Eqs.  A-13  ,  yields 


*Mxx 

'e* 

Xxv 

Myy 

=  [  B.j.] 

-  ID.,1 

Xu 

Mxy 
•  « 

ri 

» 

2X*y 

[§#%] 

Nut 

Nyy 

X*K 

Xyy 

.N«y, 

T 

Nxy 

4 

Xxy 

Nyy 

Xyy 

<  '  ; 

2Xy 

(A-20) 


where 


idj  =  [%i  [^i  - 


(4-2/; 


Note  that  [a^]  and  [d^]  are  symmetric  three  by  three  matrices,  while 
[b^j]  is  a  nonsymmetric  three  by  three  matrix. 

A .  2  A  Equilibrium  Equations 

The  equilibrium  equations  are  derived  by  employing  the  principle 
of  the  stationary  value  of  the  total  potential. 

According  to  the  principle,  for  equilibrium 

5  Ut  =  0  (4-22; 

where 

Ur  =  U.  +  Up  (4-23) 

the  sum  of  the  strain  energy  and  the  potential  of  the  external  forces. 
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From  Eq.  A-22  one  may  write 


SUt  =  +  S  Up  =0 

3**  L 

“  f  ( Nxy  &  €*x  +  Myy  &  €yy  +  A//y  §  ^£y 

~  Mxk  S  ~  Myy  &  Xyy  ~  2  Mxy  S  Xxy )  dXd )/ 

M  i  HR 

-f  ffts*dxdy  -  JT  MutSM  +  MySlT 

+  QxS<V  ~  ~  tfty  &  dy  (A~24-J 

Where  q  denote  the  external  pressure  (positive  in  the  positive  z- 

direction)  and  the  "bar"  quantities  denote  external  loads  applied  at  the 

boundaries  (N  and  N  are  in-plane  loads,  while  Q  is  applied  transverse 
xx  xy  x 

shear  load  and  M  and  M  external  moments).  Note  that  M  and  M  could 
xx  xy  '  xx  xy 

represent  moments  arising  from  eccentrically  applied  N  and  N 

xx  xy 

Use  of  Eqs.  A-2  and  A-3  for  expressing  the  variations,  in  the  reference 
surface  strains  and  changes  of  curvature  and  torsion  in  terms  of  variations 
in  displacement  components  yields 

SU  T  =  m  N  „(SW,x+w.,Sw,  +  w:*Sw.) 

+  Ny,  (StCy  -  +  W./SW.y  +  W.yJW.y} 

4  Nxy[8u>y  +  $  ^  'V-y^vV»X 

+  S  V\4y  4  lV,y  J  “  Mxx  $  M,XK  -  MyySW,yy 

-2MxySw,xy}dxdy  -  J  (frS«/dxdy 

nt.  _  °  _  _ 

(  -  A/xy  £  W  +Ny$ir  +  QyS'lA/  -Mxy^x 

-  PixySVyliyy  (A'2S-> 
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Re-writing  the  above  in  a  convenient  form  in  order  to  use  Green's 


theorem,  one  may  write 

nl  r 

•  +  Nxx  (  +  »V.x )  5 W 

■f  Nxy(V\/,y  +  W,°y)  5 IV  "  Mxx$W,x'),x 

+  (NyySlT+  Nyy(  W,y  +  4 

4  A/yy  (  W,x  +  Kx ^  5W  —  MyyS  lV,yXy 

-  [NnjSU.  4  [/Vxx  (w,x 

4  Alxy,x  +  £  ^xy  ( ** 

~  M*x,xSw  +  MyyySlf  i  C  NyyiW.y 

+  Wj)J,ySvV  4  N*y,y$U 

4  [  Nxy  (W,x  4  W,J  )]  >y§  W  -  Myy,y  $W,/) 

-Jh.sw  -  2Mv$w,ty]dxd)t 
K 

-XT  ftswolxdy  -fj-N„SU  i  Nty&ir 

iQySW-  M ,*$90,  -  M,fS PjjVy 


NxxSU  +  N«(W,»+w,»')JW-f  N<ySV  +  A/*y(W,y 

+  W,y)  +  Mkx.x^W  +  2Alx/,y<5wXx 

+  (Nyy§2T +  N,yt^y  +  'V!J;Sw  +  NxySU 
+  NxyLKx  +  Mx )  W  "  Mw 3 W,y  4  MyyjS*/ 

4  2M,,,t&w),y  -(Nxx.,SU  +  (A/w(W,x  +  W,«)},yi5V 

+  Nxy.xSV  +  t  Nxy  (W,y  4^)1, X  iW 

4MxxMSW  +  My„yS2r  +  rA/y,(W,y+W.*y)),y^ 

+  N,y.y5M4  [A/.y(Wx«  +*/.x)],g<W 

+  Myy,yy£w]  '  ^5lV 
-  ZMx^ySwjdxdy 

_  '/0  [ - Rtx^d  +  Tlxyiir 

4  &«W  -  a,8%  -  MvSVyj[dy  M" 

By  Green's  theorem,  one  obtains  the  following  equilibrium  equations  and 
associated  boundary  terms, 

Equilibrium  Equations 

Nw,*  "f*  Nxy,y  ““  0 
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(A-21) 


Alxy.X  +  AJyy.y  =  0 

h^xx,**  4  2  Mvy.xy  4  Alyy.yj/  +^p,  4  N«  ( */■«*  +  w,°* ) 
4  2  Nxy(Vs/,xy  iW,y)  4  Nfy(vJ,yy+  W.,y)  4$  -  0 

Boundary  Terms 
either 


or 


Nxx  -  ~  Nxx 

5U  =  o 

AvJxy  —  As Ixy 

S*/  =  0 

A/*x(^»x  4  */*  )  4  A4xy  (W.y  4'A/,y  ) 

+  MX«,<  +2M«y,y  =  Q»  + 

<5V/  -  o 

M»x  =  Mx, 

Sw/.x  =  o 

u~2?; 

The  first  two  equilibrium  equations,  Eqs .  a- 27  can  be  identically 
satisfied  through  the  introduction  of  the  following  stress  function 


A/xx  ~~  xy  ~~  A/ 


Aiyy  “  '  1 


AX 


»XX 


Nxy  =  -F,xy  fAlyy 

With  the  introduction  of  the  stress  function,  F,  the  third 
equilibrium  equation  becomes 

Aixx,*x  4  2  M xy.ty  4  Alyy.yy  4  h,xx  4  Fyy  (•**/> <*  4  W,aa  ) 

4  RJ0((W,yy4  W>w)  -2F '^CWiy  Wjn/j"  Aj*x(  A/,xy4  VV/,*y) 

4  2  Al*y(w/,xy4  VV,vJ  +  £  =  0  (/1-30J 

A.  2.5  Compatibility  Equation 

Since  the  in-plane  equilibrium  equations  are  identically  satisfied  with 
the  introduction  of  the  Airy  stress  function,  F,  then  the  governing  equations 
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consist  of  the  transverse  equilibrium  equation,  Eq.  A-30  and  one  more.  This 
one  mor^  results  from  requiring  compatibility  of  the  in-plane  displacement 
components  u  and  v.  From  Eqs.  a-2  one  obtains 


(Elxyi) yy  =  2^>*yy  2  W.xyy  ( **4*  +  2  N,*)  -f  y  W,*(  Vjttiy+2.  Wttyy) 

€-yy,xx  =  “  "p>~  +  2  W,y  VN/.yx*  +  2  iV,yw  Vs/,y  +  2  My  Myjf, ) 

0 

^fxy^y  ~  2^>xyy  M.jyy  vV,y  +  V\Z,x  W.xyy  +  MXyy  >A/.)f 

■+  lA/,yW,x*y  +  W,xxy  \A/,y  +  lA/,x  \A/>xyy  -3  /  ^ 

Elimination  of  u  and  v  leads  to  the  following  compatibility  equation 

£*x,/y  +  £yy,«  ~  fcy.xy  =  ~  ^  W,*y  ( M*y  f  2Mx.  ) 

—  2"  VV,j<*  (  Wjyy  2W,yy)  —  ]jrlA/,y y(  "h  2  (<4 

Substitution  of  Eqs.  A-18  [Eqs.  A-29  for  the  N's  and  Eqs.  A-3  for  the 
k's]  into  the  compatibility  equation,  Eq.  A-32,  yields 

fin  F yyyy  +  a,  F,  «yy  ~  ^?/3  F.yyyy  ^  ^//  Mxx/y  +  ■&,  ^/,ryyy  ^  2.&a]*/>xyyy 

*f  Q^ilF ,1 «yy  *f  GiiF,yxn  G&  W.nn  4  vj.ttyy  +  2 $,**/.**•</ 

~  (Xti  F ,xyyy  ~ Q.i$F xxxy  1 0.H  F*xxyy  ~  IA/,x*xy  ”  $-}i  V\/,xyyy  “.2  W**«ry 
S  ~  +  V\4xy(M«y+^  M^y)  ”2"  M><x(M/y^-2VA/,yy)— jV,yy(W,«/2  Mxx) 

<4  "33J 

Similarly,  substitution  of  Eqs.  A-19  into  the  transverse  equilibrium 
equation,  Eq.  A-30.  yields 
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•Qrw  F,gxyy  +  Sn  F ,xxxx  ~  ~Q")t  F *x*y  Ctn  W,x*xx  ~t~  C(,i  *2clts  \hl,xjuy 

+2  8*  F  .x/yy  +2&ai  F,***y  -2  Sxi  Fwy  +2cliiM  «v  +  2d.n\t^>*yyy  ■/  4  cluVJ 
+  Sli  F,  yyyy  i  $11  F,**yy  -  F.xvyy  +  ciyt  W,**/  +  dn  ^J-my  +-?  cLi  W,*»y 

+  ~  f >xx  +  Fyy  (  W»x*  +  W.\*  )  -  Nxx  (  ^>*X  +  \A/»°**  ) 

+  2  A/xy  (W»xy+  U\4*>)  ~2F,  xyL  ^>*y  *  V</>*/  )  ^  F,xx  (W,yy  j-  W/.yyJ 

+  #  =  o  (/i  -i«; 

A .  2.6  Boundary  Conditions 

The  boundary  conditions,  Eqs .  A-28,  can  be  designated  according  to 
transverse  one  (simply  supported,  clamped,  free)  and  in-plane  ones.  Since 
all  of  the  application  to  be  considered  deal  with  supported  boundaries, 
only  simply  supported  (ss-i;  i  ■  1,  2,  3,  4)  and  clamped  (cc-i)  boundary 
conditions  are  listed.  These  are  (at  x  ■  0,  L  )1 


SS-/  •• 

W  *  0  ;  Mm  =  Fly*  ;  N*x  =  '  A/M  ; 

A lyy  ~  ~Flxy 

SS  -2 : 

W  —  0  >  Mxx  -  ;  2^  “  Const .  , 

Nyy  =  Nyy 

SS-3: 

W  =  0  ;  Mxx  -  Mxx  ;  Nx*  ~  ~  Fix*  > 

'll  -  Const. 

SS"  4: 

W  =  o  i  M*x -Mix-,  U  -Coks(  i 

If  -  Ccnsl.  (f\-25) 

and 
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CC-I  : 

vv  =  0  >  Wtx  =  0  > 

A/«  =  ~  A/xx  ■> 

£ 

ii 

CC  -2  • 

w  =  0  ;  V\A*  =  0 ; 

^  =  Consi  ; 

£ 

ii 

CC-3 : 

£ 

n 

o 

^  • 

* 

a 

o 

A/yy  -  ”  A Ikx  > 

2/  -  Consl, 

CC-4: 

;  1A /<*  "  0; 

21-  Coas{,  ; 

XJ  -  Consi.  (A‘2&) 

The  above  boundary  conditions  may  be  written  in  terms  of  the  dependent 
variables  F,  and  w.  The  kinematic  conditions  u  =  const  and  v  *  const  are 
first  expressed  in  terms  of  equivalent  conditions.  This  is  shown  below  for 
each  of  the  relevant  conditions  separately. 

Note,  first  that  the  expressions  for  the  M^'s  and  N^'s  are  given  by 
Eqs .  A_20  and  a-29- 
SS-1:  W-  0 

&tF.xx  +  dn  +  2cL*VJ,ty  = 

J~,yy  ~  0  0»d  F,«y=0  (4-37) 

SS^2:  W=0  __ 

& .  F,„ + A,  + do  w,v  = ^ 

F,yy=0  Qnd  F>xy-0  C/1'38; 

The  u  =  const,  condition  is  expressed  in  terms  of  an  equivalent  condition 
by  employing  the  following  steps. 

The  expressions  for  y  from  the  kinematic  relations,  Eqs.  A-2,  and 

xy 

from  the  constitutive  equations,  Eqs.  A-18,  are  first  equated  to  each  other, 
or  .  .  « 


Jf£  =  VU+W*******} 

-  Gn(F,yy-Nxy)  +  fii}(N*y-FjXy) 

■f  Sji W,xx  +  -&n\ bjjiy  +  2  $-33  iA/>*y 


(/\-38a; 
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0 


0 


0 


0 


One  differentlon  with  respect  to  y  and  use  of  the  conditions  w  *  0  and  F,  *  0 

xy 

yields  at  x  ■  0,  L 

\j,*i  +  k/,xy  +  tV,*v\/,)fy  a  Fyyy  F  Qa  F,x»y  F  Sit 

F  (A-38A; 

Similarly, 

Eyy  -  If,/  -  7-r  +  j[  wi  +  VS4y  IV, y 

=  (  Kyy  -  A/**  )  +  A,i  Fxx  f  ^3  (/7*y  -  F,*y 

+  &i  IV, x*  +  -&a  iV.yy  +  -?  ^13  lA/,<y  (/l  "39£J 

from  which  one  differentiation  with  respect  to  x  yields 

V,n  ~  ^  VV,xy  Vs4y  =  F <xxx  ~  &*S  F ,xty  ~h  zS’a  iV.yxy 

+  Fki  \Mxxx  F  •£«  W,*yy  (  A  -  3Ql>) 

Elimination  of  v,  and  v,  from  Eqa.  (A- 38)  and  (A~39)yields  the  equivalent 
xy  yx 

(to  u  «  const)  boundary  term,  which  is: 

fl'iFiW  F  2  CXil  F,XKy  ~  (XtiF,w r  "  'V'yy  \fJ.xxx 

+  (  ’Qi  ~  2  On)  V^J, «y  ”  (■  ^1*  ”2  &si)  *V»*yy  —  0 
Note  that  because  F,  “  0  for  this  boundary  condition,  the  te^  con¬ 
taining  F,  has  been  dropped, 
xy 

Thus,  for  SS-2  the  final  form  of  the  boundary  terms  becomes 

\N-  O 

4Ft\  F.yy  +  S'hFm  +  W,x*  F  2d/sW>*y  stfxx  +  S',,  Nxx~-&3iNxy 

F,  xy  =  0 

(Xi3  Ftyyy  F  2  Qii  FiXXi  ~  &**  Fxxx  ~  "fiF  **  W,x  W,yy  -  ^hVJ.xkx 
+  (&»-2&i)\fi/,**y  -(&*  -2&i>)W'**y  a  0  (A-4a; 

SS-3 

vV  =  0 

■^ii  F,  yy  F  Fixx  F  da  Vj.m  F  2  d/s  W»«y  "^siKacy  -  Hx*  ^  A/xx  Nt/ 

F,yy  =■  0  \r  =.  Con  si. . 


10 
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Similarly,  as  in  the  case  of  SS-2  (u  -  const),  an  equivalent  condition 


is  obtained  for  v  *  const.  From  Eq.A-39a  ,  since  w  =  0  and  v,y  ■  0,  then 
the  equivalent  condition  becomes  £yy  -  Q  or 

“#/i  Nyy4  Ga  F.xx  f  (Mry  “  F.xy)  +  St,  \Ntn  4  2SaW,xy  ~  0 

Thus  for  SS-3  the  final  form  of  the  boundary  term  becomes 

W-0 

■Oxi  F.XX  d/i  lA/,xx  4  ^.d, 3  lA/.*y  “  Si t  F.xy  =  Mxx  4  Sn  Nxx  ~Su  Nxy 

F>yy  =  0 

GuF.xy  “  GaF.xy  4  SjiM.xa  4  2  SaW.xy  ~  G/z  Nxx  da  NxyiA'^O 

SS-4 

For  this  case  the  equivalent  set  of  the  boundary  terms  becomes 

W  =0 

QnF.xy  4  Si,  F .xx  ~  Fxy  4cLi  W.jot  42cl,$  VA/ixy ss  Mx<  +  $//  Nxx~SnNxy 

QuF .XX  +Q,1  F.xy  -  CXa  F.xy  -f  Si  1  W.XX  +2  St}V^.xy  -  GtiNxx  "Gts  Nxy 
G-l3F,yyy42QnF,xty-C(Xi  +Ga)  F,xyy  ~  (Xu  FxX>  ~-W|*  -  Wj»  W %y 
4  (2  Si*)  W/>jfyy  4  (Sji  \yj,xty  ~Sz\  \A/.xxx  ~  0  ~4-2) 

Following  similar  steps,  boundary  conditions  CC-i,  i  -  1,  2,  3  and  4,  are  also 
expressed  in  terms  of  w  and  F  only,  or 

CC-1 

W=  VV>X  =  F,yy  =  F.xy  =  0  {A  -4-3) 

CC-2 

W  =  W/>x  =  F.xy  =  0 

QtiF.yyy  42G»iF,j«y  ~  Gu  F.xxx  ~  Sn  W.xxx  4(Sit-2  Sa)  ^J.xxy  -  0 

CC-3 

14/  =  VV.  x  =  F.yy  -  0 

Gii.F>xx  ~G.iiF.xy  4  Sn  VV.xx  =  “  QiiNx y  ^ 

rr-4 

W  -  W,*  =  0 

GuF >yy  4  GuF,xx~  GaF.xf  4  Su  \A/,xx  ~  GuN»~G ijNjy 

Gi^F  *yyy  4  2  G.i3  F ,**y  ~  (Giz4  Qu)F ,xyy  -  GuF,rxx  ~  Qi  1  W/»xm  Q.,)  K/*wy 

=  0 
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II.  2.7  Solution  Methodology  -  Field  Equation 


The  solution  methodology  is  an  improvement  and  modification  of  the  one 
employed  and  outlined  in  Refs.  36  and  37. 

The  separated  form,  shown  below,  is  used  for  the  two  dependent  variables 
w(x,  y)  and  F(x,  y) . 

F(X,  y  )  =  Co«H-  ^  [ C;(X)Q>4  l-$  +  2);  (X)Sin-^  ) 
w(x>y) =  +  8$)  Sin-Tp ) 

where  n  denotes  the  circumferential  wave  number. 

In  addition,  similar  expression  can  be  employed 
for  the  imperfection  parameter  w°(x,  y)  and  the  external  pressure  q  (x,  y) . 
Note  that  in  most  applications  the  pressure  is  assumed  uniform  (<jq  only). 

w°u.v)  =  4’w  +  ^(AlwCosi?  +  ) 

#cx.y)~  $.w  +.£(&' (*)£,#  +  £*«  S*!?]  (A-48) 

Because  of  the  nonlinearity  of  the  field  equations,  Gqs.  a-33  and  a-34 

substitution  of  Eqs.  A-47  ,  and  A-48  into  them  yields  double  summations  for 

the  trigonometric  functions.  These  double  summations  involve  products  of 

sine  and  cosine  of  iny/R  in  all  four  possible  combinations  (cosine-cosine, 

sine-cosine,  cosine-sine,  and  sine-sine).  Furthermore,  these  products  have 

different  origins.  Some  of  them  come  from  products  of  W,^  ^’xy  ot*iers  ^rom 

products  of  F  W,  [see  Eqs.  A-33  and  a-341.  In  order  to  simplify  the 
,  xx  yy 

final  expressions  (and  use  single  sums  instead  of  double  sums),  and  in  order 
to  cover  all  possible  combinations  of  double  sums,  the  following  simplifying 
equations  are  presented.  These  are  based  on  trigonometric  identifies  in¬ 
volving  products. 
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£  Z  (<fbiCosj0]  Ch  iCosid  —  Jl  Autuo  Cc$iO 

i-  oj--0  *  £s° 

(Xil  CoSi-6  =  Z  Al’tTZ(ic)(b>&)&n.LO 
tbr-0  s  L--0 

j^£L[j  kjSinJdjG U  Sin  iff  *  Z^Amcx)(b,&)G>siff 


£_Zl[  b;Cosj& ]  (Aib  C*sL0 

Vo<j-o  & 


KK  .  I  .  ,  „ 

Z.  Ani(j<)(b>R)C&lL6 

i-0 


-  ZL  A 133  w  ( b>  (X)  Sm  i-8 

L-0 


(A-4-9) 


6os^e]  GUI  SiniB  ^Ain(jc)ibA)S'ini8 

1=0  « 

£.ZZ[biSinA8]  GliLCosIO  =  ZS  Anzuv{b,a.)SinL8 
W'0  i  m 

HCb^Sin^ausinie  =|:  Ai^ib^Cosie  laso) 


.  I  /C+i  4  t  ^ 

^  Aji\i£b>&)Cosi& 
/t,L.(t*bjCosfi62  (\iSini6  -  i 
£zqJi>S;*iMfliCos;e  =  ^ 

i.-0  j-0  ° 

y ZL U'^SinMaiSmB  - ^  AJ2UVlb,G.)Cosi8  c-A-^O 

t-O^O  a  y  t 


where 

Aijnx)(b>&) 


ji  CU+j)  bty  +( f-  fatfaw 

j--o 
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/  ^ 

Aij2iK)(b>A)  -  2.  bty  4  )H-(}I b/iyj<t  flj- 

A  +  (-!  +  bu^&Cli 

^  I  ^ 

Auiw(b>&)  -  2  bii-ji]jQ-i  (A-52.) 

J'® 

L  I  ^ 

Ai2UK)  (b,&)  =  ZT^C  k>L+j  +  (  (Xfi 

A  iim  (M  =  + o-fafabn-j  Jfat 

A  (b,d)  -  -2  + 

/I  I24(*J  (  ^  ft)  =  J  ^-e t  bi+j.  f  (  -  I  -  ^  ^  )  bli-j  I  ]<f  ^  (  /)  -Si) 

/4  J2|W  C^ft)  “  ^ +  (/'  <% 

d 

(i,aj  =  f  £.  (u+j/Ait.j  f  (-"^+£ )  at^A^.] 

Aj-MItj  (&•&)  ~  ~2  l>i$*  (~  l~%itfi)(ty)  ba-jiJQj.(/\-54) 

K>-L.  f  - 
>  ^ 


and 


o 

/ 


'/  ^  ^  o 

M  L--° 
it  l  <o 


Next,  returning  to  the  solution  procedure,  the  expressions  for  F,  w,  w°  and 
q,  Eqs  A-47  and  A-48  are  substituted  into  the  equilibrium  and  compatibility 
Equations,  Eqs  A-33  and  A-34.  This  substitution  yields  the  following  non¬ 
linear  differential  equations: 

Equilibrium  equation 

Zf/W Gs%  +  8i,u.  S;*%)  Sin % 

L=o r 

f  BiMiCos  %)  -  h*  (-#)'  U,.„Gs  W  +&,»  Snty) 

+  h  13 L-$f( A  i.xS:n^  -  Bi.x Cosig )  -tfuo^flAiC^Tr  +  8, S.rvf)] 

+JLl$jCi,mc°sig  +t>ijm,Si*yt)  +&,(%)(- c^Sintg 

+  Ccs!%  +  Vi.uS1n.7ir) 

+  g^frO.xSn^  -  ZV>  6>s%y)  +  &,  (#)4(Ci&^  +2?,SnWjl 

(Ci^Cv-f  +  -D^SmW)  +  i  ( F.ww’l 

'  £=6  ' 

"  Nxx^f^i^X +(6i,Xx  +  Bi!j0c  )5)n^J 

f-2  Nxy £(%)(’  (Ai, X  +Ai,K)$inr$  +  (3i,x+  Bl,\)Co%^] 

+  +£s«#J  =  0  (4-«; 

where 


L(F,w+wc)  =  r £  C;.„  C.s^  +Z>,„, 

-CBi+ -2(£(%) (-Ct.xSm%  4VuC~&)][£ 

+  (Su  +  8,’,;  Cesty}]  i£-(% f(CcCoS% 
■f  Vi  (Ai.Kxt  Ai.u)  4  (8,-.w  +  Bijx)Sm^  jj  ”  0  ( ^ 
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’I 


HP.  w+W)  =  (%izz  C-iUj+Ai J  Ci.u G>s 
+  ( *•)  ( Al  C°5^J  T>iMS;n!% 

4 (b'i^N‘^i8;)&n^J  C..»G,s% 

+  (%)££  +%)&n*g)  Vu,S:n% 

-  Ccs^JCU  Cos^ 

-  (|;|z  CM;.xx  +A^)C^J 

-  tsps.igjc^c*^ 

-  (fix £  C(  8,„ 

L-04-0  ° 

C,.>  i  Sin% 

42($/ff  [^4,  •/^Si^^jDi.xiCosf 
i^oj^o 

+2(|j  XX^^.* + Sj^Cos^JCi.x  i  &7nr 

4  2L%)£z  [ -i(Bj*  tq^Catyfo*  IC*% 

izOj  = 0  w 


and 


<• 


h<to  -  dn 

hn-  20(31  +2d,3 

hu-  d iii 40(33  fdi 
h.13  ~  2  da  ^  2  d*3 

h*<i  -  dn 


fro*  &. 

Qn  -  ~~-B-i\’k 2b*} 

Qh  =  ■&»  ~2&t j 

=  2  ^/j  -  -Q-ii 

~  'Gm 


<» 
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C  A-Stb) 


J 


Note  that  the  operator  L(F,  w  +  w°),  Eqs  A-56,  can  be  written  in  terms  of  a 
single  series,  which  is  the  most  appropriate  form,  for  use  in  Eq.  A-55  .  This 
is  accomplished  through  the  use  of  Eqs.  A-49-A-54. 


+  (/U+A*,C  )  +  A 

+2/1 L*  (A,m:  .  C,x) +2  A 

-Ip‘£  (a: _ 


+/4L«,  -M  J 


-2/lijj,«j(/l.»+^°*  ,Z5.x)-2Auj(»/B'«+8.»  .C.>] S:n^ 


CA-SV 


+|eW«(Ci,»«,C«^  +Vt,„aStn^) 

-7 XwC*#)  +<.-2Qn  +  aa)($)‘(-ClMSx%  -  L>iaS;„f) 


€ 


+20,» <#)’(-  G.„S M%  *  T>,,c,s%)  +  a„^)\CiCcs  f  +1 

L~f> 

i  Amm  [&+*&'•  B>JCx)-f  /4»aJAo+-M.«»4)  iAi34»)(B,^2S,xx,8) 

l  o  .  .  t  _  o  ^  i^y 

4  J2 y\ u4w  (.A>xHA,x  >  A>*)  4,2  /4  IJ/CW  S„)JGtf  7? 

-  +  /I  L<«  LBi2B\A,») 

~X-  L- o 

+  ZlLw  (/4,»+i/4“x. ,  8)  +  /4iu(w(S.»+-?Aw,/4J 


=  0 


C/1-58; 


Parenthesis 

As  far  as  the  equilibrium  equation  is  concerned,  the  summation  starts 
from  zero  and  goes  up  to  3k  [see  Eqs  a-55  and  A-57J  because  of  the  nonlinea¬ 
rity.  The  Galerkin  procedure  will  be  employed  for  this  equation  in  the  cir¬ 
cumferential  direction.  This  will  yield  (2k  +  1)  nonlinear  ordinary  differen¬ 
tial  equations  [from  the  vanishing  of  (2k  +  1)  Galerkin  integrals]- 

On  the  other  hand  the  compatibility  equation,  Eq.  A-58,  is  written  in 
series  form,  from,  zero  to  2k  Because  of  the  orthogonality  of  the  trigono¬ 
metric  functions  (4k  +  1)  nonlinear  differential  equations  result,  which  re¬ 
late  the  C's  and  D's  to  the  A's  and  B's  [see  Eqs  a-473-  This  set  of  ordinary 
differential  equations  is  shown  in  a  complete  form  in  the  pages  that  follow. 
Before  showing  them,  though,  some  simplification  can  be  made. 
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For  the  case  of  i  ■  0,  one  obtains  the  following  equation,  from  the 
combatibility  equation,  Eq  A-58. 


+  f  Q.u  Co.mn  ~  2  (If)  (  Amw  (A+^A  >A,*x) 
*f  Aj24(fi)  (S^^B  >3»xx)  ^  A  I2KXJ  (  At**+-2A)KX  >A) 

+  A12HW  (B,xx+-28,yx ,3)  +2  Aiuik)LA>xi2A*>A,x) 


+  2  A  mot)  ( B,x+i  &,x ,  B,x)  J  ~  0 


CojMxx  =■  aaJ  { "  Quo  A  -  -fcAcu  +  Jr  (r)  )  4)'*x 

i  j*(  )  %x*  Gfyxx  t2  Aj,tx)Aj  -f  d  -+2  Sj.xx)  8j 

i  Al'X  )  AJlX  (Bj.*  128#, )  8,-jj  (/t -s<?; 

Moreover,  the  displacement  component  v(x,  y)  is  a  continuous  and  single-valued 


function  of  y  (and  x) ,  therefore 


rite 

I  u}dy  =  ucx.W)  -  mx,o)  =  o 


From  the  second  of  Eqs  A- 2  one  may  write 


V,y  -  e‘„  +  w//;  -  W.,  (W.y  12M*)/z 


(/l -6o) 


{A -60 


Furthermore,  use  of  Eqs  A-18  [relation  between  C^and  ,  k  0f  Eqs  a-29 
[definition  of  stress  resultant  function],  end  of  Eqs  a. 47  end  A-48  [assumed 
form  for  W,  F  and  w  ]  yields  the  following  relation, 

r2\R  __  _ 

J  V,,dy  =  J,  C-OaNa  +Cl.>  M^Uiy 

^  Jo  C  G.u  F*yy  ~b  F >  xx  ~  Q-i}F>xy  4  xs3$  W,** 

4  -On i  lA/yy  4  2  &ii  \A /tXy  4  ~J=T  -  ^  lA/,y  ( iV.y 

+i*/',)]dy=  0  ( A-62 ) 


tXtiSlk  jJC 

J.  (-a«N «+aaNv)dy  +J.  [a*  £(%fc-c>c*% 

-ViSin^ J  +ailjr.[GMGs#  +  DiMSin%] 

-  ^f^r-GxS*#  +  axC-.f  J 

+  &,  Z-o(Au*&!, f  fBuxSin^J  +  £  c4?/f 

-AiGs-f-  Bi  Sinf]  +3  &»l(%)(-A;.,Si*W  +8uC*% J 

^7?f-(AiCos^  +  8iSin$)  --^Z  [fyC-AiSinyj; 

Lz0  iso 

+  6;Os^J  (-  (Aj  *2Ak)  Sin  %  H6ji  2  B’>  C*<?  jj 

=  0  (/I -6  3) 
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This  equation,  Eq.  A-63,  after  performing  the  indicated  operations  (integra¬ 
tion,  becomes 


+  (8jtiB‘)Bj]}dy=  0  (A-M) 


From  which,  one  may  write 


Q»=  if-  -  4?  +£>  f/UA,^)A, 

+  (^-f  *?Bj  )  SjJ  "f  Nxx  ~  d»2  NxyJ 


The  remaining  compatibility  (nonlinear,  ordinary  differential)  equations  are 
For  1  -  1,  2,  ■ ■  2k  and  cosine  terms 

fl.uCi,XXKX  ~  2dn  (^jZ^i.JOfX  “  Ci/XX  Dt/x 

1  CluC^f  Ci  +  S'lC  Ai.xxxx  ^  @3\  (%)  &i,m  ~$u 

+  &»($)%  +  4f  -  r(f/^i+^A‘vU 
-(Qfi-J Ui+ihAui+SAli)  +  <■*  ’  1i-\  Hi-jHAa-j i  + 


+  Au.jha  +  2  Aa-iijJli 

4  f  Bi+j,XX  +2  Bitj.tx  ~  %.j  (  &/iyJ,*X+2  &li-jl,*x))j 

■4  2  C  ^+j)  ( A  ty,  t  +  2  Ai+j ,  x)  -  ^  /i-ji  ( A/i.j  i,  x  4  2  A ti-ji.xi] 
j  Aj.x  +2(u+j)  ( 8^, <42  8^,*)  +  (-2 -%-L)li’j 1  ( Bli-ji'X 
4  2  B^<j  =  0  M~£6) 


For  1=1,  2 .  2k  and  Sine  terms 

QnDi^+^0.^(J^)Ci.m  -  (2<2/W  &*)($)  Di.xx 


-  2(Xk i  ($)  C,x  4  (^f^L  +  $l(@40  &i, 


xxxy 


-fa  $  Ml*-  fa  (%f  Bit„  4  UVfAu  +$•«  (9f ft 


+  -  X  (f)  t  ~  4  (Ay+tfq) 

+  ft*  ( W/(  fi'^v' B,lj,)]  Aj.xx  4  f-  (/^.^M^x*) 


4  (^“^j-iJC^^i.xx  42  4  ((Buj.t*  42  Bi^.ix) 

C  8/i-^i  ,**  42  B/i-jifx)]j  A+  ~-2  )  M i/j.x)  4  frjti'jl 

'(Aa-ji,x  •+2>4/i-ji,x)]<j’  ^/.x  ~2C-(iij)(Sit(i,x-t^Bi^tg)  +  (J-^.i)U~<il 

•  (B/i-jU  =CI  (4-67) 
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where 


0  i>  K 


l- I  0 


£  <  0 
l=o 


L  /  i  -  K  l  /  £  >  0 

As  already  mentioned,  the  Galerkin  procedure  is  employed  in  connection 
with  the  equilibrium  equation,  Eq.  A-54,  in  the  circumferential  direction. 

The  vanishing  of  the  (2k  +1)  Galerkin  integrals  yields  the  following  set 
of  nonlinear  ordinary  differential  equations. 

For  i  »  0 

hloAo.KKXX  4  Qo.XXXX  4  7?  (A*,tX  1 A o.tx)  Nxx 

-  (jz)  4(j  (fyl  ij  (Aj.** 

Q 

4  4)'^)  Cj  4<^(  Bj’Xtj'Dj  42(f 44^* )  Q ,x 
i2&(E>i.x  IB^x) ~  0  (y^-6 8) 


By  employing  Eqs.  A-59  and  A-65  one  obtains 

Ac,mxx(cL ( "^^aT, )  -  '4*'**^ 

4  4 2Aj'**)A± 

4~t 

+2(4^4]'*)^*+  (Bj*  %)%*»+ 

4  2  (Bj,x  42 %<)  %*]  4  4^4]  (Bj  +?*/ )  £jj 


-2  foy  4/4J )C^x  4  ( )  ’Zfyxx  4 2  4 
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+  2  ( Bj.x  +  8j!, )  Vj,*  +  ( 4”  +  A]m)  Q  +  ( ]J 
+  ~  SlkNxy  +  &  =  0  (>4-^9; 

,  2,...  K  (when  the  weighting  function  is  cos 

R 

c/,,4  i,XMX  "/■  4d,U%)  B  L,XXX  ~  (  2dn  4  4  da)  (  ^rj  /4  i  jXX 

-4d* B lx  4  d*x  (4z)  Ai  4  Ci>****  f  (-2  -@ni-  Su )(^)T\U) 

“(■^ii  ~ 2  An  4  8n)(^J  Ci,xx  ~ (2.  &/3  -  $j x)C^)  7 ~)i.x 
+  &ti(Tr)  Ci  +  7^G.*x  -  ($)  i 

-t&li’lUjtiAM  4  (V-JS/j  8,-jJ 

-  (^f(^£)(a,,Nv-duNv)  ~  (Al**+A  L)  A/„ 

4-2  Ajw  (-^)(B;,*  +  6u)  -j(s)  frt{du-j)  Suj(Auji4uj) 

4  C  1-^/j  M/i-ji  >4,Ly» )]  Q,** 

4  C  ^ 1+4  (  &ty  4  Bty)  ~  tyty  (  i-~j)  S/i  jl  (  B/i-j I  4  8/i-j J ) ]  T)j ,  V 
+  .2  C  Ciij.)  Suj  (Auj.x  fAtyx)  ” 

4  Atji,*  )Ji  Q,x  4-  -2.  (dtj)  S fy  ( &i #,x  4-B^.x) 

4  )  tt'jl  Su-ji  (B/i-ji,x  4  Bd-jij) lJ  <f  T^'/X  4  (Sty  ( A  L+j'XX  i  Aty.x*) 

+  (2~thi)  +  A 4 (Sty( 

~  ¥ty  £/tyi  ( B/tyi ,x>  4 B,i<ji,ix)J^ 2D; }  4$L  =  0  (A-7L 


For  i  -  1,  2...k  (when  the  weighting  function  is  sin^p) 

d„  &l,m,-4d,st%)Ai.m-Udn+4du)(£)e>u*  +  4cb&fAi* 

•f  duOfi  8i  +  Aixi  'Dx.hm  ~  (.2  Sn -  Sit)  ( 7?)  C  i,*x* 

-  U^-Axgfc,.*  + 

t  4 

4  (Bj-aBj)Bdh  (■#)’( -  Q»fi,y)- (B;m+B°«)N« 

-  2Hty (.%) (At,*  1-Ail) -  i (fflz.itw'fy ( %  +  > 

+  fyi-itSiy  (8,ivi  4  B/yi)]  C-«4 

+  S'i-jl  (4^,  Myi)]  Xjj.w  "2  t’W)  f  6^-'  +  ^ 

4  (2~?rOl  ^ n’Jl  ^  9 <j  Q-x 

-  2  [ (LW %  (/)4<4^.J  4-fc..  /i-4'  5,iV,  (4 u-jl,x+A°-ji.,)]i  Vj.x 

O  2 

+  C 5 :^  ( Bif^xx  4  ^ ,*x  )  4  ft-j  5“  i m  C  6/i-^ ,**  4  £//-</ i )  J<j  Cf 
t£*0  ^-?0 
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Clearly  the  response  of  the  configuration  is  known  provided  that  one  can 
solve  the  nonlinear  ordinary  differential  equations.  Their  number  is  (6k  4-  2) 
and  the  number  of  unknown  dependent  variables  (functions  of  x)  is  also  (6k  +  2) 
These  are  (k  +  1)  A^'s,  (k)  B^'s,  (2k  4-  1)  CYs  and  (2k)  D^'s.  Note  that  Cq 
can  and  has  been  eliminated,  through  Eqs  a-59  and  a-65  and  therefore  both  the 
number  of  equations  and  number  of  unknowns  is  reduced  by  one  to  (6k  +1).  In 
these  equations  there  is  one  more  undetermined  parameter,  the  wave  number  n. 

This  number  is  determined  by  requiring  the  total  potential  to  be  a  minimum 
at  a  given  level  of  the  load.  In  other  words  the  response  is  obtained  for 
various  n-values  and,  through  comparison  the  true  response  (n-value  and  cor¬ 
responding  values  for  the  dependent  variables)  is  established. 

So  far,  the  partial  differential  equations  are  reduced  to  a  set  of  (6k  4-  1) 
nonlinear  ordinary  differential  equations.  Next,  the  generalised  Newton's 
method  (Ref.  38),  applicable  to  differential  equations  is  used  to  reduce  the 
nonlinear  field  equations  and  boundary  conditions  to  a  sequence  of  linear 
systems.  Iteration  equations  are  derived  by  assuming  that  the  solution  to 
the  nonlinear  set  can  be  achieved  by  small  corrections  to  an  approximate 
solution.  The  small  corrections  or  the  values  of  the  variables  at  the  (m  4-  1) 
step  in  terms  of  the  closely  spaced  state  m,  can  be  obtained  by  solving  the 
linearized  differentiate  equations.  Note  below  the  way  that  a  typical  non¬ 
linear  term  (product  of  X  and  Y)  in  the  differential  equation  is  linearized. 

Xm‘ym‘  =  (X*+d)C)(Y'+dY) 

=  X*Y*+XW  Y  dx  dY« 
*xV+  YdX'+XMY’'+X’,dYl -X'Y* 


=  x"(  Y"+  d  Y")  +  Ybr+dX)  -X”Y 
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*1 

(A- 

where  X  &  Y  can  be  A^  B^,  or 

By  making  use  of  Eqs  (72),  the  linearized  set  of  governing  equations 
(iteration  equation)  is  obtained  from  Eqs  a-66,  A-67  ,  A-69 »  A-71  •  These 
are: 

1.  Compatibility  (i)  [cosine  terms,  Eqs  A-67 ] 


=XY  +YT-X7 


For  i  =  1,  2, _ K 


Q»Cu!n  -2aA%) 

+ z  a  3  +  a,  (%fcr+  sL  f  a  7.L, 

+  (■*&»-  &,) ($) b”L - ( A, f  ■&,) 

- <■  2 A* -&„)  VitB*  *  8,*(%fA7‘  +-nA*:l 

-MUCt Im,  +^r+a4l)/|"  -  4*^4 

-  & 0*£  {  jfiA+W)  A*"  i  jjM  A”  - 

+  ky'  (6  -HB )  Sj,»»  4  A;,"  (8  +-2S”)  Bj- «« -  Ay  (.8i2B)&j.„ 

t  2  C  t"  < »  +^J  4’  +  4  ( 4  f -J/T;  4  *x  -  iy  (»  W')  AjZ  J 
■*  2  C  m7(bub°)  B*<  +  rt”  (6  «s’;  8d”'  ib’j  bJ!* J 

+  AJy  (4+«4"M*  +  NyM+j/O/C’- 

+  0™(B+iB)&J  4  0"(84Jg°)5*'"-0”(84jB;8pJ 
=  0  (/)  - ; 


where 


Y)  =  ( i+<J )  %  >£  +  « '  &  «  v  /&*  y«“, 


L*  CVJ  =  C  ^  i+<j )  Yiy.*  -  j ijit-JiSjiyi  %.;,,*]( j- 

Mij  (Y)  -( u+JiSi+j  Yuj.x  +  (2~7j-i)  a  '<j  i  <r  /r^/ 

Nl(j  (V) ' t  Yty*  ^  c2  ~  7^- 1 )  Sii-ji  Ya-ji.xxJJ 


(il)  [sine  terms,  Eq  A-68] 
For  1  =  1,  2, . . .K 


Q‘i&i,x**x  +  2  (lt)Ci,xxx  ±  Cl)i)  (^)  'Di.** 

~  2fl,)  (j%)  Ci.x  +fl»(%)D*  •+SilYh,BiiX)g0l  ~  (2  YhiSu) 


*  ^J/4cxtx  “  (  ~2-$u+  Sx^C^f  Bi.n  ~H  2  $!*-&*)(%)  Ai.* 

+  Aa^ar*  -  jWX*  +  <tf«  &>aZ 

-  eX,]  -  (£/|J  O^B^S-j  ^.+ 

-  O^S+iS'Mjl  +  ^'M+^'j  B**+  Rxi<AtiA')B™„ 

-  8*x"  -2  c  S^"(8+38°>A^  +  Sy  CB«b'X»' 

-  ^ (6«8')  ^ 3  -2  r  T^rt+aO  S*,  + 

-  t;J(  4+^’;Sj7»  J  +  u*p a’J + u"  caw’)  A? 

-  U^(SU6°)A]  +V^"(.A+2A'’)B"i\/!*(4+M')B’P 

-  uj =  0  (A-74; 
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where 


Qy  (YJ  =  City )  S;t<j  Yitj  ->  Y,i-ji 

( Y)  =  -  Ci+rfj'&j Y*J + 


Si 7(Y>- 
t:J(y)  = 


U?  (Y)  - 
Vlj  CY)  ■ 


[-Crf(t  Yfyx  4  W  -  YY-JlHlSn-ii  Y,*h  J  <f 
[  (l+jl  S;q  Yi.j,x  +/i .jli'jlSlt.jl 


(Sty  %jM  +%t  Sn-ji 


(2^  Equilibrium 


(i)  [i  =  0,  Eq. 


Ac,xm(d„  -  -&L  )  )  ~hin 

4  4 „)  4  ($)%  n  ■&  tAtf*. + (^j)Aj:: 


-A}/$« 4  Aj:\ A;  4 (4?*+2A{„)/r-/\?„$*2$$ 

4  a  (A*  4  -  2  A-,  A”  4  B-"  8*„ 


,  *+l 


0  x  A**l  A«  tmlA* 


4  l8p2Bj)Bp  -  5/Sj*  4(BjZ-h3B-„)Bj 
4  (S^yfSBjY)  Bj"  -  (8j*,4  2Bj.n)B? 

+2  Bp  Bp  *  *lBP*2Bp)Bj?-2&£3 


o  ^  M 


'SsUrtinAPiA-uP-AtAi 


°  .  i  4  >*L 
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V  "  °d  °d 


-  T (t?)  C*xx  +  Mf  +^j )  C**  -  Aj  C* 


-lktt  ,  _>t  _  c  _ m*i  _ »i  _  >/ 

-f  8;  X^,«  +(Bd+B6)  T>j,xx  -  Bj  T>yJ<  J 


*  ,  _,e ,  *»♦» 


+•2  +  -2  (&  +4- )  q*'-J2  4*x  Q* 

+  2  8^  f  -2&V*, 

4  4j;«  c;  4  C  £.  -f4„;  q”"-  ^  c/ 

+  8j?»  'B’J  +  ( +  8/<*0 "Dj  -  Bj^'Dj J 

+&kN»-£iRv+£=  o 

(ii'1  [i  =  1,  2,  ...K;  weighting  function  is  cos  ¥) 

du  Aiacxx  4  4-cbi(-$) Bi.#>  -  (2 cL?+ 4-c(n) ($)  A 

-  4 cug  )C+ d»  (#/vc+ ■«..  g* 


»! 


4(2  (&11  -  Qi\)  (i£)'Di,w  ~  -  2  C  i,> 

-  u&, 3- &„)(%)’  v?:  4  4,  <#cr  +  ^  c*;: 


- ow  lb f  ~ &. ",  -  ^  (w°>az:' 

4  a,  a">c  -  i/cc  -  ^  (^r^A*Ar 

+  aili\aj\4;  j  (A-^ApA-f  a'+Al)  Up'^A'lAi 


MAi+AlHAjWl  )A?-2  (A?+Al)(Aj'-u4j)Af 

+  Mr  u: ) (S-  U 8j )  B-  +  (A?+  A\ )  a* 

+  (4?+A BJ,i-2(£+A‘)(8]+28‘)B”]} 

,  .w»+l  .0  _  _ 

.  t^)(ii ^i)ca,Na-a»Nv) 

-  (.Aim  ^Aiki)  Nix  4-2  Nti/  +  8i,x) 


-Wp'ito T?(A)  c;„  a 

-  JyM)  Cjlj  + 1  K"(S)  ©jl  +  /cg(645';r>”«' 

-  +2f  Ly  M) Q.x  4  - /.y (4 

+  ^7'W  l£«  +  Mfa+ftTg-  riqlB)  T£3 

4  A/yW )  Cj  +  A/iJ  M  cf-  <j  M )  Cj" 

4  O"l(S)Z>”i0;(BW)V’H-O”(B)7fJ+£=0  (A-76) 


(in')  [i  =  1,  2,...K;  weighting  function  is  sin  Eq  A.7l] 

K 


du  Bh.mt  -4-dn(^)  A,„  -  (idu -H-du)  (.d)  BZm 


+chi($f  Bi"f  ■6,,'pZZx  -  (2  £» -■&»)  (#;  c  ;7. 


,4  **/ 


1^1 

4%  6^1-4,  UStSD/C 
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IBr+BlHfjtzA J  )/»■  t  (B'.+BUtf'-HAiV \J 

+  (8T+sD&C+H?  &’ 


+  (e,’“+Bl)(B*+2  8* )  S “  +(g'+  B° He’* +2  Bj  )  s' 

+  <«:+&•;  ( fl?«  ^OsT-  2  (  B?+8°iKB?  +2Bb  Bp] 


-H 


“  Nxy) "  (S{,«  +^i,xx)A Ixx 

-2N,^)(A^+Ai)-j(V%  f  Qy’(a;  C^x 

+  Q^(8+s*;Q"«  - 

+  ^  wd;«  c^(SKp 

+ ^  (S  t  e;  c’x  -.S^q'J-X  t*"m ;  x£ 

4  t$<A  H'lVp-  Jp)T>J.,)  +  UpB)C*  -4  U,VX’ 


-Uy(B)Cj  +  Uj  (/? J Dj  +  l/ij  (i4+/1)’Dy  -V^-CO^J 

+  £=o  ^-77j 

Finally,  the  Boundary  Conditions  [ss-i,  CC-i,  Eqs  A-37,  and  A-40  -  A-46] 
are  also  expressed  in  terms  of  the  dependent  variables,  through  the  use  of 
Eqs  A -47  They  are: 
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A  0,XX  (dll  -jq)  -  (-  OizNxx  -tOnh\ly)+Mxx  +3n'RxJ'--£ilNx/ 

Ai--  Bi  =  0  .  . 

di\Ai,**+&*tCi,x)'-hdi)(%)Bilx=  o  -  k 

dnBi.xx  i  X}-,  tf  -2  d/3  ( %)Alj  - 0  - 

Ci-Vi.x  ~I>i  =  Ci.y  ~  0  ;  i-  !>*> . 2K  (/}- 

2  A=  0 

^ ) =  (-ttnNx*  4  A 7*)  -ftf*  +£„NXX'8}IN*> 

Ai-B^O 

4 i,*>  +  ft,  G,*X  -  &  (#/G  4.2  tf/3  (# )  &,*  =  0  l  4 ■ ■  />*? .  ■•  ■ ■  * 

dii  Bi.xx iQx{Di,x*  ~&nC$) 32±  -  2dul^)/\i,x  ~  0 

72>i.x-  G.,x  '  0  )  i  -1,2,  —  2k 

<3iiCi.xa  -2 Qu(%)ViM  +Q,i(i£fl>i+J}t,Ai.„A(lGa--&)\ 

+p  4,  -  ■&»)  (%?Ai + 4s -J?  £  f  amt  aL 
-K'-fA+fcKi-ji'A'y]  Aj.*  4-f.a+ifB^ 

+  "-lytKHiifBlj,] BjJ  =  0  ;  2k. 

CU.  Vi,,**  +2 (Xu (4?)  Ci.nx  ~  (Aiilq?)  +  Ai  5i.««x 
■(^  &ts  -  3, i)($)A i.xx  Qu - Q*.)  (1$)  Si,, 

&-&IK-  u+iUZi 


c 


+<•  i-fc+fo  B  i,a  +  (WifBq 

+  j=  0  ii-1-3.--.2K 


SS~3  4.  =  0 


© 

/I 


Vl*+ &i(i&)Ci,x 

da  Bux -f  O-u'Dijx  -iSjii^Khx  =  0 


••• 


CL--  = 


^JiG,xy~^ii(^JZ^J(+^l,/4£/XX  •+2&i(-%)8L,x-  o 


&xxT>i,xyi Gv(i$)Ci.)ii -$iiBi,xx  -2&jA;.x  -  0  J 

SSr4  4.  -  o 

/Ivx  (  $,,Nxt 

Ai=  Bi=o 

6., Cl. XX  -  4„  (%fCi  -  ■&,  (#U'.X  +  dJl» 

(hJX.xx  ~  8"  ®  +8iid^)Vi,x^du  Bi,tx~2dis(.T?) 


* 


79 


(4-7?; 


“  ^i/A/yy 


•-2/C 

(4  -^o; 


’S’sihlty 


-a,AWa  -aa  <%)  a*  +&.AU, 


-0,ityfvi+2Ci»L%)viM*(aBia»)Ci,,L%1-cutCi...* 

^1%  jtj  rfuj +  O-tiitmi)  An ji]  Aj.n 

+  ft  i+i)  ‘  8ltj  -K- 1  -  7, .  t  %  8,1.,  ]  Bj,x}  =  o 

aatyjc.  - 2 au(#)Ci, ,,+  (#»*  a,oi x,($f- 

—  -  (-2 ' $n)(-$)  Bi.x  ~  -$*i  ~  *^r~ 

^ra  ^ Aiij  ■l(l~%-i  +  t?i)  U-jf  Ahj i)  Bj,x 

’+C(Hj)1Bi+j  4  +  i-j /  Buy ]  /^'./) 


}i--l.2.-  2K 


)(A-Vi) 


/\o  ~  A o,x  =  o 

/Al  =  Ai,X  ~  &l  ~  Bi.x  =  0  ;  i-  !,  2.  K 


Oi  <X  —  G.x  =  O 


-OalWVi  13  Q»(g)Dt,„  -QUC,„,  -  e»Ai. • 


(  Sil  ~2  &3i)  (Q)  Si, XX  -  0 

(^)Ci  -2GiiL%)Ci,M  ~Qxi 'Di.xxx  ~  Q-ji  Bi, 

~  (■  @31  ~2  ■&*)  (  Tf)Ai,x  x  -  0 


i-t.2.--- ,2K 


(,A  S3) 


CC-3 


Ao  -  A*,x  —  o 


Ax.  ~  Ai>k  -  Bi  -  Bn.*  —  0  i  i  *■/, <2 .  ~  •  >K 


Ci =  ~^i  -  ChiCi,**  ~ Gh +  GiiAi.x*  ~o 

Gi».T>i,xx  +Gu(%)Ci,x  +  -QnSi.xx  -  o 


i  =  l,2.  -  2k 


(A  -84) 


Ao  ~  ».x -  o 

A  i z  A  i,*  =  Bi  =  Bi.\ 


;  i  =  K 


~Gu  (iQ)  Ci  f  GxxCljtx-  Gu(lf)Vi.X  +  Ai'*x-0 

-a, 3  (9/Pi  +-?a.i(#)T3wlr  +  fan +ali)(^)tci,x 
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L-I.2.--.2K 


~  dixCi.t**  -  -&1,  Ai.tvx  4  (  -&H  -2@}i  )Ci&)Bi.tyi  ~  0 
-  a^tDi  t  a>.-Di.„  +  a»t%)G,,+4»Bu'  -  ° 
cu(#>*c.  -laufyCi.,-,  ■>  ( a„  +a,.H%)‘ v;.* 

~  $ii  thirty  “  &it  Si.xxM  ~  (  &1I  ~2$u)(tz  )Ai.*)(  -  ^ 


(ASS 


A .  2.8  Solution  Methodology-Finite  Difference  Equations 

Before  casting  the  field  equations  into  finite  difference  form,  the 
linearized  ordinary  differential  equations  of  compatibility  and  equilibrium, 
Eqs  (73)  -  (77),  can  be  written  in  matrix  form. 

CM  jfX.TO.j  +  (M.J  {/,„.)  +  [M)][x.«.j 

+  [Mjfx.xj  +  [Mi]  {Xj  +  {Mj  =0  (A.  86 


where 


.*+1 


.ml  wtl 


MX""'  A*"'  rJ"*l  ***  _ **i  _ **' 

'  L  A.  ,  >Ak  ,  B,  ,  '  8*  ,C(  *  .C„  ,  Dt  - 


(A  -87 


is  the  column  matrix  of  the  unknown  function  of  position  x,  and  [M^],  j  ■  1, 
2.... 5  are  square  matrices  [(6k  +  1)  by  (6k  +  1);  see  Eqs  A-73-A-77]  with 
elements  composed  of  known  parameters  (applied  loads,  geometry,  and  values 
of  the  unknowns  evaluated  at  the  previous  step,  m  and  therefore  known). 

(m_1  is  a  column  matrix  of  known  elements. 

6 

Next,  transformation  equations  are  introduced  in  order  to  reduce  the 
order  of  the  linearized  differential  equations.  This  step  increases  (doubles) 
the  number  of  equations,  but  it  is  introduced  for  convenience,  because  it  is 
easier  to  deal  with  low  order  equations  when  employing  the  finite  difference 
scheme.  These  transformation  equations  are 


{tj  =  [X.xx} 

and  they  are  used  in  only  in  connection  with  the  third  and  fourth  derivatives. 


By  this  transformation,  Eq.  a-87  Ed*  A-86  becomes 

rlxJi  N)  f{*h  . 

(R]  N  +CsW  +  =  ^ 


where 


(A  -SS) 


03-- 


[T]  - 


[0] 

CM.]' 

[S]  = 

(I] 

Co] 

0*] 

H  {sj  • 

Co) 

-CD  ; 

4 

(H)  (M.) 

(o)  Co] 

-CM*)) 


(/l  -91) 


The  governing  equations  (linearised  ordinary  differential  equations)  shown 
in  matrix  form,  Eqs  A-8&  ere  next  cast  into  finite  difference  form.  The 
usual  central  difference  formula  is  employed  and  the  equation  become 

ft  1 


(Ktf^csn 


,  M 


+  Clf  ) 


ffjfi  V''1’ 

j  <<lj  (4-90) 


where  j  denotes  the  j  th  node  of  the  finite  difference  grid.  At  each  end 
(x  =  0  and  L)  one  more  fictitious  point  is  used.  This  requires  (12k  +2) 
additional  equations  at  each  end  [the  total  number  is  (24k  +4)].  These 
needed  additional  equations  are  the  boundary  conditions  at  each  end,  Eqs 
A-78- A-79,  (whichever  set  applies  from  SS-i  or  CC-i)  and  their  number  is 
(12k  +  2).  The  boundary  conditions  may  also  be,  first,  expressed  in  matrix 
form  and  then  cast  into  finite  difference  form. 
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at  either  x  =  0  or  L 


tNj{x,„)  +  CN,](X.»  j  +CNJ  +fA/r]  =  0  < A-fO 

where  I  N  ~j,  j  =  1,  2,  3,  4,  are  matrices  [(12k  +  2''  by  (6k  +  1)]  with  known 

L.  J 

element,  and  'j  in  a  column  matrix  [(12k  +  2)  by  one]  with,  also,  known  elements, 


Use  of  the  transformation  equations,  Eq  A-87,  yields 


N  jg 


i  ’  =  (^G-] 


where 


(8S]  =  [CNJIIM,]] 
[BTJ  *  [HKN,]] 


(>1-92) 


[8C,]  =-  |NS] 


(.A-?}) 


Note  that  [bs]  and  [bt]  are  square  matrices  [(12k  +  2~)  by  (12k  +  2)].  In 

A  V 

finite  difference  form,  Eq.  A-92  ,  becomes 


^rss/jgf' 


rsTi* {^H*-  -t  rssi*  fas,}1 

m  r  j  lf?il  ■*  3  1ft)]  i 


<,A-W 


where  j  in  the  node  number  at  x  =  0  and  x  =  L(1  or  N) 


A .  2.9  End  Shortening,  Average  Shear  Strain  and  Total  Potential 

Before  outlining  in  detail  the  numerical  scheme  of  the  solution  methodo¬ 
logy,  it  is  necessary  to  write  the  expressions  for  the  average  end  shortening, 
average  shear  strain  and  the  total  potential  in  terms  of  the  dependent  varia¬ 
bles,  A^,  Bj,,  C  and  D^. 

The  average  end  shortening  and  shear  strain  are  defined  by 
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© 


© 


* 


<• 


<• 


<? 


T-§£<^ 

-'O 


'o  ✓  o 


^  =  2TLRL  £  /#  (If  4  4%)dxd*f 


In  terms  of  the  variables  w(x,  y)  and  F(x,  y) ,  the  above  expressions 
become : 


_  , _  i  /W/t  _ 

=  &n  Nxx  ~  Q-H  h Ixy  ~  JJlfiL  J0  J„  £  'V/  4  ^'x*  ~  ^/J  ^'*Y 

-h$„  vy/,w  +tSj  w>yy  +  2-6*1  K*y  V*  +2<)]  dxdy  {A  -96 ) 

*  -A*N*+foNv+imf.  i  t#»F'yyl&iF>*x~fi*F>v 

*9*yy  +2-Q& v^.»y  -  j’h (x(W.y4-2VV.y) 

t  j^(v^+iw,x'>JflUc(y  (4-97; 


Finally,  if  the  expressions  for  w  and  F  are  substituted  into  Eqs.  a-96  and 
A-97  these  equations  become: 


Qw  -  On  Nw  -  B,}  Nx/  "a  /.  i  4.  ~  A  */R  +  OnN-u 

(s$£i't(Ai*2Aj)Ai  +  CS;  +-?8j)S^J 

+  4."-  2  M.'+2 4.’)  4 s; (a;  tiB-yjdx  , A - ) 

ftv  a*Nt,  4  fli5  n«,  +Tf‘{-§fJ-  A.  4."  -  4.^ 

4  (?,j  Nxx  -  C?o  Nx>  +  (a?) (4j  4^’)  4j. 


+  +2  fip  8j]j  +  •&.  K  -  £  (iJ/C  ^  (Bj+Bj) 

~  Bj(A^Ai>  +*!'%  -  Wl dX  U-K) 

Similarly,  the  expression  for  the  total  potential  is: 

~  2~f0  £(Nxx£*x  ^  NyyQyy't  Nxy'fty*  ^X/U  ~Myykyy 
-2Mok„)dXdy  gwdxdy  -ft -RuU 
+  Nv  If]  l  dy  4  f  ( pja  vj,x ) /‘ dy  { A  -too) 

where  M  =  -  EN  and  E  is  the  load  eccentricity  meat.  ~ed  positive  in 

XX  XX 

the  positive  z-direction  and 

4=/;^  *  v£*iigdx 

Thus,  the  contribution  of  the  in-plane  loads  to  the  total  potential  becomes 

-nr-  NxxU  +  n„  2/J  /*  dy  ’  -  [M(-  Nn fgf*  +  fyl'&dtjdy 

In  terms  of  w  and  F  the  expression  for  becomes 

JT  Jg  J  [fill  f~>* y  ^  F>xx  ^^33  F.xy  +  -2  QnF,xxf’,yy 

~  3  fill  F ,yy  F,xy  - 2  fij3  F,m  F,*y  ]  dxdy  -  yf9  J*  (ch,  iv,* 

4  da  ^Vyy  4  4  da  ^4xy  4  ■?  ^/j  ^/<xx  ^4  yy  "4  4  d/$  lV,xx  W,  xy 

4  4^aj  W,yy  IV, xy )dxdy  -  NuJt  [  (fill  Fyy  4  F,  XX 

”  fiii^,xy)dxdy  4  A/xy/#  (fiaF,xx  ~  fiuF, xy 
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+  a II  F.yy)dxdy  l  Zvdxdy  4  7Efit 

•l&iiNxy  7LRi  (&4\,  Nxx  Y#j  Nxy  )  ~  2  TlRLQ/i  Nxx  A/*y 

✓  —  —  it 

~  j6  (  H  Nxy^'x)  jo  dy 


{A -so  i) 


Finally,  the  expression  for  the  total  potential  in  terms  of  A^,  and 


D.  becomes 
1 


Ur  =  7iRj;(±{- &»Al-A'M  + 

4(8 j+2Bi)Bj]  40„N»-  G.nH^+2{a*N»~Ql>Nn) 


"Zfcl'&iA*  ~A'/r  f  [(Ai+2Al)Ai 

4  (Bii2B’)Bi]  4  Q,,Na  -  CL}N„ J  -  d„  (Ao)‘ 


4  iz.{a„  (ffiCi  4  d'  j  4  a,  fccr/4  c  d;/J 

4“  l 


4a»{#/f(Ci'/4<BkYJ  (j?  JCCfCj  4  #Z>£) 

-2a*(%i(-  gd;  4  D£ci')-^<3ac-g)cci*23i-t>;,ci)J 

-i£fd„  f  (4.74  (8i)‘J  i  clU%f(Al  iBi) 

4  4  ^33  (-^/f  (/4  i  )l  +  (  5/)*J  -idnij})  (AlA i  4  5* ) 
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LG 


~4dn(Tl)(-Ai8i+BiAi)~  4du($)  UiBi~Si  4.  )jj 

-  Kl&'Wh  +  -WRU-enN,, 

4  ]fcvb]ty)  +  TIPI  (a..  '^GuNtyNxy  4  G.nNxy  ) 

-  47[gN„/?4.' 


(4-/02) 


Before  leaving  this  section,  it  is  important  to  give  the  expression  for 
the  modified  potential  an  expression  needed  in  the  estimation  of  dynamic  cri¬ 
tical  loads.  As  explained  in  Ref.  39  the  modification  is  associated  with  the 
deflectional  response  of  the  system.  When  an  axial  load  is  applied,  an  axial 
motion  will  result  (with  some  related  transverse  motion).  If  an  instability 
of  the  type  described  in  Refs.  40-43  and  37  is  to  take  plane,  under  sudden 
application  of  the  axial  load,  it  should  not  be  expected  to  occur  through  the 
primary  axial  node,  but  through  the  existence  of  transverse  deflectional  nodes, 
unrelated  to  the  axial  node.  Because  of  this  and  since  the  governing  equation 
for  dynamic  buckling  is  (though  conservation  of  energy) 


U 


X  ,  +  T  =  0 

Mod. 


where  T  is  the  kinetic  energy  (unrelated  to  transverse  deflectional  modes) , 
then  the  modified  potential  must  not  contain  in  plane  node  terms,  when  sud¬ 
denly  applied  in-plane  loads,  N  and  N  ,  are  considered.  In  the  case  of 

xx  xy 

lateral  pressure,  the  modification  is  different,  therefore  the  expression, 
given  below  for  the  modified  total  potential,  applies  only  to  in-plane  loads. 
This  expression  is  obtained  by  excluding  strictly  load-dependent  terms  and 
those  terms  related  to  F(x,  y) ,  [cD//],  which  correspond  to  in-plane  motion. 
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ft 
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ft 


ft 


UT„*i  =  UT  +  krl  Cn«  (a..  -  at /a,,)  +  N„  ( 

-  ai,/a„) 4 2N«N,y(a«  M -'®« 

A.  2.9  Solution  Methodology  -  Numerical  Scheme 

A  computer  program  has  been  written  (see  Appendix  A  for  flow  charts 
and  Program  Listing)  for  data  generation.  The  linearized  finite  difference 
equations  are  solved  by  an  algorithm  which  is  a  modification  of  the  one  de¬ 
scribed  in  Ref.  43.  The  modification,  which  consists  of  a  generalization 
of  the  algorithm  of  Ref.  43  is  fully  described  in  Appendix  B.  The  solution 
procedure  used  for  the  problem,  herein,  is  based  on  the  algorithm  described 
in  Appendix  B. 

The  field  equations,  Eq.  A-90,  can  be  written  as 

tcjflj  *  (6.J  IA-hq 

where  K  =  1,  2 . N  and 

(5j*  -kin)’'-  ;  raj*-aiprtt*+CT3* 

CA,]- Tel*?* &  CSJ*  ;  }f*}  =  (Wj*  C A-,oi) 

Note  that  there  are  (12k  +  2)  elements  in  the  vector. 

In  addition,  the  boundary  conditions,  Eqs.  a-94  can  be  written  in  a 
similar  [to  Eqs  A. 1053  form, 
at  x  -  0  (k  -  1) 


'  [ c.3 f  *.] +  C  s, ) a,}  +  (A ) [  1}  -  /86,j 


and  at  x  *  L(K  -  N) 


-CC„]{§4  +  +  flJf*"'] =  f8G"j 


(4-/07; 


(4-/OJ) 


whefe 


CM  -  [bt]L;  CAJ'izC Bsf 


(A- /or. 


i-  I." 


Note  that  j^Z^J  and  \ZN  +  jJ  denote  the  vectors  of  the  unknowns  at  the  fic¬ 
titious  points  (k  =  0  and  k  =  N  +  1). 

By  properly  arranging  Eqs.  A-105,  A-107  and  a-108  for  the  entire  cylinder, 
the  following  matrix  representation  is  obtained. 


m 


Ci  IBilAi 


Cw  Bj*i|  A 


Ch  Bm  Ah 

Cv  Bh  Ah 


Eq.  A. 110  can  be  put  In  the  torn  of  Fig  C.l  (Appendix  C)  and  it  will  be  a 
special  case  of  this  form,  by  the  following  changes.  First,  there  is  no 
common  unknown  vector  and  thus  all  the  {d^}  vectors  are  zero  (tridiagonal 
matrix).  Next, 


(-24**4)  by(HK+4) 


(-24*44)  by  07H 


(-34*44)  by(/JK+l) 


CCJ  «  CCOJ  CC.3] 


eg  -  c  <g 

CSi)  «  ts*) 
C  Ail  * 


(-24*  +4)  by  071a. 

(tt*42)  by  (24*44) 

i  =  3,4.  —  .  Af-  I 

j  =  .2.3. - .*  Af-2 


I 
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M-  '[!»] 


l  Ah.,)  =  CC^J.fo]] 


CCJ  = 

r  bn)  = 


rc„) 

(6] 

-  4 

®j  a„} 

,(§J  til 


(U*+2)  by  (24*  + 4) 


V4K+4)  by  (M +1) 


{24KW  by  (24K+4) 


(24K+4)  by  Om 


(.2AMH)  by  0/U 


Note  that  m^  «  *  24k  +  4,  while  m^  *  12k  +  2  for  i  =  2,  3,  4,  ...» 

N  -  1. 

Note  also  that  Eqa.  A-110  represents  equilibrium  and  compatibility  equa¬ 
tions  in  which  displacement  components  (A^,  B»)  and  stress  resultant  components 
(Cj,  D^)  [see  Eq.  A-86J  are  the 
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unknown  functions,  while  the  geometry  and  the  loading  (taken  in  increments) 
are  taken  on  known  parameters  (assigned  everytime  the  equations  are  solved). 

Thus,  this  special  case  of  the  algorithm,  Eqs  A-11Q,  is  employed  for  find¬ 
ing  pre-limit  point  response.  When  approaching  the  critical  load,  the  in¬ 
crement  in  the  applied  load  parameter  is  kept  small  and  the  sign  of  the  deter¬ 
minant  of  the  coefficients  [d  in  Eq.(C  -  19)3  must  be  checked.  If  convergence 
fails,  the  load  level  is  over  the  limit  point.  But  if  convergence  does  not 
fail  and  the  sign  of  the  determinant  changes  from  what  it  was  at  the  previous 
load  level,  then  the  load  level  is  also  over  the  limit  point.  Desired  accuracy 
can  be  achieved  by  taking  smaller  and  smaller  Increments  in  the  load  parameter. 

It  is  also  observed  that  by  employing  this  procedure  (special  case  of  the  al¬ 
gorithm  in  which  the  load  parameter  is  known),  no  solution  can  be  obtained 
past  the  limit  point.  Because  of  this,  the  more  general  algorithm,  described 
in  Appendix  B,  is  employed  at  this  point  of  the  solution  procedure.  The  new 
and  more  general  algorithm  simply  changes  the  role  of  one  of  the  displacement 
terms  with  that  of  the  applied  load  parameter.  By  so  doing  the  form  of  the 
equations  changes  and  the  matrix  of  the  coefficients  of  the  unknown  ceases  to 
be  tridiagonal.  Depending  on  the  position  of  the  particular  term  that  replaces 

the  load  parameter  [which  one  of  the  (6k  +  2)  terms,  and  at  which  node  (x-posltion)3 

f  ' 

column  matrices  appear  all  along  the  column  corresponding  to  the  vectoT 
and  the  new  equations  assume  exactly  the  form  shown  on  Fig.  C-l.  Thus,  at 
some  level  before,  the  limit  point,  the  procedure  is  switched  to  the  more 
general  algorithm  (Appendix  c)>  in  which  one  of  the  displacement  parameters 
(Aj  or  Bj)  at  some  specified  node  is  taken  as  known  (specified  increments) 
and  the  load  parameter  is  the  unknown.  This  solution  procedure  is  continued 
until  the  desired  portion  of  the  post-limit  point  response  is  obtained. 


Finally,  in  generating  data,  numerical  Integration  is  used  to  £ind  the 
values  of  the  total  potential,  the  average  end  shortening  and  the  average 
shear  [see  Eqs  A-102,  A-98,  and  A-99]. 

A. 3.0  The  u«  v.  w  -  Formulation 

The  geometry  and  sign  convention  for  this  formulation  are  shown  on  Figs 
A. 3  and  A. 4.  Note  that  for  this  case  the  x-axis  (and  therefore  the  transverse 
displacement  component  w)  is  taken  as  positive  outward. 

In  this  formulation  two  distinctly  different  kinematic  relations  (dif¬ 
ferent  shell  theories)  are  employed.  One  is  due  to  Sanders  (Ref  34)  and  one 
due  to  Donnell  (Ref  33).  In  the  case  of  Sanders'  equations,  it  is  assumed 
that  the  reference  surface  strains  are  small,  the  rotation  about  the  normal 
is  negligibly  small  and  the  rotations  about  ln-plane  axes  are  moderate. 

One  of  the  reasons  for  expressing  the  governing  equations  in  terms  of 
u,  v,  and  w,  is  that  it  is  not  possible  to  define  a  stress  resultant  function, 
in  order  to  satisfy  the  in-plane  equilibrium  equation  identically,  when  using 
the  Sanders'  kinematic  relations.  The  case  of  using  Donnell-type  kinematic 
relations  is  a  special  case  of  the  Sanders  case. 

A.  3.1  Kinematic  Relations 

The  kinematic  relations  derived  by  Sanders  assume  a  perfect  reference 
surface.  These  kinematic  relations  (Ref  34)  are  modified  to  Include  the  ef¬ 
fect  of  an  Initial  geometric  imperfection  w°(x,y)  as  shown  below. 


£xx  =  £xx  +  k*x 


where 


Gxk  =  +  2  lV,x  ■+ vv/,*  vv/,* 

--  */„  +  #  4  jw,;  +  w„k/,;  4(j!  ("**•%» 

^*y  -  ^  +wij^  +  w.k  v.ji  +  iaA*  w,y  -  S  i^(*Ax  +  j 

<P,  •  -  w,x  J  ^x-W.y+S.^ 

X„  =  -W,xx  ;  *„  =  _  +  J,  Ut 

x,y  *  -  w.«y  +  i s;& 


{A-"S) 


where 


Sr 


r  i 


•/•A  SmJ*K'  k;nt+njtic  ajlLL 


(A-"*) 


0  "Dumtll's  kifumtu  aJL&** 

A.  3.2  Stress-Strain  Relations 

The  constitutive  equations  are  the  same  as  in  the  w,  F- formulation.  Be¬ 
cause  of  the  different  sign  convention  the  relations  between  the  stress  and 
moment  resultants  on  one  hand  and  the  reference  surface  strains  and  changes 


in  curvature  and  torsions  on  the  other,  these  equations  are 
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where  the  expressions  for  A^,  end  are  given  by  Eqs  A-14  and  A-15. 
A.  3.3  Equilibrium  Equations 

Following  the  same  procudure  as  the  one  described  in  section  A. 2.3, 
the  equilibrium  equations  and  associated  boundary  conditions  are: 
laui librium  Eouations 

Nx*,x  A/» y,y  =  ^ 

N«.t  +  N„,  - 

+s,-^  +  $,%-’  =  o 

C  Mn(w«<+tv.',)],,  +  CN,r  (w,,+ w,‘ )]iX  +  [  Nyy(w,x+w.J)],, 

+  ^  -  -^C(Nwu),»  +  (Nnu),,) 

+  M*,.ior  +  -?Mxy,xy  4  Myy.yy  *  $  =  0 
Boundary  Conditions  (at  x  -  0.  L) 

Either  Or 

Nxy  =  Nxx  £  ^  ^ 

+  Su*  0 

N*tWx  +  »/'x)4-  NqCW.yM'y) 


-  S.lfZ/  4  M*r,X  =  Q*  +^y.y  Ew  -  0 


Use  of  the  first  equilibrium  equation  in  the  third  yields 


( a/,y  f  vs/.y )  (  A/«y,x  4  ^  ^  W,**)  4  2Nxy(W»«y  +  W.xy) 

+  Nw(w,w+w/,i)  -  ££ 

i  tf  JU.t*  2Mxy,xy  f  Myy,yy  f  =  0 

A. 3,4  Solution  Methodology-Field  Equations 

The  solution  procedure  for  this  formulation  is  as  follows:  assume  a 
separated  solution  for  u,  v,  and  w;  express  the  known  (assigned)  parameters 
w°(imperfection)  and  q(pressure)  in  a  similar  form;  find  expressions  for 
reference  surface  strains,  changes  in  curvature  and  torsion  and  stress  and 
moment  resultants;  substitute  these  expressions  into  the  equilibrium  equations 
and  use  the  G  alerkin  procedure  in  the  circumferential  direction  (this  changes 
the  nonlinear  partial  differential  equations  to  a  set  of  nonlinear  ordinary 
differential  equations);  use  Newton's  method,  applicable  to  differential  equa¬ 
tions,  to  reduce  the  nonlinear  field  equations  to  a  sequence  of  linear  systems; 
finally  cast  equations  into  finite  difference  form. 

All  of  these  steps  are  shown  herein,  in  detail.  Then,  once  this  step  is 
completed,  the  solution  scheme  of  Appendix  B  is  used  to  solve  the  final  set 
of  equations. 

The  dependent  variables  are  the  three  displacement  components  u(x,y), 
v(x,y)  and  w(x,y).  a  separated  series  form  is  assumed  for  each  of  them 

UOc,y;=  Z  (  U,iix>  +  Wit*) 

i  *0 

inx.y)  =  £  C  ir,i(x)  Cca^;  4  Unix) 

1*0 

W(*.y)  =  +  WxiWAltll?]  iA-119) 
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Thus,  the  number  of  unknown  functions  of  x  is  (2k  +  2)  for  each  variable. 
The  total  number  is  (6k  +  6)  subject  to  the  condition  that 


W,  o=  U.  =  Wao  =  0  (A-'*0) 

Note  that  the  true  number  of  unknown  functions  is  (6k  +  3) . 

Similarly  the  expressions  for  w°  and  the  pressure  q(x,y)  are: 

vs/'xj) r  JL  C  W ^  +  Wififanig)  (A-tJi  ) 

1=0  ^ 

£(  2*  (4-/22) 

°  i*®  O  V 

In  this  case  also,  the  condition  W£q  =  =  ®  *"s  ^P08^* 

In  order  to  express  the  equilibrium  equations  in  terms  of  the  parameters 

of  Eqs  a-119  "A-122»  one  needs  to  first  find  the  expressions  for  the  stress 

resultants  and  therefore  reference  surface  strains  and  changes  in  curvature 

and  torsion. 

Use  of  Eqs  a-119  and  A-120  in  the  expression  for  andn^,  Eqs  a-112 
and  A-113  yields 


6xx  =  S 1  Uli-*  +  ifli  +  til  ) 

4  ( Si. Mu,x  +  ill  +  il i)  (A -Hi) 

where 

t  l  =  A  l  OO  (  /  Mi,x)  f  ^  •  VJi.n  ) 

t»i  =  AiUoWi*.  W,.*)  ■+  A3W  Mi'X.Wi'x) 


th  -Ha  kkjCMi.x,  W,,y)i  AijWi,*,  Mi,*)} 

~txii  ~  ~2  l  Auk)(M,,x,  bU,x)  i  Ai(p(  Wi,*/  W//Xj] 


£y y  r  f  t  ^  Wi  ■+  *^)  £  i  4  C04  Mjjf 

-tCC^Vlil^jSi  4  ^y2i  4  -^y2i] 

where  / 

V=  (W1f/ll‘j«(*)(W/-'V.)  + 
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1  L  J  * 

Tyai  *  "(tf-)  4  >4  ZJ3(/g  (W.'/tV*)] 

4"  ~jp  04  J3(k)  ( <^» >  tT, )  -  ^  jjk;  (W»*,l/i)J 
^y/i  =  J  ^  i?  )  64  ij4«u  (  Wn  W4 )  4  >)£,  tk)  ( iva  yvVoj 

+  +Ai<*>(Vi,UO) 

4 (AwwW'.Vt)  -Asm  cw«,tr#)) 

j./n\x,  a i 

^ni~  '2U/  + 

+dw*  (AtooiVi,v,)i  Axo^cUnV,)) 

^  $  (A  33(X)  (-Ml,  IT, )  "  A1100  (Wt,  1/i)  ) 

K  ■  + mu. t  iL,rn  t 

4  C  C“  i/,£  *f  Uai,*  )&.  4  f  "fctyi'il 
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4**»  = 
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(A  -/<?£) 
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-  AUvo^*'*  sv.)  i 

”  -7^  (  4i<jy  (Vi ,  W^j, )  4  A4W  (V**  Mi,*)) 

~t*yii  =  "A?  r^32(w  '  Mi,*)  -  Ajiot)  (Mi  ,  Wt,x) 

^  • 

“  (^>#x  /  Mi)  *f  4i  j  (K/!*,* , 

"  7?  [^3  Cfcj  (  /  M,,* )  f  y4i,*»  (Vi ,  Wi°x  )] 

/  7| 

~ty yi'i  z  2t^(.^JUic)(Mi,Mi,x)  ~  Am(fc)(M< 

"  Al4W  (Mi,*  ,W,)-b  AlUk)(Mi,x, Mi)] 
-$(A  1  Ik)  ( IT,  ,  M,x)  f  Aide)  M  /  M,x)) 
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(4-/2^) 

(4-/20 

(4  -/27; 

(4 -/.?$; 


102 


Note  that  6^  and  T]  ^  are  the  same  as  before,  or 


L  >K 


1*0 

l  >0 


iA-uv 


The  symbols  A(j)(k)  (i  =  1,2, 3, 4),  A(^(k)  (j  =  1,2, 3,4),  AIj(k)  (j  - 
1,2, 3, 4),  Ajj  (j  =  1,2, 3, 4)  and  Aj2j^  (j  =  1, 2,3,4)  result  from  the  use 

J  00 

of  trigonometric  Identifies,  which  are  employed  to  change  double  to  single 
sums  [similar  to  Eqs  A -49  **  A-51and  symbols  defined  by  Eqs  A-52  -A-54;  note 
that  some  are  common! .  The  needed  trigonometric  identities  and  definition 
of  symbols  are  given  below. 

jLr  (^iCrti^GLiCrtiO  =  15*  Ait*)d>,*>  C«Ud 

L-O  9  0  £s s 

|i:C °-iSinL9  =  J|  AlMa.a>  Smie 

K  L  *H  I 

x.j8}  QiCstif)  *  Ahk]  (b.Ck) Si  niff 

jiLlfySiniel  Qi  Si> tiff  * 4?  Atw(l.aiC*4i6  (A-Vo) 

jLtlifyCMjffJ  a,c.i  iff  = A (6. a) Co^iS 
g'feQqcsjriaiSixie  =  £„ 

V»vdaiG« **  (b.a) sinis 

4,4 Q-fySinjelfliSmiff  =  (4-/31. 


(4-/30 
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(A-'31) 


^  ^£#4^0}  (XiiGtfiB  ~  ^  Anwt'&'A)  &410 
£.£.{  'Sj  C#<jd') CUL  ^niB  s  |?  Ajufc)  (4.&)  ^ 16 
aiiC**iO  sf&A*3i*)(4-a)  'tonW 

t£  C Sj /Un^Clil  ^iniB  r  feAz4*j(£’a  >  °*ll8 

lzb4~0  ^  • 

'  Si  Ax jiwC^a)  cnie 
&£Mfy  C*si$]ad*ni0  =C  A^wC^aJ^iO 

&J.  '&*ty6]  CLilCwiO  >4ij3flo  (•&•&)  AiniB 

=5  An 409  <*'*)<*& 
^feU^oHjejCh  ceaie  s%  Anl(K,{8A)GxLe 

U.O%C^B)  CkAiniB  K  |?  AjnwL4>&)  <4*nio 

0.iC*$i6  =1^  ASuw  (■&.&) dM-LO 

£a£Q2-8js4Mj6]  Q-idtolLB  *£j>Aji40')(4>a) 
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Atm  i*  a)  =  i£c&v+('-&ffc)  6,,.;,)  a< 
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Q*"  V 

=  j  +f-/  +  ^+ft)  <f 

4»»^.«)  =  i  +  (- 1  ■%  *K)  4'*j>] i  ai  a'm) 

ALaiM(S.a)  =  f  £(  (W>4trHi-fc+fi)ii-4i4ltyJjai 

A  Loo  CM  =  ft*  +  fuyJjQ- 

aLo»cM  =  f|0f<t  ^  A*+H+^+fc;/w/M<fe 

AL»Mm  *  j  +<'-&  *%)(.i-i)&'H>)&i 

ALooCM^^c^^i 

Ajm iM  -  +fa a-ii4*#}(b  o ww) 

In  order  to  write  the  strain-displacement  relations  in  matrix  form  the  fol¬ 
lowing  definitions  of  column  matrices  (vectors)  are  needed. 
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Note  that  t^  and  t~  elements  are  given  by  Eqs  A-123a.  A-124a  and  a-125*  while 
the  and  elements  are: 


G.XXX  -  5i  l(xi,X 

^U*  =  ^4?  +  7?)  Si 

€m  = 

)fai  = 


l  kn 


;  Kwi  *  ** 


;  Hiyl  -  C-  ^  Wu  +3f  UuJ  Si 

;  (A-UV 


106 


I  L 

Note  that  the  an<*  t*2iJ  vectors  result  from  linear  portion  of  the  kine- 

matlc  relations;  the  |t^jj  and  ^t~^y  from  the  coupling  between  the  Imperfection 
parameter,  w°,  and  the  displacement  components  v  and  w  (thus,  in  a  sense,  non¬ 
linear  relations);  and  the  end  {t^jj  vectors  from  the  nonlinear  terms 

of  the  kinematic  relations  (v  and  w  coupling). 

Substitution  of  all  the  derived  expressions  into  the  equilibrium  equations, 
Eqs  A-116,  yields  in-plane  equilibrium 


4-  (Ttmi*  -  t  4  Tlftu,*  ~  4  %aui.x  -  j$7tiyn)dtol'7£ 
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«  • 

*  A  1  00  (Mi,*  >7l*n)) 

4  Q'ftnu  +  TUfll.x.  4  4 


(A-IS3) 


=  -§^(  >4«<»  (W/,^)  -  /Ui<*>t^  *yy, » 

4-^r  C /Ijtw  (  ^<x  ,W*y)  4/1  iatf  {^i'X  /?Z* y*0 

-  +  #*y*i,x  4  -J?#?xyji.*  ~  %  S'  W,yn 

Sm  =  4  C  4  -k^V/Ln  +  %*!*»* 

+  /|i#  ->4j4wtW,,#yyiO 

-  |  Anu(.V,,n„,)- 

f  jj/CA  i«)(W'.»,Ayi)  +  ??xyi)) 

§!a  =  -  Tfftyyil  +  ?l*yiL,x  +  ^ftlx,n.x 

4  ~  Ajic qC.'M'i  flyyi)] 

-  ^  +  A<»(^». ^ry*)3 

•^(/'fic*'/-.*  .Ay.) 

5,1  =  /4/««,(W,*.^,) 

+  f  [  A<«>  ( W,;  X.  >  + 

SJ  =  M[-  Aixox,  ( W/,  <W,*.<JJ 

4  Al on  ( w/,» , )  +  Aw 

5a  *  ^rl  Ajn>v(*/>w!-  *tm)  -  Aj*oh(*,-*w< »*W) 

-  V-.  )  +4iw(VV.»^0 

+  jlCAtumC'’**"*  ’*1,)  4Aiu»lW"*V>) 


M-/S-4) 
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■+  ^ ,  W#i)  “  A jmno  W'* Krtifi 

-  %l  A) mi  * V''  *ii\)  4  Aiiito  ( Vi  •  fyy>  )J 
+  ^  u>w  C  ix  i  Tttyt* )  <4<i  t**>  ( vi*,*'Wtyi )J 

^i0l  =  A ji<j*j  ( **,+*/“ .  7?#,)  A  ji ay (Wi*-*?,  7}fy>)) 

jpt  ^4j««ju  £1^ '  ftft*)  Aituo  tu*' 

4  'fifCAitm)  C^.X  +  ^M  >M*yi)+  A  it»v (*'•'*+ W,‘X '%*/*}) 

X  7t  L  L 

**  f”  ^  JJUJD  /4ji£»iy  (lV*,  Tfyyi  )J 

-  ^5  (  A  *mj  ( Vi ,  7l»>)  i  Aii»)  ( Vt , 

■+  ■§  r  >4  Jt»*/  (W*,x .  #1 w)  +  All*)  Wt.x ,7Ufl)]  (A-lSS) 

transverse  equilibrium 

&(  t'c**  + 1'^^) + gW'L+ti)c*% 

*!&«<*% +2**  f)  =  0  (/l  -/«j 

where 

Y«  =  **««*. + -2*w»  (£)  -  $j»wi  -  ■%> 

7,i  =  #?«!*.»«  -2%tqii,tO$)  -  (tJ)  Xi*!'  (A-K7) 

fc«=  M«a.»  +2ttwi --Sfu 

4  j?(  A ji  (*j  ( W*#/  #*/,* )  -  A  j4( q  (  */,*,  flW)J 

*  * 

4&)(  Axjiuo  ( W*  / fyyO  t /4 j ( W,#,  ^yy/-)J 
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4  £  4  IW  (  >  ft/HI  )  ~t  A  too  (  M/XX  ,  7ll(X*)J 

4^f  f/4  ji«w(*Vi,\,;?Vi)  -  AwqW't  'Tlty*) J 

*(%K~  Am otjW,  fyvi)  "  A  jm <*j  ( w*  >  )J 
■+(  ^ )  f  A  j  200  ( *)  -  Ajj  m  Km* ) J 

4  C  7?)  ["  Aljiut)  (.  i?lyyi )  “  A 133  00  (  ^*,  ?lyyX) 

4  Amo  ( 4  #*„ ) 

(Allot)  M*,X  Jit'll)  -  A  731/tjC  Jtiyi )) 

*(i)\-  ^  Jl3  (W»  i^yy/jJ 

tylin  ~  y%jnn,*X  4-?  Mxyuji}}^)  -  (&£)  7rtfjyL  -  “^p* 

‘)j$C  AjiotfW*>%Krij)  ~  A  ( w/( ,  ytxji.n. )] 

+$)'(A,7W(«2.n  vy*  )  *f/4  IJ4  (*;(*/»<  ??yy»  )J 
4$f  f  "A <ic)  ( Vt,7?ty,,x) -Ain  ( &  '%w,x)] 
4^C-AiIU)(v„x„t)  +Aitwiif., ??„,)] 

4  ^  too  (  ^.XX ,  %f )  4  4#(W  ( tfi.xx ,  7?n> ) 


4  7?  ( AjIOoW**  ’^W)  ~  Aj 4Q,)  (W,X,ff»yi)] 

"  R  IA?W(V,.K  >%y,)+jhlOolU.fni»)] 

4  $  /  (-  A»,w  ( w, ,  J  “4wa  M  (**.*»,)) 

4  i^s?  ("  >4  J'  ’  (  - Xyyi )  4  A  J4  wlVi , fyt)] 


-  2  Wx,tljx  (.$)-(*$)  Wym  ~  ^P* 

*(1)(A  ja  <*j  (  v^*'  W*y2it )  ~  >4  JJ(«  c  v'/l  >X*ybx  )J 

C ~R  )  t ~  A  IV  (  >  ?lyyi )  ~  Avi  ty  t  ,  97yy/)J 

t  ~R  l"  A  UK)  (  Vt  ,?l*yitX)  -  A'3w(Vi  ,9?*yi,x)] 

■+  /f1  l  A  2 2 ik)  Wi  ,9lyyi)  "  Axjoo  C7/xt97yyi  J J 

*  » 

A  a  uo  (  W*xi)  *f  j  fly  ,  97*xi) 

+  ^  (^Jidc;  C^V/^rya)  *■ 

-f  ^  (" /4i(»o( ^ix,#xya)- /)ioy  (Ut,x;Mxrt)J 

+  (j?)  (~  ^JaauyO^'Tfya)  -  djvwiWi, 7lyy/)J 

1  -7=? r  £*“  f 

Y'H~  tAji  <W Wl  '  MxyixX  )  -  Aj/t(2k>  ( w',  Wxyi.x) 

2  A  3HXK)  (  ,9?x/i)  ~-2  >4j4fuy  C  W/jX  ,  )  J 

**  ($)  C  »  #/ya.)  4  Aijauk)^’  '^7»  ) 

**  4j2K3X)  Wt  ,7?yyi)-  4*nuyW*  *^V»)} 

^  /(a*)  (  ^XX  *  ^W/ )  ^  A  4im)(  */*-XX  ,  Ttxiu  ) 

Vait  =  ««*>  (^i  ,97  Xyj,x  )~  A  33U*i  (  W,*  -  #*yM  ) 

*  “2  ,  ?7gyi)'m  2Aj Jl**)(Whn  ,97xy,)J 

+  C/?J  (“  AlJid*)  (Wa*,  JtJJ(aarj(W*,^yyi) 


i 

J22 


(\A/»°,  Tlyyt)  ~  A 


i 

Ji  J  UXJ 


( V*.  ??£,) J 


4  A  2  {.ikj  (  f  >  ??*,)  4  AiC**J  (  V^2,X*  ,??XXi) 

tyi*  ~  7£  £  A  jmk)  (  MxiW*  ,  7?xyi,x)  -  ^  J4£2*;(VV,  +lv/,  %y*  <jr ) 

4  4jk**J  (  V\/a,  %iy/,y)  -  A  J4dk)(^ij  y[xfji,x  ) 

^  2  A  JUW  +  Vi,*  ,77yyl  )  ~2  A  J4 U#J ^hX^i.x ,  Tltyi ) 

4  2  A  ji(2aj  ( ^K)tj  7l*yt)  -2  Aj4m(W,,x  ,7?^)) 

% 

J  ^ 

4  fjj)  M  W C w*4  'v»°  >  Tlyyi )  4  AiJ4(2AJ  (  #yy, ) 

•f  A  J3M*)  (*ti  ,7?yy2)  +Al J4W9  ^ •  rfyy,)  “  A  7lU*)LW,,7lyy,) 

-  AjMQJoWz.Ttyyj)]  4  7?C"*4i  «u*(#  #  %yi.*)  "  ^iow  ^ 

**  4  /  (utj  (Uhx>7?xyi)  ~Aaok)  (Vi,*  , Ttxyi )  “  A i  <.iv  (V\  / ^xy/,x ) 


—  Aim*)  *•  Vx  >7lxy*,x)  -  A  iiucf  ^  -  Am*) LU*.k  >71*yi)^ 

4 ^C“4ll(»J  CVt,7ly yi)  iAu<Mj(^i>7f^,)-AjU2Xj(^^yyO 

4  4j4LVi)(^,,  Xm)-Aiu*)  W,,7tffi)  iAx4(2k)(Vz,7lyyi) 

£  * 

-  Tlyyi)  4  A  j<iUK)(.U‘'%yyi)) 


4 /(a/y  (^f,*x  4  yTlxxi)  *f  4*{Uy  (  M,*X  41V*,**  ,  T?***) 

.  i  L  l  A 

4  ni(jx)  (>V»,xx  4  >7lxm)  (A-l£jc) 
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tyiin  =  (7?)t  4  -  A  jsottjC^^  '^yi,*) 

4  2  4j2oa)(^j/X4  ^ i,x  t^lxyi)"^A  >7lx/i) 

4  A  nay  (W*/7Zjfy*)j)  “/4jift*;(Wi,^xyWJ  42/4  jam;  W**tJUy*) 

-  2  An  OK)  WhuTlxyi  )3+6  4  X3i  £1X)  (  )  -  A I  Jitot;  ( Wi  ftn,) 

\ 

~  -4  C^'j^yyj)  ”  4  Jil  (2/y  ( Wl>%yyi)}  ( s’] 

4  (]?)  C“  4jJiOW^1+  ^  >  ^yy»i  ~  4  uicuj  (Wt+W,0, 7l?n) 

~  4  J320itj  (W,  +  IV,  ,  #yyt)  "  4  Jaj(2*j  (W»+W/ ,?7yyi)J 

4  C~4au«  (UtfTlxyi'X+ftgyijx)  '■^j(vej(l/i,Pt^i,xi^xyi.x) 

”42  OK)  (Lh  >x  >  TLtyZ  4  Tftcyj  )  ”  4j<*j  (  /  %fl  4  ^?*y)  )J 

■4  7p  [/4xzwc)(  Wh  ffyyi  4  ??yy/ )  -  Ali  (J*;  (  Vi,9tyyi  4  /fyyj  )] 

4  4^«*j  ^  ^/»J0(4  IV|,x«f  >%oc.O  4  4  3<J*)  ^ >7lx*\) 

*  « 

4  4it*j  (W,,*  )  4  Am*)  (WisXXyTl#,)  (A-iSfd) 

According  to  Eqs  a-119  subject  to  the  constraint  of  Eq  a-120,  there 
exist  6k  +  3  unknown  functions  of  position  .  These  are  the  displacement  co¬ 
efficients  u^x),  v^x),  w^x)  for  i  =  0,l,2...k,  and  u2i(x),  v2i(x), 
w2i^x^  *or  *  =  1*2*3... k.  Note  that  if  one  can  solve  for  the  displacement 
components  the  response  of  the  system  1b  fully  characterized  (deformation 
approach) . 

Next,  the  Galerkin  procedure  is  employed  in  the  circumferential  direction. 
The  vanishing  of  the  Galerkin  integrals  leads  to  (6k  +  3)  nonlinear  algebraic 
equations  in  the  (6k  +  3)  unknowns.  These  equations  are: 
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si  +^u+§,-^C*C  =  o 

C  +f«i  **£■*&  *1L  -  -2,1 

4*  i'-o.  1.2 .  •  •-•,* 

^Wi,x”  ^R^xyti  4  ftx*n.yt  ~  ^7lmi  =  -  P?XxiLy 

s*  +  %L + Kl  4  &1»  "  Sir  =  0 

C  +  K  +  ?4  +  +  {54  -  -^. 

^  i*  /,<?,  •--■  ,  *  (A-/6o> 

Next,  the  generalized  Newton's  method  (Ref.  38)  is  used  to  reduce  the 
nonlinear  field  equations  to  a  sequence  of  linear  systems.  This  procedure 
is  similar  as  the  one  in  section  II. 2.  Because  the  final  set  of  equations, 

Gqs  A-160,  contains  n's,  g's  and  T|'s,  and  because  these  are  in  turn  functions 
of  other  parameters .then  Bq  A-72  will  be  applied  to  all  of  the  elements,  need¬ 
ed  in  deriving  the  iteration  equations.  In  so  doing,  only  the  nonlinear  terms 
need  be  considered.  Thus, 


{£*li)  55  A'iOO  (  *  Wf.g  )  4  A^dt)  (  h/»,<  ) 

("kizS ~  A2uo(Wi,t,*fht)  +  Am » W/.r ) 


M-^0 
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( i#i)  s  $  C  Am*)  ( W,*,  Wi  )  4  A  ijiuoirt? ,  *T. )J 

+  jr(4L(v*  vr)  i  AiwivT.iT)] 

4  ~f['  £  Aj4(V  (W?, ,14  )  +  A  X4Ut)  (  J  W§  j 

-di,  <*"0  -aU(»w>] 

"  2ft*[  AjJ4»)(K*>  K  )  +  AlJ tw  (W>  /  W,J] 

'  &  [Au«  M”  V?)  +  4i« 

**  “/p1  CAj4l/tj(M,  jV»)  -  A  Jl(W*,lS,  )J 

t&r- 

+AU>(v,*>vr')) 

^  ^§il  AsitytWi  >  &  )  +  AnW  (Vi  >VJ,  ) 

-  A  3109  (Wi ,  uj  ;  -  /4  *^4 

4  iAwCM,"*?)] 

-  I#*Crfsw(  V**,Z/f”)  i  A  i  uo  c  v,*,  wT. ;  J 

-  ^3  f  dj30c}(W,  jtf  )  -  AjJ(Kj(W**  U  )] 

(  iyii)  =  ^  f  /4  J»  uo  (  W? ,  whx  )  -  A  Jfflcj  (W,"  Wi.g  ) 
lAituoM*  ,wD  - Ainw^Z >*/!**)] 

-  %  lA\mCV"  w%)  +aLiv?M Z') 

+  aL  c<  .  vr;  *f 


~  C  Ajiw  (*V»  j  Witn)  ~Ajw  (&,*,  W*x) 

-  Al409^k*X  jrt* )  4  Aziotj  CWltx  f  Wa  )] 

■*  if  f  >4,«j  ( w*  <)  t/U**  ( *>*.  w* )] 

did  *  -  ,4  iw  <<,<?) 

c<  ,vo  ^  4Jj» 

"  f?  IAl»  (  V? ,VJ™ )  +  /^LIW  (V*} 

4  4 ally  (.^.if ,  Z/»  J  4  y^iOO  (W>,X,1fi  )] 
~  jj^r  Aj/wCw*,  YJhX  )  -  /4«w(*v»,,'iv^x  ) 
-Aui(tt)W*xj  )  iAziwM«M)) 

+%  (A  ,lw  ( U*  <x)  +  ALj  t  v?>  w'*ti 

( till.*)  =  A  loo  ( v^i.xi fj  Wi,T)  4/4  Aw  <  ,  tV*«  ) 

-*Aiv)(wSu, wT,x)  i A i(JtJ (*/?*, w2<x) 

"  2  C.  A  iov  (.  ^t,xx  >  Wi,x )  4 

■+  Aiwt.  M/*  >  v\//^x )  4  A <00  ( *C  >  >J 

,  I  *f  J 

vCxjk.jj}  -  Wa.x  )  -f-  A  xtv  (W/,]r,Wi,v() 

+  As{K)(rt*-.xx>Wi,rt  )  4  AiltofVi.x  ,^htx) 


~5  CAiOif  (  Wl,x»  >  Wi,x)  4  Ai(Kl(^''**M*,Xx) 


,*>  ..." 


n  .*» 


4  )  -f  AiUU  (.^i.x,  Wi.x*)] 


n  .» 


“  )  W,  )  i"  A  ZJ4H0  V^kx) 

Amoy  (»/£,)*/?* )  4  AiJtwi  W?, */”')] 

4  /^(AioyC^bytVhx)  +A*totVi*x,  V?’) 

4  AiooC  v?,vfi!)  4  A*i*)tv?>  w**  )J 

S  L 

4  TP  w"  ,20 1  Anoo(*l?>V~) 

4Aut*>  Ifi*  .^")i  /jwuj(V?,M,Zl 
'  4j«*j( -  Anat  itfiVh?) 

~  A  ink)  ( tfl > wT* i  **  Azi oqLV  *>  Wi"?)] 

~  ^  j?  (  Aijioy  OV,* ,  w,*)  4  vi/,*; 

^  /4  2Jt(kj(^K),,WiJ  f  Au/w  (w*jWltX/j 

-  &C Aw  (rt*,  v?) + AiM(v*t,v?) 

% 

^AioyC^i ,  irhX)  4  A* (V  ( Vk  >  v2/x)J 

C  Aj400  (  *'*  >  U*M  /4  J4oy  CM  /  VXit) 

* 

4  Aim)  (v£,  O  +  Ai^ivM 

"  Aji  oy  C  Mi,x  ,  If*  )  -  jjfly  (ft// ,  U,tX) 
“  >4  iwo  C  ^/x  >  ^  ^  ^  itlkj C^  *  W*,m)] 


(  =  '(%)(  A  3J WO  Wutj  t/,**)  4  A 23 it*)  (  Hi* •  */•'*  ) 

4  A  xjioy  ( ^ht }  w*  ;  -f  Azn  w  Otf »  ^x  Jj 

4  (  A\juq  (  2A,x » )  4  AtotjCUi,*  »V*  ) 

4  yji*  C  tt"  O  -t  At*)  ( V*  US) >J 

+  ffr  >i  w.»-  #**;  4  ^iw  c  w*  O 

4  if/TV  4 

^  j*(*j  (M*x >  “  /4  ( \Va"  14*"^ 

~  ^  X3(*/  (  2/y,  V\/^x  J  (.  Z/j  /^.x  )] 

~  3R*  C A 300  (  *W  J  4  (14,7,  lO 

4/4  3(jy(uTj  24>x)  4/4ioy(l/*>  Ui.x)] 

"  7^  f  ^/>x  '  24*)  +  A 3S(toW*,^l‘*) 

+ AJiw{xri,jttW,n)  +  AiiiftWi ‘Wt,*) 

-  ^  jaW(W,*,i/"j  -  24,xJ 

“  /4l  i  oc/  tt,x ,  Wa  )  ~  /4 13  cw  ( 24  /  Ut,x)J 

4  5”(^j  [  /4iJ*c*j(M'*>  W,  J  4  y4  ua  (»;  ( W*  V\Z,#)fj 
«  * 

4  /4  ijsf/y  ( » W4 ,)  4  AuwjW,*  Wt,x)J 


I  c  Ajm  ( 4  Ajtoo  i 

"  Aj4(X)  (  ^J,Ww  )  -  Aj4(V  '  Wl'*») 

+  A  J((K)(  j  )  4-  Aji (fij(.M/)j(  tVJi,x) 
-  A n  oo  ( v^x*  >  ^z**)  ~  *4  24  &>  ^  ,  vj,,t  >] 

k(AL(  v>T,.vC)  +  A!,',  ivr.  <’*) 

s  , 

~f‘Ant*>(  J  4  A400  (U*>Wt,tx) 

4~A  uk)  (  M/btx  )  -f/jtCWi.x  M.x  ) 

4  /44(*J  (k/i.*x»Ui  )  ’i’AuvoM*-*^**  ^ 

Jj?  C^J/w  ^2/  fc/ox  )+Ajiw(W*<  W/hn) 

-  y4j400  O^'.X  *  ** 

■f  /4 ji oc)  C ^htx )4  A iiou (,^**x*  ^hK ^ 

-  ~R  l  A  I  (*>  (  ^/,y  >Mi,x)  4  At  W  (  U,  ;  W//Jry ) 

4  ^4«j  c  is*,  x  >  ys».* )  4  Atm  ( vl  ,VJk*x) 

4  Ai^oy  ( Vv/ilif /  V?)  f  /3  ,ocj  ( > V/. x ) 

4  ^4tx;  ( ^*.w / 24* )  4A4OV  (■ 

*  5i?(  ^4i4w  (Wi,xx*Ui )  4 


|f  C  Aj 200  (  ■»  J  4  ^KKx) 

-  4  jjoo  t ^3  > M*)  ~  A  jj(fy( M,* ) 

-  aLk)  c^,»0 -aL^*^ 

4-  A  l23ut)(^i‘ix,)^1 )  4  Antk)  (M*»w/i,x  )J 

'  j?  ( A»w  ( ^3 » m!7;  t  Asw  t ^ ) 

4  /4»w  (  ^*>x  ■»  ^  4  >4*au  ( xjt  j 

4  Auk)  (Mxx >2/i  J  4 
4  A  3<*)W*'**>V »  4  4/4»fly(^X>^r^J 
~B  >V*CxJ  4^ 

” >4  J309  JJOU^' 

"  /J  7iCw(  ^xx » W,  ^  ^  ii0c>  ( M>,x>  W,A) 

4  /4i3W  (  wZxxrf?)  4  /4  wc*; ( Wi*^*x)) 
4  jq?  04jw  W*  j  ^/,x j  4 iU*,Wi,x*) 
i  A  1(KJ  C  U* ;  )  4  /4i{iy  ( 

4  Anv  (  ^hn>  U*  )  4  Wjw;  (t V/iX>24,x  ) 

4  4 joy  ( I^IJU )  l/*)  4  /4jV>  ( M,jc  >  V/"j] 


=  A  ioo(.  ,  v\//.T )  -f 2  wi!** ;  +  4  too 

■+  /4a  w  (  >  ^j.x  )  ■f’  -2/4w  (  ^,x*>  ^*m)  ~t"  Awo  ( W*.% ,  Mi.***) 

~  2  [A I (*)(  ^'***»  »  *^1.1 )  +  -24/CKj  (  y  wC,)  "M 1 ttj  (  W/,!J  ,  Wl*«  ) 

AiiwiWijut/V'Jt#)  ^  2^400  (Mi.iXtWi'ig)  f  Ail*)  (M/i,x  sWim)] 

,  j  it  w  .  i  ^ 

itxija)  *  /Uiw(W,  .*«,*/ 4,XJ  +  2Aluo(W*t,  ^i)4^^^vvljW,v^Jj 

+  ^3  (W  (Wi.ttX  y  )  +  2  /Cc*,  (  i  Aw(*t?lK ,  WtX'x) 

"  xC/)  JW  (  ^/,ax>  M*.x  )  +  2  >4 i  <*3  ( M/XX  >  V\A, X*)  ‘f /4  wo  ( «*) 

+  4jty^,'w'^'*)  +  2  >4  3(*,  ( M** ,  w,*,  j  ^  m*  .vVmojJ  (>4-/74) 

I  I  71  At)  -  l  y(  ■ 

+ §>M  I* w  ( Wi ,  UD  *3  /4  ,‘w(  W* .  2^»  J  +  Am  ( W*<  Wjnr  j 

■*  A  Uv>l,  v?)+2  aL  «*2 .  «*”)  +4«  «£  Cj] 

+  >4  Wj  ( z& >  *0  1 2A™P&  w™)  iA^  ( i£  <  j 

-4W«A )  u,")  -2A‘mwZ  y  u,T)  -4«K  C; 


"  /4  r»W  (^>/>x  j  Wa  J  '2Azi(ig(Ui,x » ^,x)'>4iftei(^  Wi.x*  )J 
~  J"?iP  (4  I3409^'*,K>^'  i ^  ^4mw(l'/i!xy^,x)^4lJ46y(^|,,»  ^i,w) 
4  «  4ijKkM*  .  h/^j  f  /J  i 
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”  5r*  r  All*J  (  Ulltt,Ut  )  4  2  Al(t)(Vl,X  •  Vl.x  )~^A KiQlJfi*  ul*jx  ) 

+  Aw  ( a*w ,  if;  +  2  ^  ,  u* )  M  U,  ( i£  wJ)J 

~  ^  ( >4 J4(jt/  w'.*»  '&")  +MjiOg(  W,*  ,  u*x )  +  Autt)(W*  viJU) 

"t  AlQfH)  (&J,*x  ,Wi  )+ ^Al400  (  M*.x  ,W,,x)’i‘  Az4(y(U*  •  ^>m) 

-  A ji oo  (  m*-** •  is, )- 2Aj, m( n/3 , z/,* )  -/4 ji oo ( ) 
-AimlvS*#™)  ~ 2 Ait hUtZm") -AiatflS* V>%)] 

~  ~in)  )n 

A  UiUtjW'-**  )  +2  AxjiOO  (  ,Yji,x  )  4  A  X3l9)kfiJ?MlM)] 

2Jt.x  )  4  Aity^U*  >  Utxx) 

4  AnK){Pi.n2JT)  1 2A*ik)(U,x >Ui,x )  4  AiW(  If, ,  Ui.tx) J 

W.x  )+^j3(n)^<  M.u) 

AuwiUi.tt,  M  )4  ^^iJwCL/?,x,lV/,x;^xi(*)(^n! 

~  A  JLi(K)(v^,.*x»Vi  )  ~-2  Aj2iii)[Wt*,Ut,t  )~Aj2w(W*,Ui,tx) 

~Ai3w(U*,xxM*  )  ~2  Aimi^x,  ^4,xJ-Ax30tf(lfi'M,j(>)] 

4  dfiftAzjwto.*  .  Jt^AijtotfiWixjIA/i.x)  4  A  hi  (*j  (W*,  iVi ,  o  ) 

<4  w3w  (  a//,w  j  4  2  4  >4  wtf.*u )J 

"  w.*;  4  v~?,u,,*x) 

4  Ah*)(Vi>xx  tlf*)  +  2A*ut)Wt‘*,V*  *)  ^Auki^i  >l/i.xx)] 


4  4  L.(  w.«  ,iv*j  +^<«( v"  ^ 


-  Aji(AjW***>V?)  -2  Ajm)M** *  ^ 

-  /4li (*>(^.XX  *  ^  ^  ">2  ^XJ(W( y»!x  I  */■>,*)  “/4ww^"'  ^  A'1 7  ^ 

(  T^xy/ijor)  =  Tyf  /4  •  ^4.x)  4.2/4  Jjbcj  (Wx,x  >  ^'-wj  4  4ji(jy(^J  Mi.w) 

~  >4  W-*.X  )  4^^  j4(jy(tV/jJ(,lV^w)  ~AjU(t}(^i  •  WJj.wt) 


4  /4l/0g  ^»,JU0r  »^i  )  i ^Aiuv(Mi*x* M*x  ) 


*•  /.,,**  ^  ,...*  .  ,*+*> 


-  /1x4W< ^ ^/4l40^4'** > >VJi.u)J 

**  ^  (  4 1  oy  (  i//!xx  ,  W/^Tx  )  4-2/4  iwlVl*  *  ^,'a  )  4  Ai  w(  if*  ^M<x*  ) 

4^^(jy(Ui.xx  ,  Wa,y  )i2A<WM.„K,„)-t 


n  .  **i 


4  AmM***>Vr')  4  )  +AU*'  >V,«) 


4 2  i2  AAoq(*U.»*>Vi.K  ^  4^4W( towUHn)] 
jk ^  A  jIU()( 4.2  ^jiwK!  ,  w ,"*)  4  4xi(Wj*  w/,^; 

“  Ajwj  >4 jw M»x,  Vi,** )  4 ^J4fu(^>  » W^w0 

4  ^  ii(fi)(.WuBo,,Wi )  -f2 ^iuw(M»w  >Wt,*)i  Anv*M'‘*>^*>w) 


125 


4  ^  (  A  m  (U  ** ,  vV/.T )  4-?4  U  C  U.7 » M,xx )  4  4,  ^(  U*!  ^,x«) 


4  Auw  (u*x  I  vV*J  4  2^4w(  ^ '  *^0  4/44(x)^^ 


^  i  i  ^  >1  *■  /■>  u  /  ^ 


4  Alt*)  (  Vi  )  4 -2  4  ,  yy  (  V<,«  »  )  4  A/ (/cj^hx  '^'*0 


1  /.  ,*  ,r"*  .  *i  Of"' 


4 4<*(« (IA/j,wx»Uj  )  4 2  Aqua OMi.x* , Ui,k )  4  AjfwCWt.x tUi.x*)'] 

||i| 

( 'tyiin)  ~  %  r 4  (Aj(  v^1‘xx  •  Wj,x )42 A  4 >4  jiw 0^*  > 


4  J3W  ( *  ^/.x ')  ~2  AijmMi*  ,Wt.xx )  ~Aji  uyM * 


"  AziVi)Wi,m,W,  )  - 2Aiw\(Wim ,^hx) -Ai2M(W* >Wi,x») 


4  A u do ^■l<«x*v'4  )  4^1JW(IV1A>,Vy/j)^)  4  4n(x;(W*x>Wi,x)t>)J 

~  #  f  4s(w  (Wi3x» W4.x  )  424jw (ZT,*,\a/^)  4  4ww(^VV*« ) 


4  Aiucfa'**-*  ^»X  )  4  -2  4^|,'i^* ,  Vy/j^)  44j  w  ( WiT » 

4 “  /l2(*j(^l,X<X»2£*  )  t2All*AWl,M)it}fj'X)'l'Aiw(Ni,x,Vi.K,) 


4  /4j(jo )  4 4<4j(K>^a-x,^,xx  1) 

“ip  C  4j2(k/V'^*x  »  ^x)  4^4  jhkjWi*  ,Wl,w)'l'Aj2(rt(Wi 

~  4t3K)£M'x*  » V\^x) -2 Wi,xx^  “* 4j3<w(^'  »W/,x«) 


”  /\z20qWi‘***  ,  )  -2  AxiU<)(W  h*X  fVJi-x)  ~Al2(Jt)(Wi'X'*/t.*x) 

+  A  73 w  ( ^i?xw,  vy/*)  ■+  2  rfnoyM-**  >  Mkx )  +  <4  iiuj  ( 

+  §'  (/4,oy(  ^<XX,V\A*x)  +  2/4  3(*J  (  U>*  ,  \a/^)  +  /)  i(XJ  (  Un  •  ^'W*) 

*t  /^i(Kj  ( Vim ,  Wj,x )  +  2  >4i(«(  l/; ,x  >vJmx)  ■+  /4i(w(2/»* ^v»*) 

i  AmWw'V*)  +  24jm(MZ.l/S)+A^MJi’tf£)] 


(4 -/7£> 


where  m  is  the  number  of  the  iteration  step. 

Substitution  of  Eqs  A-161-  A-178into  Eqs  A-149  and  A-150  one  may  obtain 
the  iteration  equations  for  the  nonlinear  part  of  the  stress  and  moment  re¬ 
sultant  vectors  ^  \n2^i)'  In  80  doinS»  new  symbols  are  Introduced 

and  defined.  The  part  of  the  t's  or  nn*s  that  is  linearized  (linear)  with 

respect  to  the  Iteration  parameters  (containing  u^1,  v™*1  &  w”*1)  is  denoted 

!  nL 

by  superscript  L  next  to  n,  i.e.l  tj^  .  The  part  that  only  depends  on  the 
value  of  the  parameters  at  the  previous  step  (um,  vm,  wm),  is  denoted  by 
superscript  n  next  to  n,  i.e.  \tl^J . 

^  ®  1  /  f  /  *M*+I  f  L*n\n  \ 

(Ha)  +ftcj  ) 


K  = 


B  D 


V!  - 

A  B 

rf  - 

8  D 

\  4 

(»i  -  w 


(A-iJV 


(A-iSo) 
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lA-iti) 


J 


,wl 


(A-ite) 


c  a  -m) 


In  a  very  similar  manner,  the  nonlinear  terms  of  the  equilibrium  equations  are 
also  linearized  by  Newton’s  method: 


~  7P<)l(  ftlttli  t  Wmi)  +  J  Wjyt  )  f  ,7tyri) 

'  4m«j(  Wi ,  ?t„i) - AituAwr'.Tl^, )-A}ilMAw?.%n>)+Aj>wW**,’r‘j) 

-  &iAUv"X')+AUrf'.’C) +A'Uvr.7i;,) 
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Auw  ( V*  .  7lfti )  i  Auih  (  Tlyyi )  -  Aiuo  ( U*‘  7?m)] 

^~R  [  At iw^1'*  >  V.t'fx)  *f  Ai(n)i^i.x,^lxyl )  -  Ahk)  W*,7byi) 

*  ^4WC>  *■  ^  ^  A400  (wll » Tim)  -  A  ( W4x ,  7?** jJ  ( » - 

=  ~7?(T?> yn*  TlMl)  y  4  7lmiiX  + 

7Z$t{  W-yiii  )  4  -—?-[ y) jw; ( ^4  ,  +4n«6v" ff™) 

~*  >4ji(u(iV,*#w)J  -  ^[A3ucjivT\  Tim)  iAsutji  U?,7l™i) 

-aUvTmS,)  -i-AUur.Ti;,)  tAiuo  wrxzuuv'Kfi  m 

/?*  f  4jia«  »  Tty,,)  <W/,  -+  >4 i, IV,** >>^1^ 

-AitwMxS)  -Aj4W9i<*?j-4^w:x,i) 

-  Ttyy, )  ~t  An  t«y  £  ^  <  #yyi  )J 

4  §[Au*j  ,  9tw)  +A,LUX)CW,l{  Am)i  A,\uo^> Z ,  7lw ) 

Amuo  1  >ftvt)  -M 4I1AJ ( *  flip )] 

«w,T)  -/^jw  (V, ■",»-) +^w; 7jj 
-{■l  A  n»,  c  vr,n£  >  i  4n»  ( v',^j  -  aLau:.  X) 

*&•»>  iv?>*£)  uU  <£.  o 


-Hi) 


+%  1 A  L  < W,7.  nz, )  +  /l  L,  ( <  Xu>  -  /!  ,U  ( <  ■  <u  > 
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4  Aiuw(  -%yi)  4  tju)J  (.A~f$7) 

^>2bi  ~  ~R*  C  ~  A  J3(JKJ  (  W,*,  7?yy,  )  -  /4j3(Jjy  (Wl0,  ??>/,)  -  4jj(JXJ  (  'Nf>)7lyj> ,  ) 

4  Aj2ukj  (.  V\f?,7lyyi  )  4  4  4 ^uivjW1  >%yi)] 

-  %  (  A  i>*i  ( V?,  %yyi)  tAu*j(V?, 

4  £  [  An*)  ( ,  T/jjr/)  4  >M*yj)  i  A}tiy(Wi,x,?tvi  ) 

4  /I  j  (  A/,.x  jTlxyz  )  4  /I  »«*;  ( ^i,x  >  7l*yi)  4  4  2(«;  ( M/yX , 7?xyt)J 

4  |p  [  ~  4  jlt*y  (  W,*7  )  -Ajh*)W?'  )  4 4w t*j  (  W*  ??/*J 

4  >  Wyyi)  4  /jjiCiA;  »^yyi  )  -  A  JH1KJ  (W*ttyyi)J 

"  ~fcx  Miwsi  (  ^  '  >  %/i)  ^Axuxj^A>^yyj  )  ~~Azu*j  Wt  rfyyi) 

4  aL)  ( 2^:  JO  4  - 4A*  (^,^7;; 

4  A?  C  Aj(Ucj(  Wi,y  ,  Tlt/t)  -f  /liCiAj  IV/*X  J?liy,')  ~  Aj(i*j  (Wt*  t7?v/J 

4  4  itutj  (  ^Ixyl  )  4  4-2UO  0V(7  ,  %}■!  )  ~  AllW  ( K//.X  >?7*)U  )]  (A  -m) 


ft*  =  #fc?/U  +*C»  4  2  c4?3  («va\.x  4  ^i.y  )  -  4  *Cl J 

-Riftyyk  4  ?lyyn)  4  ^  (Ajity  ( W»  >7l*y\.x)  44j/w('v* 

~  Aji(V  W*  }  Axy}  x )  "4j  4^;  (  w,  \?lxyilX  )~A  liAjilwiW^Ww.KjJ 

4(?)  M  XJI  wWfxtyy*)  4  4ijW  ( w/,"*  J  )  ~  /4w«W  f  V^>^yy,  )J 
4  (W,  '  >W  +4i34vW*X»i)'-dxxw  (W,*7tZ)J 

4  R  A tw  (  ^  '  Wx/i,x)  -  A  it*.)  (U*,  Wxyi,* )  4 /I  i\*;(  Z/-7 
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"  Aaou  ( ~  Atoy  (■$?•  Wxyj.x)  'bAi(it)(Ui't??xyi.x)J 

+  Au»  (if!  1&)  f  Aiu>lU?.*Z)  -  Alm(vr.*m) 

+  Am  ( <«X)  +A  miwZ.  <')  -  AmjCW'Z  J& ) 

"M^W  Wi,xx  ,  7lxi (2)  +/4»W  Wi,**  ,  %ij  )  *  ?l*L  )J 

4  T4  J»(«  (  Wi,x  i  )  ~t  AjUkJ  (  TJi,%  yTlxyi  )  -  4  J| w  (  V\/i?x  ,  %?I ) 

~  Ajiity  (  IA A,*  ,9lty))  ~  Ajtf  (*)  >7lx/i  )  i  Ajuo cl  ( IA/|!5  ,Tlxyi 

Ajuou^1  *  Myyi)  -  A  auk)  (^1  , T^yyiJ  +  Aji/cv  C'^'*  Wyyi ) 

A  Ji*w  (Wi  >Tlyyi)  "  ^JWWcjClVx  >Tlyyi  )  'hA  ZiU(k)  U 

^  Rl  £  AjiwiVi  •%#()“ Aj  100 (V*  >%w )  ^  AsiwiVf’ftjwJ 

4AjhUc}LUi  >7lyy>)  +Aj4wLV*>??#i)-A3ii(v(V'*>Wyyi)]  {A-iW) 

=  *W,x*  -f  ^x«i/XX  -2  L%)(M?yn,i+ ftwi)) 

-<#/<*& +*2u)-ii*Z+nZ) 

~*(%)CAj2w(wrl,?ixyi.)i)-i-  ASm  (w?  ,Tt?y'ix)  -  Aj2ud  w?.Hw,x) 

~  A w<  * %/«)  - /4j3fw ~t  Ajm (w" 

4  ^  ^  “  >4  a  «&;  (Wx  >7lyyi)  '"AijiwW*  >Wyy\ )  A  zj iw^* >Wyyi) 

-  A1J Jiw < w,  -  /4 1«w  ( ha"  -f  4  i JiW  <*T. ?C Jj 

4  >W*y2,x)  ~Aiw  ( &*>Ttxyj,x)  "f’Aiw  •%*y*.'t) 


“  Ai  uv  (V*  '  ytxyitx)  ~  Aim) (  vT,7ltyi.x )  +  As<*j  ( V*.  %yi.x )] 

+  ]?[  A 1 2ikj  (.U>  >  W»i)  +  A  uov  ( &  >  ftyyi )  -  Azov  ) 

■+  A  H{K)Wi  > 72 yyj  )  A  AiitoiQ  >  W-iyi  )  —Ali(K)  )] 

A  /^2w(lV/,xx  i Tlxxi  )  A  A*m)  (Wi,m  i  VI w  )  -  Ai(ju(^h*x .  Jlxxx) 

A  Auto  (Wi.xx  •71v*l)  AA3(*)(W*.xx  >7ltx\  )  -  A3(k)M*!*x>?tiui) 

A  lAjiuo^Wi.x^xyAAAjjuioiWi.xjflfyij-AjimjMi.xjTlxy) 

~  ^  J30CJ  (.W'.x  >  ?lxj\  )  -  Aji(K)  Xx  >  %yi  )  A  A  JifaM** 


4  -§  i>  /O  i£  Xn)~  Am  ( vZ ,  n£)  aaL»  ivz  .  n?y> ) 

I  J  W  htti  U  ygm 

~  A  3W  (  V>,x  J  7lyy[ )  -  A  J(*J  %*yi  )  A  AiLkj  %*//)] 

4  (% )  C  AjixfkJ  C  w!* Tlyyi )  -  \ y\yn)  i  Ajnw ( W  *»  9?^) 

~  A  jisw  ( W*  >  %ty\)  -  A  J2iw  (vs/i*1 [  ft* , )  f  Ajuotj  ( *  %yyi)] 


J2(*)  1 


4  ASmilC.nZ,)  A  Ajm(V?,n£)  -AjmW'ffiZ,)}  (A-IVO) 
^un  ~  ~r  C  AjmyWi,  Tlxyi.x)  AAjnmji^',  ^lxy,.x)A  Aja^i^^Xx) 

~  Ajnllk){.W?,72A)niX)  -AjuwCW,,  Tliyj.xj-AjtnivW,0, 


~  2 Aj4 (hoLW^x  ,7Ltyi)  - 2  Aj*u*)LvJ°x  Jltyx)  ~3Aj<u*j J 
A  Ajl(k)(-Wi  iV.xy),x)  A  Ajtw(W*  ,72xyl,r)  (  W* ,7Uyi,t) 

~~  AjUW  yix)  ~Ajm(^l  tfloi.X  )  A  A J4()ej  W,  ,7Uyi.t) 
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+  2  AjIOCJ  ,  7Ul\  J  2  A  JUK)  (  Wj.A  >  fttyl  )  ~  2  Ajw(Wl,x  >%*yi) 
~2  Aj4(JA)(  ^*X'ft*n)~^Aj4(iV^,x,%xyj  Ajtuik) ^Xju )] 
*K$J  i  AlJU2K)tW>  Jtyyi)  (  Wj°,  yiiy)  )-f  A IJ I  OK)  (• ) 

■f  A IJ462W  ^?)jr 3  "M  A  1J4qx/N>°>  7Iyyt  ) 

-  A  JJIOAJ  (W I* >  Tlyyi  )  ~A  JltUtyCW?  ftyyi  )  ~AjiiUioCW,  >7?yyi  ) 

-  J.24(<ytf  (Wt  'Myy)  )  ~  A  JiiUik)^1  >??yyi)~  AjiUtUO  (W»  iftyyj) 

+  A  UUiK)  (W™1*  Tlfjx  )  t  A  UlLiK)  )  ~  AlJiutoWitflyyi) 

^  A  u4 (in)  ( *V»  /  ^  ^  X340K)  ( W/  » ??yyj  )  -  A  moy  (IV/T  ?/>yi  ) 

~  >4jafW  (  W/  I  Tlyyl)  ’'AlXuiKjV'^1  >%M  )  iAj>U*>W*  Wyyi) 

■iA  won) - Amm^r.Tly™')  •+Aj#k#)(h"%»)] 

47?(~Ai ilk)  ( 7lmx )~4iuaj( 2/*,  7l*?i,x)- A ,\nv  ( U**  > %*») 

"  A 4WJ  C  Vi, X  1  Tlxyj  )  ~  Anijiji.  V? \  ??tyhx )  ~  A  Uik)  (  V*,  ??*yl,x  ) 

-f  Auik)  %W*)  "  A 4iimj(Vi,a  ,p?xyi)  '-4<U**A^i  ’  ^*yi) 

-f  aLj  (V?,  ri")-Au«)  ( UZ\  n%)  -  A1U*)(V?*  >riw) 

i’A  Umj  (Vi,x  >%yi)  ~/^4(«U  (  Vj,x  ,  %yj)  ~  A  4  L*k)(V>,X.  >?hy>  ) 

^  A  40*)  (V?>X  ,  Tlxyi  )] 

-*  ^[-aL*j  (vr,n£)+Ai4**)  (V?>  Kl)  -AjK*>(v?>rf') 
-tAjUW(vr,?i£) -a  i  l(M<vr\n£)-Aiu*,  c  uvC> 


4  ^i<K4K)  4  4i4fc*J  ^  ^ 

4  Ai<u*0?>ri?i)  -A^0"[nZ)-A^m0r,nL^) 

4  ^j/(«j  (&  > ??yy/ )  4  4j4ojy  (W  , )  4  >4 j«(*>  (Z//*  ?Zyyi  j 

4  4j4ti*/  ( ^yy-*  )J 

4/4<o/o  OV/!xx  *^xif|  J  4  4/(40  4/^Kajy  Wi'txJVti) 

~tAtiixK)  (W***  >7?*u)  +  A m&)  M*.n  »7lm)  -bAuM  (W!mx 

4  4 1 au  C^i,x»  t  ^Ixxi )  4  4  /  (x*y  ( >V<,)u  > ^  ~  4»t**j  (vv/mcx *  J 
4  44(U0  W-Mx  j  f  4  4 (0/0  (M'*x  •  %xx>  )  -  At(W  Ww*  •%*«  )  (A'tf  I) 

$u*  ~  j*r  [  Ajiuk)(  ^  )  4/4jiuxj(v^i  >  )  4 4  J  i (wj  ( W* ,  7fo«.0 

-  4j3t^  (w*%y/,<)  ^JJc^WA^xy/.x)-"  4jWm'.  ,?7 »y/.*j 
*+«2  4  J2Ca*)lW£x  ,  7l*yt)'t2A  JUikpti.X'TltyiH  lAjaugiullA’Myj) 

-^4 1 -zAjiiuyW1'*  •  Mw)  ~2A  wcuot^*  flx/i) 

4  4  J  2£h>(^  >  4  4ji(ycj(  IV*>  7lx>i,*)  ~  '  'Wty*'*) 


-  4j3(*>(*A  >%yi,*)  -4jjo«(W*  %wx)  4^jjowW 


+2/hU*?*-*S>)  ^aLhJWZ-Xv*)-^*^'7^) 


'2Anm(<VM,7lw) 

+  ($)  { ~  Ajjicw  )  ~  Auxij(Ml  ’flyyi) 

"  4ij3(ixj  iw,  »  ^yyi  j  ~  4jj3(yc;  (VA/i°,  7l^t  )  “4wi«4  (W»  >  ^yyj ) 
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>^»i)  ~Ajhuio  >7lm) 

~~  A  ja nuo  W*  >  ^>yi  J  “  ^  ja^tuo^i  *^yyi )  ~  A  jmnoi^*  *^yyi ) 

”  ^IJiWcJ  (wn  ??yy/)  “  ^xja(«j('VAW.^?y/i  )  "  /^WJttifJ^* 

“  /4 i3iv*)  i  W|  >  ^yya)  “  »^yya  ) 4‘4ij3uw^'  • 

”  /Ajitoti)  (*V*  1  /^yyaj  “  Ajnwji  hf,  ,  TZyyj  )i  A  jn<.uQ(W?,4lyy3) 

~  ,  ?ZyyJ  ~  A  jJiw ( W?,7l„ ,)  -f  Aj&mW*  > ^yyij] 

(~  Aiay  (  Z/<  ;  92xyz,t)  -  mo  (  ^4*  >  9lxy*.x)  ~A  iu*J  C  2//*  ??>0«»x  ) 

4  A  uvj  C  2T(  >  Tfryv)  *~/4i  wj ( Z/i ,  %y/,  x>)  ~A)liK)  l 2^1 7?*y**) 

-  (ZTt  i  #*y/,y)  +Aniy(U*>  ftw,x) 


-  Ai<j*j(&t*’Mw)  -  AutvC^ ~  4  ,%*/?) 

4  /4itlKJ  (U/,*,  7?*#  )  -  ^  3<iW  ( lf*t  ,  92xyi  )  -  AiVUgiU*,*  ,  Tlxyi ) 

"  AjUxjWJt,*  ,%xy)  )  4  A  j  ia*;  (  V*.x  >?lxy\ )J 

Aij(jtj  (  IT*  ?7yyi  j  4  /4laUW  (  ^  >  %W-)  4  ( IT,"  2?yy/  ) 

-  4  x«W  (■  W*  ^/y;  )  "Aui**J  (  *  fyyj  ^  “Aum^Lli T#yy4) 

-  Ano*)(jX*  ,7lm  )  4  Wyy*)] 

4  A  -i«W  (  Wbxyxflfti  H^auw  (Wt.xx  tAtn  )  4 


A  3  l»J  (  ^-xx  » Wxtt  )  4  (  M,** »  Wxxi  )  4  Am*K)  (  >7lxx\) 

4  AwitiiK.to  ,?bai)  4 /\uw(.W*»  ,71m)  “  Aittyf^/'JucWxKi) 

4  Auik)  (  Wi.jo, , ?£**, )  4  /jitJW ( lVa,x* , 91  xx i  )  - 4jW(  *y*!w , )  (A'Wi) 
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After  linearization  the  equilibrium  equations,  Eqs  a-160 ,  can  be  written  in 
matrix  form 

[Cuj(n,« I  +  cc„j[«m]  +  rc.]U]+t  mnL\ +CE„](«i,} 

+  c  E, ']{%„]  4  C8„)[cj  HCu]{#,m | 

4  c +  [C,]{%J4(E„]{^1„}  +f£j{^]  +  (Ej{«i,' 

<  ao  [nZ J  +  c  ajf^j  +£  ej^J  4  (4,0  {xj  4  c A,)(X, 

+  CA.][Xj  =  [#J  (-*-/?*) 


where 


{^•j  ~  L  »  ftyylij  ^X/ii.  'Mxxi .^tyyn,  tyxyilj 

f^1}  '  L  1  ?lyyii ,  ^lyyii  >  n„>, .  7%xyiiJ 

[  ~  l  2X 1  i ,  £//  l  ,  iV/i ,  2^21,  Wii ,  Wii  j 

f  „i)T  L  L  L  1  L  l  '**** 

>  ftyyii  .%wu  /fthxu ,  fflyyn  ,Wwij 
[^ij  =  iP^txn  »  Tlyyii  >  flxyii  fflxm  >  'ftyyiit7ftxyi'i j 

[nf^  jiZ ,  , 9£tl.*£*.*G» 

f<f*  u  C.  .7C.riS1.vC2 .JC.jCu*' 

In  Eqs  A-145  A-150  {*},  {t1}  and  {tn}  can  be  written  as: 

[6,i}=[kij{X„jf[kl.]{x) 
le.i}*[Jaj(X.<}  +  W‘){x} 

tf}<*fflxj+[ki:)  (xj 
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c 


c 


c 


e 


r 


i  i'j  -  (xi:jfx,j+f«:jw 

{*;*}■  c*/»]Du}  Hti'nUU +u/.*,)(xj 

{t*h[xu]{X,»}  H*£]{X,«)+(  k  I ,*«]  [X.,\ 4CK  Cj  (xj 

(ti jmj'  [kl lit]  (X.urj  t[K-21„]fX,x*j4  +fr/..][Xj  ( A-l% ) 


Substitute  of  Eqs  A-145- A-150,  and  A-194  into  Eq  A-193  yields  a  matrix  equation 
which  only  contain  the  vector  of  unknown,  {x} 

fR«{x4  +f«3]fx,»,j+cw]^j+rRf]{x.j+r??o]/jfj=f?j  (a-hs) 


As  in  the  case  of  W-F  formulation  transformation  equation  are  introduced  in 
order  to  reduce  the  order  of  the  linear  equations. 


{^j  ~(X‘Xx\ 

By  this  transformation,  Eq  A-195  can  be  written  in  the  fbllowing  form: 

OTft)*w(a*w(;j-/6i 

A. 3. 5  Boundary  Condition 

Boundary  condition  A1117  can  be  presented  in  the  following  form 


Either 

Nxx  =  A/xx 
A/xy  *  /Vxy 

Q*  =  5  4  AW 
Mx*  =  Fiv 


Or 

=  Conat . 
y  -  Onut . 

0 

lV.x=0  (^-/I7ft) 


f 
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Q  —  Nxx(V^iX  tMx)  4  AJxy  (VS/.yfwAy) 
4  Mxx,x  4?A7Xy>y 


Obviously,  the  boundary  condition  can  be  written  in  matrix  form  (at  x  =  0,  L) 


[121] 


+  [Al] 


=  H 


(A-i  n. 


where  the  form  of  [-OJ]  and  [xfl  depends  on  the  type  of  boundary  conditions. 

The  stress  and  moment  results,  and  the  displacements  are  represented  in  series 


form. 


£[ T&Oityi i  UnAtn^) 

C-fiO  . 

£(QCe?Tr+Q  I  Z.  +  VJaMnn$] 


( Wn, xc#n( + 


After  applying  the  Galerkin  Procedure,  the  boundary  conditions  can  be  written 


as: 


Olxxi 

’  U,i 

—  *1 

iTa 

??*yi 

~*i 

Ok 

1*Al 

W.«i 

'  +  f/V]  - 

VJa.x 

—  1 

/liyL 

V« 

OOvyi 

Q 

Wit 

Oft  mi 

lA/il,)t 

* 

=  f-e] 


i-  o."'  .4 


where 


—  i 

01 


n1 


XXL 

**i 

xyl 


ar 


—  M2 

Qi 


onU 

Oftxti 


=  OOxya  4  Otxxii  4  otxxii 

-  7ltxii  t  Olxxii  t  O^xxu 

-  Kli  +  t  +4lWw  +  *t* i+Ka) 

*  4  0?xyiL  f  flxyii -t -jijr  (fllxytit  fflxyii +Wxy*i-) 

•f  (WitlVi"»  %ry)  “*/4  jmwi 

"  ^  (A  it*>  (•  V*  >  Otxy)  4  (no  ( ^  >  ^*>-0 

=  Oftxxi.x  '~*2-1*fe7Kxy  4  /Umv  (Wi.v+W,*  ,71*)+ Ailw0rf*.t+W>*t, 
% 

4  #  £  -4  j,  «*,(*.♦*.•.  *i;J 

“"■^r^iCugCyoTtxyJ  4  /43fKj(Z/i,7?Xy)3 
~  fflxxti  4  Mali  f  OPIxul, 

=  OJlttii  4  OHkxiI  4  Oflxxti 


(A  -100, 


xx  ) 


n») 


(A -20 


Using  the  similar  procedure  as  used  in  sectionll,  Eqs  a-200  can  be  linearized 
and  written  in  matrix  form: 

4  (AJ  jx‘j  =  (AJ(/M,1 

rAj(fNlj  +  fA/,‘i+/’N,“}J  +  CA]fX‘j  *  M-CAlM  (A-iot) 


Substituting  of  Eqs  A-145-A-150,  A-194  into  Eqs  a-202  yields  the  following  form 
for  the  boundary  conditions 


--  M 


(4 -2  03; 


I 


I 


m 


<• 


A. 3.6  Solution  Methodology  -  Finite  Difference  Equations 

The  linearized  iteration  equations  (equilibrium)  assume  the  form 


<A-iU) 


Note  that  the  true  number  of  unknown  is  (6k  +  3) .  These  are  u, . .  v  ,  w, . 

li’  li  li 

(i  =  l,2...k)  and  u 2i,  w2i  (i  =  l,2...k)  [see  Eqs  (119)].  For  conveni¬ 

ence  though  the  number  of  unknown  is  treated  as  (6k  +  6)  with  u2q»  v2o  &  w20 
existing  for  the  count,  but  subject  to  the  constraint  u^  =  =  w2Q  =  0. 

Thus  with  the  transformation,  vT|j  =  \  X.  _  J ,  the  number  of  unknowns  is  (12k  +  12). 

The  equilibrium  equation,  Eqs  A-196,  are  next  cast  into  finite  difference 
form,  by  employing  the  usual  central  difference  formula.  Thus  at  each  node 
point  j,  the  equations  become  (in  matrix  form) 


(x>o +  &sf) 


(A -104) 


At  each  end  one  fictitious  point  is  used.  This  requires  (12k  +  12) 
additional  equations  at  each  end  (j  *  1  and  N;  the  fictitious  points  are  de¬ 
noted  by  j  *  0  and  j  ■  N  +  1) .  These  additional  equations  come  from  the 
boundary  conditions. 


9 
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Paradoxically,  the  number  of  boundary  equations  is  (8k  +  8)  at  each 
end.  Note  that  these  are  either  natural  (8k  +  8  through  the  Galerkin  (pro- 


4 

t-.- 


S 


cedure)  or  kinematic  (8k  +  8,  u^  =  U£^  =  0,  =  0,  w^  =  =  0  & 

w^  =  =  0  for  i  =  0,l,2...k).  This  necessitates  the  requirement  of 

’x  *x 

(4k  +4)  additional  conditions  at  each  boundary. 

The  additional  boundary  terms  are  given  below  and  they  only  involve 

uli  ’  U2i  »  vli  »  v2i  at  eac^  boundary.  Their  existence  deriva- 

’xx  ’xx  ’ XX  ’xx 

tives  with  respect  to  x  of  the  displacement  components  u  and  v  in  the  equi¬ 

librium  equations.  On  the  other  hand,  regardless  of  whether  or  not  the 
boundary  conditions  are  natural  or  kinematic,  they  do  not  contain  second  de¬ 
rivatives  of  u  and  v  with  respect  to  x. 


£  1)  + 1 

\  wu.(3  *1  lh 0 


I.A-20S) 


142 


Where  the  constant  1^,  T]2  and  T|  are  assigned  to  achieve  certain  goals 
(in  generating  some  results  =  1  T|2  =  -2  and  ff3  =  1  are  used,  which  implies 
that  a  derivative  at  a  boundary  is  obtained  in  a  forward  manner). 

Note  that  Eqs  a~205  are  the  additional  (4k  +4)  boundary  terms  and 
that  these  equations  are  incorporated  in  the  matrix  form  shown  in  Eqs  A-203. 
This  means  that  [db]  and  [dc]  are  square  matrices,  [(12k  +  12)  by  (12k  +  12)]. 
These  boundary  equations,  Eqs  a-203,  are  also  cast  into  finite  difference  form. 


I 


-f  iva 


where  j  =  1  or  N. 

A. 3. 7  Total  Potential  &  End  Shortening 

The  expression  for  the  total  potential  for  a  supported  (ss-i,  cc-i) 
cylindrical  shell  is  given  by 

/  /i*£/ L  f  o  0 

“  2  jo  J0  C  +  N  yy  € yy  ^  A1 xy  Ifxy 

■+  MxX  k*x  +  Myy  Kyy  +2Mxy  Hxy]  dXcly 

•+J,  C -N„u-Nvv+Ffx*M.x}l,  Ay 

'  j.  J  Jwdxdy  {A-20D 


or 
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-  ~2  f0  9lxxo  £.**»  4  €yy°  ^  fcy®  ^  ^K*° 

4  ^xy*  +  -2  #xyo  A'xy®  4fjL  (  71ml  &ml  4  ?2yyi  £j 

i*  i 

4  ^xyt  foi  4  ??jg(i  6xXi  4  ??yyi  6yyi  4  fttyi  Yxfi 
4  /%rxi  /Gxi  4  ^Tlyyi  Km  4  297lx^iKtfL  4  ^x*i  /Gt, 
4  TtyyiXyyi  +  -2  971  tyi  K.xg  £  )]  ^ 

4  *27£/2  (-  7lx%  0  U«  4  %(*»  l/o  4  ^x^e 

4  WxL  W'  -  ^xxo  ,) 


«i 


yt 

yiy 


Ui 


4 


=  Ml  4  M.  M 


c 


€y/i 


=  fe,;}  +  {t‘}  +  {taj 


| 

l2>Cvi. 


K* 

.-?Xxii 


=  *(£’ j 


and  n^\,  n1  ,  u*  ,  v]^,  wf^are  the  values  at  x  ■  i,  n*°  ,  n10  ,  m10  , 

“xxi’  xxi  xxi 5  1  1  *  i  xxi  xxi  xxi 


10  10 
V  V  i 


are  the  values  at  x 
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APPENDIX  B 


COMPUTER  PROGRAM 

B.l  w,  F-Formulation 

B.2  u,  v,  w- Formal at ion 

Flow  charts  and  program  listing,  for  both  formulations, 
will  be  made  available  upon  request.  (Write  to  Professor 
G.  J.  Simitses) . 


APPENDIX  C 


MODIFICATION  AND  GENERALIZATION 
OF  POTTER'S  METHOD. 

The  behavior  of  several  structural  configurations  is  often  fully  described 
by  a  set  of  linear  algebraic  equations.  In  general,  when  these  linear  equa¬ 
tions  are  put  in  matrix  form,  they  can  be  partitioned  as  shown  in  Fig.  B-l. 

The  blank  spaces  in  the  coefficient  matrix  are  zeroes  and  [c^l,  [b^]  and 
[A  3  are  matrices  of  orders  m^by  m^  by  m^  and  m^  by  mi+j  respectively 

Z .is  the  vector  of  unknowns,  each  of  order  m^  by  one  and  there  are  N  such  vec¬ 
tors.  Let  Z  be  the  common  unknown  vector.  Moreover,  g  is  also  a  vector  of 

L  i 

order  m^  by  one  and  d^  is  a  vector  or  order  m^  by  one,  which  includes  the  co¬ 
efficients  of  the  common  unknown. 

Note  that  the  presence  of  vectors  d^  make  the  whole  coefficient  matrix 
nonbanding  and  irregular.  If,  on  the  other  hand,  the  d^ -vectors  do  not  exist 
then  the  coefficient  matrix  is  identical  to  that  of  Ref  c-1.  In  this  case, 
the  matrix  is  a  banded  tridiagonal  matrix  with  zeroes  everywhere  and  with, 
at  most,  three  submatrices  banded  along  the  diagonal  as  shown  on  Fig.  c*l. 
Therefore,  the  present  case  la  a  bit  more  general  than  that  of  Ref  c-1.  The 
solution  procedure,  though,  is  basically  the  same  on  that  of  Ref.  c-1. 

C.l  Description  of  the  Algorithm 

The  explicit  form  of  the  system  of  linear  equations  of  Fig.  C.l  is  given 
by 

[bj{&] +[ao(z]  +  {4}  **<<)«  (gj 


[Gita  [Bjlftl  +  C  41 1  +  W  M)m[% ) 

(.  J  I  J  l  J  L  J 


W iik  L~  2.3, 


N-t 


CCi){£i-J  +CBi][?i  +  [4][4.}  =[<? i] 


viiVi  i-  L  - 1  ,  L  ,  L+  I 


CCj[£«.,J+f8„][^)+[^)2t(i)={^j 


(c-o 


Note  that  Z^Q)  is  one  element  of  the  common  unknown  vector  (see  Fig  C. l). 
A  short  description  of  the  solution  procedure  is  next  outlined. 

By  using  Gaussian  elimination  for  the  first  (L-2)  matrix  equations,  one 
may  find  the  equivalent  set  of  equations,  which  is 


(2i) +  CHlfiJ  i{eJ  z,u)=  (x,J 


where 


[p,]  -  CB]  C/1.]  s  {Ej-CBj'W 

CCBi]-cci)c?i.,)]'CA] 


(C-J) 


(C-3) 
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{E;j  =  CCBil-CCOCP,.,)] 

{XiJ  =  fCBJ-fCOCp^rfe-Cc,)^..]] 

Jn  i*2,i.  .1-1 

Note  that  the  order  of  the  various  matrices  is  as  follows: 


CCil 

Ay  Mi., 

fej 

Mi  4f  Mi 

(41 

nu, 

(Pi) 

Mi  %  Mu, 

Ujh  {g^,  UjL  {)<£}  and  {E^  are  all  by  1 

Next,  for  i  =  L-l,  L,  and  L4-1  the  equivalent  equations  are: 


{2k}  4  «  fXi„] 

-fa  i*l-l,L,L+l 

where,  for  i  =  L-l 

(Pj  =  [CBJ  -CCJCfajj'fa,]- CC,]  [ 


{Xij  -  CUM  -  ICO  cm"  ( f?ij  -  £C,]  ffc.} ) 


C  Ei-»]  =  [  0  [Enj  0  ] 


( C-7J 


!  _  '1 


Note  that  is  an  m^_^  by  m^+^  matrix  (defined,  as  shown,  for  convenience) 


and  for  1  =  L,  L  +  1 


(PJ  -  fCBi] -  ccj  raj]  (Ai) 


{XiJ  -ccdrJuKfaj-fcjfoj) 


(C-g) 


Finally,  for  1  -  L  +  2,  L  +  3 . K,  before  writing  the  equivalent 

j'  ^ 

equations,  ydj  Is  eliminated  from  each  matrix  equation.  The  elimination 
Is  accomplished  by  multiplying  jd^  with  the  appropriate  terms  of  matrix 

'  This  leads  to  a  matrix  with  only  one  nonzero  column  (vector) ,  as  shown 


below 


f fdsj  •  CPjJ  - 


diW?i(L.v 


diWP^i.xo. 


fC-9) 


Note  that  the  symbol  ©  is  Introduced  to  define  the  operation  that  leads  to 
the  matrix  of  Eq  (C*9). 

Similarly,  the  symbol  O  is  introduced  to  define  an  operation  that  leads 
to  a  column  matrix. 


rv.ci)  vxo)\ 
V,(i)  \AU) 


(C-IO) 


With  these  definitions  one  may  now  write  the  equivalent  equations  for 
i  =  Irf  2,  l/f3,  .  ...N-l.  These  are 

{^}  +  CRl{ii«}  *  {x‘j  (c-ii) 

where 


CPJ  =  (C BJ  -  Cdd CB-3] ' 04d  (c-«) 

{Xi}  =  [CBiJ  -  CCiUPjfttfr]  -  (Cj{x,.,j )  (C-'SJ 

with 

CCJ  =  (Ci)  -  CPJfc,,] . OiJ  (C-/4) 
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and 


ffl  --  W  -  M  ©  ( W  -  m]  (x,,}  +  c  pj  ctu  to«| 


Finally,  for  i  =  N 

(£*}  =  (x*j  (C-/^) 

where  '^X  '  is  given  by  Eq  (C-13)  with  i=N.  The  recurrence  formulae  for  back¬ 
ward  substitution,  in  order  to  calculate  Z^_2»***  Z2’  ant*  zi  are 

I4.j  --  fX„i 

[?ij  =  {Xij  -  tPi]  ;  UH-I.N-1.  ~  -  -l-l 
{?ii  r  [Xi]  “  CPJ  (£,»,]  'ff,j  (*)  ,  t=  L-i.i-s.  —  .i.  i  (C 

G2  Deternininant  Calculation 


(C-/7J 


In  each  step  of  the  inversion  process,  one  must  calculate  the  corresponding 
determinant  e^,  namely 


42,  =  M  (8,] 


=  MftBiJ-CcjCP,  ;)];  Ul.3, . . 


«i*  ;l -L+*  ,L+  3,---  /V 
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V 


Thus  the  determinant,  D,  of  the  entire  coefficient  matrix  of  the  system  can 
easily  be  computed  by 

^>-  tt  ( c-n ) 

i-  I 
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Appendix  d 


INSTABILITY  OF  LAMINATED  CYLINDERS  IN  TORSION 


by 

D.  Shawt  and  G.  J.  Simitses+t 
School  of  Engineering  Science  and  Mechanics 
Georgia  Institute  of  Technology,  Atlanta,  Georgia 


Introduction 

A  Galerkin-type  solution,  for  the  buckling  analysis  of  a  perfect 
geometry,  laminated,  circular,  cylindrical  thin  shell  subjected  to  pure 
torsion,  is  presented.  The  torsion  is  applied  through  the  reference 
surface,  which  is  the  midsurface  of  the  laminate  and  the  boundaries 
are  classical  simple  supports  (SS-3).  The  analysis  is  based  on  Donnell* 
type  nonlinear  kinematic  relations  and  linearly  elastic  material  behavior. 
It  is  assumed  that  a  primary  state  exists  and  that  it  is  axlsymmetric . 

This  primary  state  can  be  obtained  by  solving  the  field  equations.  Through 
perturbation  of  the  governing  field  equation  a  set  of  (linearized)  buckling 
equations  is  obtained,  along  with  the  related  boundary  conditions.  A 
Galerkin  procedure  is  employed  for  solving  the  buckling  equations.  Thus, 
the  problem  is  reduced  to  an  eigen-boundary-value  problem.  Critical 
torsional  loads  are  obtained  for  several  Boron/Epoxy  configurations  of 
symmetric,  antisymmetric  and  asymnetric  stacking.  Tn  addition,  approxi¬ 
mate  buckling  modes  are  established  for  both  positive  and  negative  torsion. 
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Governing  Equations  and  Solution  Procedure 


The  geometry  and  sign  convention  are  shown  on  Fig.  1.  The  torsion 
is  positive  if  applied  clockwise  at  the  right  end  (x  =  L)  and  counter¬ 
clockwise  at  the  left  end  (x  =  0) .  The  governing  equations  for  a  general 
laminated  circular  cylindrical  shell,  with  or  without  orthogonal  stiffen¬ 
ers,  without  geometric  imperfections,  and  subjected  to  a  pure  torsion,  con¬ 
sist  of  two  coupled  partial  nonlinear  differential  equations  in  the 
transverse  displacement  component  w(x,y)  and  an  Airy  stress  (resultant) 
function,  F(x,y).  One  of  the  equations  characterizes  transverse  equilibrium 
and  the  other  in-plane  compatibility.  These  equations  are  taken  from  [D.l] 
by  setting  N  =  q  =  w°(x,y)  =  0,  where  N  denotes  the  uniform  axial 
compression,  q  lateral  pressure  and  w°(x,y)  an  initial  geometric  imperfec¬ 
tion.  The  two  equations  are 
Equilibrium: 


Wi,r  ,  I  -u_  ,  r  ,  T  u,  _  w .  -r  u,_w,  *r  *.u,  0w  - 

11  yyxx  21  ’xxxx  31  ’xxxy  11  ’xxxx  12  xxzz  13  xxxy 

+  2b.,F,  +  2b0_F,  +  2  b__  F ,  +  2d„w,  +  2d..w,  +  4d,,w, 

13  ’xyyy  23  'xxyy  33  xxyy  31  xxxy  32  ’xyyy  33  xxy 

+  b..F,  +  b„„F,  -  b00F,  +  d_.  w,  +  d__  w,  +  2d  ,w, 

12  yyyy  22  xxyy  32  xyyy  21  xxyy  22  yyyy  23  xyyy 

+  |  F,  +  F,  w,  +  2N  w,  -  2  F,  w,  +  F,  w,  -  0  f  P- 

R  xx  yy  xx  xy  xy  'xy  xy  xx  yy 


Compatibility: 


a.,F,  +  a,_F,  -  a,«F,  +  b,,w,  +  b,»w 

11  ’yyyy  12  ’xxyy  13  ’xyyy  11  ’xxyy  12 

+  a,„F,  +  a„_F,  -  a00F,  +  b„.W,  +  b0-w 

12  ’xxyy  22  ’xxxx  23  ’xxxy  21  ’xxxx  22 

-  a,  ,F,  -  a„-F,  +  a.-F  -  b.,w,  -  b-„w, 

13  xyyy  23  'xxxy  33  xxyy  31  xxxy  32 

w, 

xx  . 

-  — r—  +  w,  w,  w,  w, 

R  xy  xy  xx  yy 


,  +  2b. _w, 

yyyy  13  ’xyyy 

,  +  2  b..  w, 

xxyy  23  xxxy 

-  2  b,_  w,  1 
xyyy  33  'xxyy 


(D-2) 
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where 


-  [A^r1  ;  [byl  -  [A^r1  [Btj] 

[dij]  =  [b^]  -  [d^] 


and  [a^],  [b^I  an{*  ^Dij3  are  the  extensions!.,  coupling  and  flexural 
stiffnesses  appearing  in  the  usual  lamination  theory. 

The  expressions  for  the  simply  supported  boundary  conditions  (SS  -  3) 
are  given  below  in  terms  of  w  and  F  (at  x  *  0,  L)  . 


b  -F,  +  d  w,  +  2d,_w,  -  b-.F,  =  -  b_,N  ; 

21  xx  11  ’xx  13  xy  31  *xy  31  xy 


a00F,  -  a_0F,  +  b0,w,  +  2b,,w,  =  -  a__N 

22  xx  23  xy  21  xx  23  xy  23  xy 


(D-4) 


where  N 

xy 


is  the  applied  torsional  stress  resultant. 


[C.l]. 


For  more  details  see 


It  is  assumed  that,  under  the  action  of  pure  torsion,  a  primary  state 
exists,  which  is  axisymmetric  (all  three  reference  surface  displacment 
components,  u,  v  and  w,  are  independent  of  the  circumferential  coordinate  y) . 
Note  that  for  symmetric  construction  (regular  angle-ply  or  cross-ply  with 
odd  number  of  plies,  for  example)  a  membrane  state  exists  and,  therefore, 
the  above  is  not  an  assumption.  How  reasonable  this  assumption  is  depends 
on  the  nature  and  magnitude  of  the  coupling  stiffnesses  [b^].  Primary 
state  quantities  are  denoted  by  tilda.  With  this  assumption,  the  field 
equation  becomes 


b-.F, 

21  xxxx 

A 

*22F’xxxx 


+  d,,w,  +  F,  /R 

11  xxxx  ’xx 


+  b-,v,  +  w,  /R 

21  ’xxxx  xx 


0 

0 


(D-5) 
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Moreover,  the  expression  for  the  reference  surface  hoop  strain  c° 

yy 

is  given  by 


*yy  '  - 


a0(,F,  +  a0_  N  +  b.,  v,  f 

22  ’xx  23  xy  21  xx  J 


(D-7) 


These  three  equations,  Eqs.  d-5,  D-5  and  d-7»  are  employed  to 

A  ' 

eliminate  F  and  thus  there  is  only  a  single  field  equation.  This 
resulting  equation  is: 


(A  b  21 N  -  .  1 

Id-.-  Z - )  w»  +  2  — 

\  11  a_„  /  xxxx  a. 


„R  ’xx 


23  — 
~  N 


..„R2  ‘22R  Xy 


(D-8) 


^  A 

The  general  solutions  for  w  and  consequently  [from  Eq.  D-7J  for  F, 


become 


v  «  sinh  Xj(x  -  j)  sin  >2(x  - 

+  B2  cosh  Xj(x  -  cos  X2(x  -  -R  a23 

-  7^  (b2lV*l  -  X2X  +  2  b2lVlX2  + 


(D-9) 


°2  L  l 

+  — )  cosh  X3(x  -  p  cos  X2(x  -  p 

'  <b21Bi<Xl2  -  X2>  *  2b21B2XlX2  + 


B1  1  1 
+  — )  sinh  Xj(x  -  p  sin  X2(x  -  p 


(D-10) 
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-Ufa  R2,l21_d  .-T*  ,  1  f°21  .  f.  D2l'l\ 

1  b  ^a22R  <a22  dll>;i  +  2(i^R><dll-  7£  1 


%  ^  % 
*22 


_  i I  r.  b2/^21  .  lf  b21v, .  b21 

-  \2  La22R  (a^  -  dn)]  -  2(a22R)(dll  -  a22)  j 


(P  -ID 


The  constants  Bj  and  B2  can  be  obtained  by  making  use  of  the 
boundary  conditions,  Eqs. 

Next,  the  buckling  equations  are  obtained  through  a  perturbation 
of  the  nonlinear  governing  equations.  The  dependent  variables,  w  and  F, 

A  A 

are  replaced  by  the  sum  of  the  primary  state  parameters,  w  and  F,  and 
small  additional  quantities,  w*  and  f\  necessary  to  represent  the  buckled 
state.  Moreover,  the  related  boundary  conditions  for  the  buckling  equations 
are  also  obtained  in  the  same  manner.  Note  that  since  the  additional 
quantities  can  be  made  small  as  one  wishes,  only  the  linear  terms  in 
w*  and  F*  are  retained. 

The  buckling  equations  and  related  boundary  conditions  are: 


b21F’xxxx*  (2  b23~b31)F,xxxy+  <bll"  2  b33+  b22)F’xxyy+  (2  b13’  b32)F’xyyy 

+  b12F’yyyy+  dllW,xxxx+  (2  d31+ 2  d13)  W,xxxy+  (d12+  4d33+  d21)  W,xxyy 
+  (2d32+2  d23)  w,xyyy+  d22wlyyyy+  -R—  +  F,^  w,yy 


+  w,  fJ  +  2  N  wl  *  0 
’xx  ’yy  xy  xy 


(D-12) 


a22F,xxxx”  2a23F,xxxy  +  ^2  al2+  a33^F,xxyy"  2  a13F,xyyy+ allF,yyyy 


+  b21W,xxxx+  (2  b23-b31)w,xxxy+  <bll"  2  b33+  V  W,xxyy 
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1  1  *  xx  A  1 

+  (2b.,-b,,)w,  +  b,  0w;  +  — r—  +  w,  w,  =0 

13  32  xyyy  12  yyyy  R  xx  ’yy 

w  =  0  ;  b^F*  +  b,,F*  +  d,  ,w*  +  2d, .w*  =  0  ; 

at  21  ’xx  31  xy  11  *xx  13  xy  * 

x  =  0,L 

F’yy  ’  a22F’xx+  a23F’xy+  b12W*xx+  2b23W’xy  =  0  ' 


•  D-13) 


<b-14) 


The  Galerkln  procedure  is  employed  for  both  equations.  The  following 
approximate  series  is  used  for  generating  the  Galerkin  integrals.  Note 
that  the  boundary  conditions  are  satisfied  by  each  term  in  the  series. 

N  M 

1,1  ‘  Ij  “  +  BinslnT?>  ,inir  •  (life? 

pl  -  Ij  +  v1”?1  it  ,lrnr  ■  (i«)S  *lrrU±P2iJ 

Substitution  of  the  above  expressions,  Eqs.  d.15>  into  the  buckling 

equations  results  into  a  set  of  systems  of  linear  homogeneous  algebraic 

equations  in  A,  ,  B.  ,  C.  and  D.  for  each  n  (decoupled  with  respect  to  n) 
in  m  xn  in 

Assuming  that  the  lowest  eigenvalue  corresponds  to  the  critical  load, 

N  ,  a  computer  program  has  been  written  to  this  effect.  The  Georgia 
^cr 

Tech  high  speed  digital  computer  CDC  -  CYBER  -  170/760  is  used  for 
generating  data.  Note  that  a  minimization  with  respect  to  n  is  per¬ 
formed  in  order  to  find  the  lowest  eigenvalue. 

Numerical  Results  and  Conclusions 


The  geometries  considered  in  the  Investigation  represent  variations 
of  the  one  report  in  d.2.  Each  lamina  is  orthotropic  (Boron/Epoxy; 

AVCO  5505)  with  the  following  properties: 
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Eu  -  2.0690  x  108  kN/m2  (30  x  106psi.)  ;  p  =  0.21  ; 


E„  -  0.1862  x  108  kN/m2  (2.7  x  106psi.);  R  *  190.5  cm(7.5in.); 


G,,  -  0.04482  x  108  kN/m2  (0.65  x  106  psi)  ;  L=  381  cm  (15  in.); 


h  ,  =  0.013462  cm.  (0.0053  in.) 

ply 


<D-16) 


(hply  M  \  ~  \-l  for  k  °  1»2*3’4  •  four  Plies> 

Five  different  stacking  combinations  of  the  four-ply  laminate 

r 

comprise  the  various  geometries,  I-i,  i  =  1,2,  --5.  These  are 


I  -  1 


450/-45°/-450/45° 


1-2 


:  45°/ -45°/45°/ -45° 


1-3  :  -45°/45°/ -45°/45° 


1-4  :  90°/60O/30O/0° 


(1>— 17) 


1-5 


:  0°/ 30°/ 60°/ 90° 


where  the  first  number  denotes  the  orientation  of  the  fibers  of  the 

outermost  ply  with  respect  to  x,  and  the  last  of  the  innermost.  A  pure 

torsion  is  applied  through  the  midsurface  of  the  four-ply  laminate. 

Some  of  the  generated  results  are  shown  on  Table  D.l.For  each 

geometry,  the  critical  torsion  (for  both  positive  and  negative 

application;  clockwise  and  counterclockwise  at  the  end  x  =  L) ,  the 

minimizing  value  of  n  (full  number  of  circumferential  waves),  and  the 

values  of  the  coefficients  A^n  and  B^n  (normalized  with  respect  to  B^n) 

are  shown.  Note  that  the  A,  and  B.  when  substituted  into  the  first  of 

m  in 

Eqs.  D. 55  yields  the  buckling  mode.  It  was  concluded  that  M  =  5  suffices 


for  determining  critical  loads. 
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Table  d-1.  Numerical  Results 


Note  that  Geometry  I  -  1  is  symnetrlc  (with  respect  to  the 
midsurface) ,  Geometries  1-2  and  1-3  antisymmetric,  and  Geometries  1-4 
and  1-5  asymmetric.  For  the  symmetric  geometry  (1-1),  the  positive 
direction  critical  torsion  is  6987  N/ln  (39.9  lbs. /in.),  while  the 
negative  critical  torsion  is  13,220  N/m  (75.5  lbs. /in.).  The 
respective  reported  D-2  experimental  values  are  4640  N/m  (26.5  lbs. /in.) 
for  the  positive  direction  and  11,508  N/m  (65.72  lbs. /in.)  for  the 
negative.  This  suggest  that  the  geometric  imperfection  in  the  tea teg 
cylinder  D-2  is  such  that  the  configuration  Is  more  sensitive  to  it,  when 
loaded  in  the  positive  direction,  than  in  the  negative  (the  ratio  of  the 
experimental  to  theoretical  value  is  0.664  for  the  former  and  0.87  for 
the  latter).  The  difference  in  response  is  understandable,  because  of 
the  anisotropy.  The  antisymmetric  geometries,  1-2  and  1-3,  yield  the 
same  response  when  loaded  opposite  to  each  other.  Note  that  the  positive 
direction  critical  load  for  1-2  is  the  same  as  the  negative  direction 
critical  load  for  1-3  (the  same  is  true  for  the  buckling  mode).  Also, 
observe  that  the  two  (+  direction)  critical  loads  are  very  close  (9534 
N/m.  and  9454  N/m.).  This  is  due  to  the  fact  that  the  extensional,  [a^], 
and  flexural,  [d^U,  stiffness  have  the  same  form  as  if  the  shell  were 
isotropic.  The  difference  from  isotropy  is  the  existence  of  some  small 
(in  value)  terms  in  the  coupling,  [b^I,  stiffnesses. 

Finally,  for  the  asymmetric  configurations,  1-4  and  1-5  the  response 
is  completely  different  when  each  geometry  is  loaded  in  the  positive  and 
in  the  negative  direction.  Although  the  [Ajjl  and  [Djj]  stiffnesses,  for 
the  two  configurations,  are  the  same  and  only  the  signs  are  different 
in  the  [b^]  stiffness,  the  geometries  behave  (radically)  differently. 

The  only  similarity  is  that  the  number  of  full  waves,  n,  is  approximately 

the  same  (12  and  13). 
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